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!  .  l.N'l'KODUC'J'lU.N 

The  conference  on  "ilF  Radio  Propagation"  was  held  at  the  Aiierton 
House  Con fe  fence  Center,  MonLicello,  Illinois,  on  June  2,  J,  and  A,  l'^u1,/. 
ihe  con 1 rente  was  principally  concerned  with  the  propagation  of  Hi- 
radio  signals  and  the  effect  ol  the  characteristics  of  tlie  propagation 
on  RDF  systems.  Planning  and  operation  of  the  conference  were  under¬ 
taken  by  the  Radiolocation  Research  Laboratory  of  the  Department  of 
Electrical  Engineering,  University  of  Illinois  at  Urbana/Champaign .  ihe 
Conference  was  sponsored  by  the  Office  of  Naval  Research. 

The  objectives  for  the  conference  may  be  summarized  briefly  as 
follows : 

1.  To  assist  in  delineating  those  areas  of  investigation  and  devel¬ 
opment  which  should  be  emphasized  for  the  future,  if  significant 
improvements  in  the  operation  of  RDF  Systems  are  to  be  expected. 

2.  To  examine  the  present  status  of  RDF  Systems. 

3.  To  examine  the  limitations  of  these  RDF  Systems  which  are  a 
consequence  of  propagation  characteristics. 

4.  To  examine  the  understanding  of  the  ionosphere  and  HF  radio 
propagation  which,  it  is  believed,  can  contribute  to  a  reduction 
of  the  limitations  noted  in  (2)  above. 

Attendance  at  the  conference  was  by  invitation  and  included  repre¬ 
sentatives  from  several  agencies  of  the  Department  of  Defense,  ESSA, 
research  organizations  conducting  investigations  in  the  area  included 
in  the  conference,  and  the  University  of  Illinois. 


The  agenda  of  presentations  and  discussion  are  given  on  the  pages 
following  tills  introduction.  Those  presentations  for  which  papers  are 
included  in  this  Proceedings  are  indicated  by  two  asterisks  (** )  pre¬ 
ceding  the  title  of  the  paper. 
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AGENDA 

CONFERENCE  ON 
II F  RADIO  PROPAGATION 

June  1! ,  '3,  4,  19G9 

Allerlon  House  -  Mouticello,  Illinois 

MONDAY,  JUNE  2 

BREAKFAST  8:00  -  9:00  a-m. 

CONFERENCE  OPENING  9:00  a.m. 


Introduction:  Survey  of  RDF  and  Propagation  Research 

at  the  University  of  Illinois.  Chairman,  E.  W.  Ernst. 

"Early  History  of  RDF  Efforts  at  the  University 
of  Illinois."  E.  C.  Jordan. 

**"Recent  and  Current  Status  of  RDF  of  the  University 
of  Illinois."  A.  D.  Bailey. 

COFFEE  BREAK  10:30  -  10:45  a.m. 


SESSION  1. 
9:00-10:30  a.m. 


SESSION  II.  Direction  of  Arrival  Studies  over  a  Medium  Range 

Path  (452  KM).  Chairman,  A.  D.  Bailey. 

10:45  a.m.  - 

12:15  p.m.  **"HF  Direction  of  Arrival  Studies  over  a  Medium 

Range  Path  (452  KM)."  A.  D.  Bailey. 

**"Short  and  Medium  Range  Results  with  a  DF  Inter¬ 
ferometer."  P.  E.  Martin,  R.  B.  Mathews,  C.  Dodge. 

LUNCHEON  12:15  -  1:15  p.m. 


SESSION. III. 
1:30-3:00  p.m. 


Direction  of  Arrival  Studies  on  a  Long  Range  Path. 
Chairman,  A.  D,  Bailey. 

**"Direction  of  Arrival  Studies  on  a  Long  Path  (Sum¬ 
mary)."  A.  D.  Bailey. 

*A"The  Coordinated  HF  Propagation  Experiment."  L.  H. 
Tveten . 


COFFEE  BREAK 


3:00  -  3:20  p.m. 


t) 


SESSION  IV. 
j  -  _  0  -  E)  p.m. 

SOCIAL  HOUR 
DINNER 

VLSL'L'  TO  rUL  U .  OF 

TUESDAY,  JUNE  3 
BREAKFAST 

SESSLON  V. 
8:30-10:00  a.tn. 

COFFEE  BREAK 

SESSION  VI. 
10:30-12:00  noon 


LUNCHEON 


Propugut  lull  '.indie:.  Us  i  up,  Direction  oi  Arrival  Data. 
Cliatrman,  N.  N.  Ran. 

**"Ray  Tt. icing  Simulation  ol  Ionospheric  Effects  in 
High  Frequency  Radiolocation  "  N.  N.  Rau . 

3:00  -  b : 00  p.m. 

b : 00  -  7:00  p.m. 

t.  FIELD  STATIONS  OR  OTHER  ACTIVITIES  7:00  -  10:00  p.m. 


7:30  -  8 :  30  a .  m 

RDF  Research  by  University  of  Ij.lj.noas  Graduate 
Students.  Chairman,  W.  R.  Flowers. 

"Auroral  Backscatter  Studies."  F.  Fahlstng 

**"Short  Range  High  Angle  Studies-"  E.  A.  Walton 

**"Aperture  Size  Effects  in  a  HF  Interferometer 
RDF  System  "  C.  R-  Talbott - 

"Tracking  Mode,  CDAA  RDF  System-"  H.  M.  Rice. 
10:00  -  10:30  a-  m- 


RDF  Propagation  Research.  Chairman,  A.  D.  Bailey. 

**" Interferometer  Instrumentation  for  Propagation 
Research."  D.  N.  Travers,  W-  M  Sherill,  J.  D. 
Moore  - 

**"The  Use  of  Reflectrix  Techniques  in  the  Analysis 
of  FM/CW  Oblique  Ionograms,"  M-  E.  Coffey,  R.  P. 
McConville  (presented  by  S.  M.  Bennett). 

**"Review  of  Small  Aperture  Investigations  at  Georgia 
Tech."  H.  Jenkins,  R-  W.  Moss- 


"Automatic  Radio  Phasemeter 
Frequency  Direction  Finder 

12:00  -  1:00  p-m. 


Technique  and  High 
Concepts."  R-  Russell 


TREASURE  HUNT 


12:00  -  2:00  p-m- 
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SESSION  VI 1 . 
2:00-3:30  p.m. 

COFFEE  BREAK 

SESSLON  VLLL. 
3:45-5:15  p.m. 


SOCIAL  HOUR 
DINNER 

VISIT  TO  U.  OF  I. 

WEDNESDAY,  JUNE  ^ 
BREAKFAST 

SESSION  IX. 
8:30-10:00  a.m. 
COFFEE  BREAK 

SESSION  X. 

10:30-12:00 

LUNCHEON 


Digital  Mu i ho ilb  tar  RDF.  Chairman,  I  .  N  ■  Ernst 

**"Digilal  MeLhuds  i<>r  Radi'.'  Direction  Finding. 

E.  W.  Ernst 

**"Digi Lai  DF  Applications."  W.  M.  Sherrill,  1- 
Green,  R.  1 orenz. 

3:30  -  3:45  p.m. 


RDF  Systems  and  Instrumentation  lor  Propagation 
Research  Chairman,  E.  W-  Ernst. 

**"Elevation  Angle  Measurements  on  the  University  of 
Illinois  Wul lenwebe r  Antenna  Array."  L.  J.  Miller. 

**"0blique  Incidence  Lonosonde  Amplitude  Data. 

W.  W.  Wood. 

**"Wave front  Testing  and  its  Applications."  G.  A.  Smith 

"Azimuthal  Noise  Measurements."  W.  Flowers. 

"System  Testing  and  Maintenance  Procedures.  W.  W. 
Wood. 

5:15  -  6:00  p.m. 

6:00  -  7:00  p.m. 

FIELD  STATIONS  OR  OTHER  ACTIVITIES  7:00  -  10:00  p.m. 


7: 30  -  8:30  a.m. 

The  Shape  of  the  Next  Generation  of  RDF. 

10:00  -  10:30  a.m. 

Future  Directions  for  RDF  Research, 

12:00  -  1:00  p.m. 


CONFERENCE  ADJOURNMENT  1:00  p.m. 


preceding  page  blank^-  not  filmed 
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1 1 

A.  D.  Bailey 

-  i’i ’  i j*r c 
ivju:': 

RECENT  AND  CURRENT'  STATUS  OF  RDF  OF  THE  UNIVERSITY  OF  ILLINOIS 

Hi  is  is  the  third  conference  on  RDF  to  be  held  at  the  Allerton  House. 

The  first  conference  was  held  2-3-4  June  1959,  exactly  ten  years  ago.  The 

second  was  held  26-29  May  1964,  five  years  ago. 

1  hope  this  can  be  continued  and  quite  possibly  we  will  be  meeting 
together  again  on  or  about  3-5  June  1974.  Dr,  Jordan  opened  each  of  the 
previous  conferences  and  he  is  hereby  respectfully  requested  to  start  plan 
ning  his  opening  remarks  for  the  next  one  in  five  years. 

Our  emphasis  at  this  meeting  is  on  the  use  of  the  RDF  systems  that  we 

have  developed  at  the  University  of  Illinois  in  the  HF  Directional  Propa¬ 

gation  Research.  Such  studies  are  of  academic  and  scientific  interest 
to  our  University  and  its  Department  of  Electrical  Engineering.  The 
results  of  these  studies  are  often  of  interest  to  the  sponsors  of  the 
program  in  RDF  at  the  University.  We  hope  that  this  arrangement  may  con¬ 
tinue.  The  research  efforts  are  currently  sponsored  by  ONR  and  USAECOM. 
There  are  two  or  three  other  agencies  within  the  DOD  who  supplement  the 
basic  ONR  contributions. 

The  previous  conferences  were  concerned  primarily  with  the  work  done 
at  the  University.  In  this  conference  we  will  hear  several  papers  pre¬ 
pared  by  other  research  groups  with  whom  we  have  cooperated  in  one  way 
or  another.  This  is  as  it  should  be;  in  the  end  it  has  to  be  some  kind 
of  a  team  effort.  Any  directional  propagation  experiment  has  to  have  at 
least  two  participants,  one  at  each  end  of  the  path.  Many  of  the  tasks 


1J 


appeal  to  us  to  be  as  enormous  as  "seeing  to  the  lop  ul  I’aul  Bunyun's 
upside  down  mountain,"  and,  as  in  that  ease,  one  needs  cooperation. 

There  will  be  two  visiting  periods  to  the  field  stations.  We  will 
not  have  much  to  say  about  Lhe  details  of  the  system  per  se ■  Many  of  your 
technical  questions  on  our  systems  can  be  answered  at  the  time  of  the  field 
station  visits.  Tonight  the  L'DAA  and  the  Interferometer  sites  are  open  for 
"knob  twiddling"-type  inspections.  Tomorrow  night  we  will  have  a  propa¬ 
gation  experiment  going  and  one  can  hopefully  see  real  live  system  opera¬ 
tion  and,  if  he  wishes,  make  a  comparison  of  the  two  systems. 

Our  program  schedule  outlines  pretty  well  what  we  hope  to  cover. 

I'm  particularly  interested  in  what  other  speakers  have  to  say  so  I'll 
get  on  to  what  I'm  scheduled  to  say  and  let  the  others  have  their  chance. 
Accordingly,  I'd  like  to  begin  talking  about  the  second  session  material 
now  and  get  as  much  of  this  covered  as  possible  between  now  and  the 


coffee  break. 


A.  L>.  BALLHY 


..  c  :  ‘ l  ”.l  1  z .  t.  oi  i‘  i  ou  -  !' f i  t'  •  c .  t  t* i  n  j 

i>: .  :v  nt :  :  >.v:  a,  1  /<•  ,(/;  i,  no:  J 

11 F  D l RUCTION  01-  ARRIVAL  STUU1KS  OVKK  A  MLD1 UM  KAN  UK  PATH  (4!>2  KM) 


1.  PHILOSOPHY 


:  jajl.v.;;  jl'&c'j,  ii'tir,  uu-wokeeiu  \  .t  a;:l  ihn:::^r. 

Complementary  to  t lie  statement  of  the  four  general  objectives  wiiicl 
was  carried  in  the  letter  of  announcement,  1  would  like  to  add  a  state¬ 
ment  to  the  effect  that,  "There  is  a  time  for  rethinking  and  that  time 
is  now!"  i  hope  that  each  of  us  has  come  to  this  conference  in  this 
spirit.  Rethinking  is  so  necessary  at  this  time  when  Universities  and 
Institutes  for  Research  are  looking  for  ways  and  means  to  continue  sup¬ 
port  of  our  continuing  search  for  the  truth  and  the  advancement  of  the 

« 

state  of  the  art  and  science  in  our  disciplines. 

Before  taking  up  the  subject  of  this  particular  session,  it  is 
necessary  that  we  review  a  few  definitions  and  concepts  which  will  pro¬ 
vide  a  basis  for  arguments  and  later  decisions  which  we  might  reach. 

Certainly  each  of  us  has  at  some  time  or  another  thought  about 
the  "two-point  problem."  Thus,  given  two  separated  points  A  and  B,  or 
T  and  R,  or  X  and  Y,  etc.,  one  might  ask  the  following: 

How  t  ■  get  from  A  to  B. 

How  to  talk  between  A  and  B. 

How  to  deliver  something  from  A  to  B. 

How  to  exploit  B  from  A. 

How  to  destroy  B  from  A. 


How  to  deny  to  B  certain  benefits  at  A. 


Ln  Figure  L  1  have  imbedded  puiiiLs  A  and  li  in  a  coordinate  system 
L'lie  displacement  vector  A 15  is  a  iunction  oi  Lhree  variables  such  as 
x,  y,  x  or  r,  c.  thus,  a  triplet,  oi  values  is  needed  to  define  Lhe 

position  of  B  relative  Lo  A.  liiese  could  be  Lhe  range  All  and  Lhe 
two  angLes  / 1\  and  respectively. 

W<‘  will  stress  that  in  KDF/RPL  disciplines  at  11F,  range  and  the 
two  angles  are  needed  for  completeness. 

U.  DEF1N  1 1  ION  UF  SOME  GENERAL  TERMS 

Turning  then  to  a  practical  aspect  of  Li F  directional  propagations 
studies,  we  next  review  some  definitions  as  follows: 

DEFINITIONS 

DIRECTION  FINDER  —  "A  direction  finder  is  a  radio  receiving  device 

which  permits  determination  of  the  line  of  travel  of  radio  waves 
as  received."  American  Standard  Definitions  of  Electrical  Terms. 
(AIEE) . 

SENSE  FINDER  —  "A  sense  finder  is  that  portion  of  a  direction  finder 

which  permits  determination  of  direction  without  one  hundred-eighty 
degree  ambiguity."  ASDET/AIEE. 

RADIO  DIRECTION  FINDER  —  "A  radio  direction  finder  may  be  defined  as 
a  device  for  determining  the  direction  of  arrival  of  radio  fre¬ 
quency  energy.  It  is  a  receiving  system  and  operates  on  the  energy 
that  it  extracts  from  the  passing  radio  waves."  -  VHF  Techniques, 
RRL  Staff  of  Harvard  University. 

NOTES:  (1)  The  word  bearing  does  not  appear, 

(2)  The  words  direction  of  arrival  as  received  do 


appear. 


Lo 


{  J)  S  u>_  1 1  systems  could  be  at  live  in  ihc  belike  of  the 

del  ini  t  ions,  although  at  first  reading  one  might  infer 
that  RDF  systems  are  passive. 

(.4)  Direction  of  arrival  has  two  components,  e.g.,  j  and 

y  • 

Figure  2  ls  one  representation  of  a  radio  direction  finder  system 
in  functional  block  iorm.  An  observer  is  necessary  to  complete  the  pro¬ 
cess.  it  may  be  a  human  being  or  a  computer.  in  any  event  the  human 
being  is  still  important  --  he  is  the  master  --  he  is  the  truth  seeker. 

We  next  explore  briefly  some  of  the  applications  of  the  device  we 
have  defined. 

A  radio  direction  finder  has  application  in  at  least  the  following 
fields : 

1.  Navigational  Aids. 

2.  Military  Surveillance  and  Countermeasures- 

3.  Location  of  Sources  of  Radio  Interference. 

4.  Directional  Radio  Wave  Propagation  Studies. 

5.  Remote  Environmental  Sensing,  eg.,  OTH  RADAR  at  HF. 

Next,  we  should  note  that  the  first  three  are  of  vital  concern  to 
"management."  But  this  latter  disciplinary  group  wants  the  bearing  of 
the  signal  and  not  much  else!  Hence,  we  now  need  more  definitions,  to 
wit : 

BEARING  —  "The  situation  or  direction  of  one  point  with  respect 
to  another  or  to  the  points  of  a  compass." 


IS 


GCb  --  I'lu*  groat  circle  bearing  ol  something  as  it  applies  lu  a 
spherical  eat  Lit  is  laciily  understood  Il's  what  the 

managers  want 

Our  sponsors  are  alien  very  practical  people  who  are  willing  to  sup- 
pan  direction  oi  arrival  studies  buL  they  really  want  bearings! 

In  this  connection  one  may  recall  the  classic  story  of  "The  blind  Men 
and  the  l-.lephants . "  l' he  meaning  of  the  word  elephant  had  different  sen¬ 
se  rv  interpretations  to  six  different  men  depending  upon  where  the  in¬ 
dividual  touched  the  elephant. 

Managers  are  not  blind;  of  course,  they  need  and  want  the  GGii: 
they  buy  devices  with  this  need  in  mind  in  order  to  measure  the  bearing, but 
the  very  best  that  the  device  can  ever  do  is  to  measure  an  apparent  angle 
of  arrival  as  received.  It  is  up  to  people  like  ourselves  to  resolve  the 
dilemma.  Also,  the  bearing  can  be  very  elusive  in  any  case  other  than 
controlled  propagation  experiments  and  cooperative  transmitters.  When 
one  thinks  it  through,  the  bearing  of  a  noncooperative  radio  source  can 
only  be  known  in  the  sense  of  being  most  probable. 

A  ship  at  sea,  an  airplane  in  flight,  an  automobile  heading  west 
on  Highway  66  are  examples  of  situations  where  the  exact  location  at  any 
instant  is  hard  to  come  by,  particularly  when  there  is  lack  of  cooperation. 

Of  course,  our  motivation  here  is  to  find  better  ways  for  inferring 
what  is  the  best  estimate  of  GC13  given  a  set  of  measurements  of  angle 
of  arrival  as  received.  Again,  the  point  of  all  of  this  was  to  emphasize 
the  need  for  more  careful  use  of  the  word  bearing. 


i llTilrtftTilViVm  I  II 
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Wo  need  still  more  definitions  und  1  nexL  need  to  lalk  about 
RAN  G  i'. . 

It  appears  necessary,  Lo  me,  to  delineate  certain  cnaracteristic 
ranges  in  radio  position  iinding  work.  These  will  he  "hail  parks"  as 
opposed  lo  the  hard,  lasL,  and  rigorous. 

For  reasons  Lliat  wiil  become  evideiiL,  it  appears  to  me  that  one 
should  delineate  at  least  three  ranges  termed,  SHORT,  1NTLKT1LD1ATL ,  and 
LONG.  I  propose  to  do  this  by  joint  use  of  Lhe  Earth-Ionosphere 
geometry  shown  in  Figure  3. 

First  off,  we  note  cerLain  given  information; 

next, we  add  two  lines  to  delineate  the  three  sectors; 

we  discuss  the  properties  of  each  sector; 

we  indicate  the  antenna  element  types  that  are  appropriate  to 
each  sector; 

the  surface-wave  sector  is  interesting  but  not  very  meaningful  in 
the  sense  of  RANGE. 

Finally,  we  stress  clustering  of  rays  at  extremely  high  or  low 
angles . 

A  significant  contribution  to  bearing  deviation  at  the  shorter 
ranges  is  the  effect  of  lateral  and/or  longitudinal  displacement  of  the 
apparent  reflection  point  in  tne  ionosphere.  For  example,  assume  that 
the  reflection  point  is  laterally  displaced  by  5  kilometers.  The  ap¬ 
parent  deviation  in  azimuthal  angle  of  arrival  which  occurs  as  a  func¬ 
tion  of  range  appears  in  Table  1. 
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Locul  ions  close  Lo  those  ol  Lhe  unknown.  Thus ,  we  use  additional  in¬ 
dependent  measurements  which  should  reduce  the  uncertainty.  We  close 
this  introductory  discussion  wiLh  Figure  A,  which  illustrates  Lhat  fur 
very  shorL  ground  ranges,  the  slanL  range  line  of  travel  can  be  long 
but  alinosL  constant  and  not  much  different,  less  than  5  ■  percent  from  the 
effective  height  shown  on  a  vertical  incidence  ionosonde. 

ill.  DESCRIPTION  UF  AW  INTERMEDIATE  RAN  (IE  EXPERIMENT 

The  geometry  for  the  452  kilometer  path  K-W  directional  propagation 
experiment  is  shown  in  Figure  5.  It  was  an  example  of  an  intermediate 
or  medium  range  path.  l’he  important  subsystem  components  were  the 
pulsed  transmitter  located  at  OSU  Columbus,  Ohio;  the  midpath  v.i. 
ionosonde  just  North  of  Anderson,  Indiana;  the  triple  interferometer 
receiving  site  at  the  Monticello  Road  Field  Station,  S.W.  of  Savoy, 
Illinois;  and  the  receiving  end  v.i.  ionosonde  at  the  Bondville  Road 
Field  Station  3.5  miles  N.W.  of  the  receiving  site. 

The  transmitted  signals  were  cosine  squared  100  p  sec  pulses  having 
a  repetition  rate  of  approximately  200  pps.  Because  of  the  pulsed  nature 
of  the  signals, the  several  gross  propagate  modes  could  be  identified, 
separated,  and  the  angle  of  arrival  of  each  such  mode  could  be  deter¬ 
mined  from  the  differential  phase  data  obtained  with  any  two  of  the 
three  interferometers. 

All  of  the  reduced  data  obtained  over  the  18-month  period  of  ob¬ 
servation  are  contained  in  an  aTLAS.^ 

1,  A  Study  of  High  Frequency  Directional  Propagation  Over  a  450  Kilometei 
East-to-West  Path,  R.F.  Donnelly,  D.C,  Detert,  A.D.  Bailey,  RRL  Pub¬ 
lication  No.  254,  November  1965. 


L  should  R-cjgiiizt.'  the  cunt i ibut ions  ot  several  people.  Messrs. 

Lyle  Rocke,  Chris  McClurg  and  Gerald  .Smith  de  veloped  and  operated  the 
systems.  Graduate  SLudents  Hubert  Crush,  Richard  Monnelly  and  Mavid 
Detert  did  the  data  reduction.  There  were  many  others  who  also  served. 

Figures  b,  7,  and  8  are  examples  oi  the  reduced  data.  In  each  1 igure 
one  sees  live-time  series  oi  data;  Lhe.se  are: 

a)  Azimuthal  angle  oi  arrival 

b)  Incidence  angle  of  arrival  as  measured  by  one  pair  of  inter¬ 
ferometers  (0-120°) 

c)  Incidence  angle  of  arrival  as  measured  by  a  second  pair  of 
interferometers  (0-240°) 

d)  Predicted  incidence  angle  of  arrival  obtained  from  the  midpath 
ionosonde  data 

e)  Predicted  incidence  angle  of  arrival  obtained  from  receiving 
endpoint  ionosonde  data, 

Thr.se  data  permit  one  to  make  assessments  of  the  performance  of  the 
system  as  a  radiolocation  device. 

IV.  SUMMARY  OF  THE  FINDINGS 


There  were  several  important  findings  from  these  data  that  may  be 
listed.  The  attempts  at  radio  range  finding  often  came  out  short.  Mr. 
Winkelmann  showed  that  by  the  use  of  simple  ray  tracings  using  the 
so-called  Booker-Seaton  method  and  taking  into  account  both  the  under¬ 
lying  E-region  and  the  F-regions  he  could  increase  the  range. 


interferometer  measurements 
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Figure  6.  Columbus  Propagation  Tests,  January  14,  1963,  17:00 
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Figure  7,  Columbus  Propagation  Tests 
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Ml.  D.'V  l  cl  Di-liil  i  n  Vi  ,'i  l  i  g.il  t-il  llu  list  oi  muUil.s  u  J  Lluvclling 
JisLiiL  bailees  la  explain  l  lie  variations  in  i  in.  i  Jem  e  alible  ol  aiiival 
lie  achieved  veiy  consistent  lesuits  wilii  his  models 

i.  D.  Bailey  showed  Inal  Llieie  was  a  Vt  l  y  sLruiig  eo  r  re  i  a  L  i  on  be- 
Iwet  11  slight  bul  lloLiceable  cnanges  in  llie  elieelive  iioiglils  oi  Llie 
tuiapaih  ionograms  and  llie  lateral  deviations  in  the  indicated  uzimul.ial 
angle  oi  arrival  l  rom  GC1>.  Ibis  result  made  a  strong  impression  in  the 
mind  ot  having  an  ionosunde  aL  the  midpoint  oi  Llie  paLh  in  any  cun- 
t  rol led  study . 

the  studies  pointed  ouL  that  Lhe  so-called  random  iaLeral  devia¬ 
tion  error  could  be  controlled  after  all.  There  are  ways  and  means  of 
geLting  a  handle  on  iL. 

Experience  with  the  system  proved  out  its  potential  promise  as  a 
versatile  radio  direction  finding/ radio  position  Locating  device, 

The  separation  of  the  modes  in  the  vertical  plane  was  very  evident 
although  the  azimuthal  separation  was  small.  This  means  that  the  wave 
interference  fading  effects  have  broader  interpretations  and  an  alter¬ 
nate  method  for  determining  range  now  appears  feasible. 

The  lateral  deviation  error  was  in  evidence  in  all  data  at  all 
times  of  the  experiment , 


V,  SIGNIFICANT  CONCLUSIONS 

The  significant  conclusion,  in  my  mind,  is  the  fact  that  the 
triple  interferometer  system  as  instrumented,  for  example,  at  the 
University  of  Illinois  has  proven  to  be  a  powerful  research  tool  as 
well  as  a  practical,  true  radio  direction  finder  It  is  a  wide- aperture 
system  measuring  differential  phase  over  a  large  base  line  and  it  deter¬ 
mines  both  angles  of  arrival.  Our  contribution  is  not  unique;  of  course 
several  other  investigators  have  used  the  interferometer  system  at  HF 


with  similar  success. 
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Of  special  interest  at  this  Lime  is  Llie  article  by  J.  W.  Wright  el 
ESSA  entitled  "Some  Current  Developments  in  Radio  SysLems  ior  Sounding 
ionospheric  Structure  and  Motions,"  Proc  IEEE  VS/,  N4  ,  p  481.  lie  dis¬ 
cusses  a  K1NES0NDE  which  permits  one  to  measure  additional  daLa  from 
a  v.i.  ionosonde.  it  appears  Lo  me  Lhat  Lhis  is  in  essence,  or  can  be 
considered  to  be,  an  application  of  interferometer  principles  and  appears 
to  be  one  partial  solution  to  an  urgent  problem  Lhat  we  are  attempting 
to  solve  here  at  the  radiolocation  research  laboratory. 
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data  .sets  are  dumped  to  paper  tape  (only  150  sets  can  now  be  stored  in 
the  51u's  8k  memory  due  to  programming  requirements ) ,  The  results  to 
be  described  will  be  from  off-line  processing.  *’  At  present,  it  requires 
about  100  milliseconds  to  acquire  and  process  a  frame  of  data;  that  is, 
the  four  baseline  phase  angles  to  compute  the  angle  of  arrival- 

The  digital  system  has  been  operated  on-line  for  approximately  one 
y.ear;  however,  the  circuits  providing  threshold  and  plane  wave  editing 
were  only  recently  installed.  The  results  herein  were  obtained  in 
April  and  May  1969  and  include  data  before  and  after  the  operation  of 
the  data  editor  for  medium  an  1  short  range  targets.  Two  medium  range 
targets  are  located  near  Houston,  Texas:  KLC  transmitting  A-l  at  8.666 
MHz  and  KA2XTZ  transmitting  continuous  carrier  at  8.00  MHz,  The 
transmitters  are  about  15  miles  apart  (bearing  and  range  are  noted  on 
the  histograms) . 

( 

The  remaining  three  targets  are  short  range:  Mew  Berlin,  Texas 
located  at  89°  azimuth  with  a  range  of  48  km;  Comfort,  64  km  on  an 

/ 


x  -  — - i^,t- 


azimuth  ol  i-t  /  °  ;  and  Handel. i,  55  .au  mi  an  azimuth  of  iUO  ^  . 


Transmissions 


we  iv  nL  8.5  and  b.  i  Mlli;  with  modulation  either  A-U  or  A-l .  Transmitting 
antennas  were  re.somniL  i/d  wave  length  horizontal  dipoles  oriented  to 
minimize  ground  wave.  A  linear  I'M  LW  ionosonde  transmitter  was  colocated 
at  each  oi  the  shorL  range  targets.  Hie  ionosonde  receiver  was  located 
at  the  DF  siLe. 

Li.  RESULTS 

The  capability  of  the  system  when  the  propagating  medium  is  stable 
is  shown  in  Figure  2.  The  target  is  RLE  at  346  km  range.  Mean  devia¬ 
tion  is  0.18°  while  the  standard  deviation  is  0.37°  for  a  sample  size 
of  150  frames  taken  over  a  period  of  2-3  minuLes.  Subsequent  runs  show 
the  same  stability.  Near  vertical  incidence  ionograms  predicted 
E-layer  propagation  well  beyond  8.0  MHz  for  the  path  in  early  afternoon. 
The  elevation  histogram  is  shown  in  Figure  3.  The  mean  elevation  is 
60.4°  with  a  standard  deviation  of  0.27°.  Note  that  we  have  been  using 
mean  deviation  and  standard  deviation  very  confidently  in  these  examples. 
However,  as  suggested  later,  these  normal  statistical  parameters  must  be 
used  with  care  when  short  range  propagation  is  subjected  to  ionospheric 
tilts  or  other  path  disturbances. 

Figure  4  shows  the  azimuth  histogram  for  the  55  km  target  at 
Bandera  at  8.5  MHz  about  25  minutes  later  than  the  KLC  data  above.  The 
mean  deviation  here  is  4.19°  with  a  standard  deviation  of  3.56°.  Only 
threshold  editing  was  used,  the  plane  wave  option  was  off.  The  dis¬ 
tribution  appears  normal.  The  elevation  histogram  in  Figure  5  shows  a 
mean  elevation  angle  of  86.13°  with  a  small  variance.  About  8  minutes 
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prior  to  Lli  is  da  La  sol  ,  a  similar  run  was  r.iado  using  the  plane  wave 
Lest  option  as  well  as  Lhe  threshold  so  Lliai  data  were  taken  only  when 
the  signal  amplitudes  were  equal  aL  the  spaced  antenna.  Figure  0  shows 
distincL  separation  ot  wIiul  appears  to  be  two  modes.  Although  the  data 
sample  is  the  same  length  and  ol  approximately  the  same  time  duration 
as  the  last  set,  an  overall  mean  and  standard  deviation  do  noL  ade¬ 
quately  describe  what  we  see  here.  (For  the  record,  the  calculated 
mean  was  0.4°  with  a  standard  deviation  of  5.5'.)  lhe  important  charac¬ 
teristic  is  what  appears  to  be  the  mode  separation.  The  elevation 
histogram  in  Figure  7  shows  very  iittie  difference  from  the  last  sample 
although  the  mean  elevation  is  now  84.8°  rather  than  86  lc  which  can 
be  significant  for  target  location  by  vertical  triangulation. 

The  next  data  sets  show  performance  characteristics  for  the  two 
medium  range  targets  near  350  km  and  the  short  range  64  km  target  at 
Comfort.  All  data  were  obtained  on  the  evening  of  13  May  1969  between 
9:00  and  10:30  PM  CDT.  Only  threshold  editing  was  used  since  the 
plane  wave  test  option  was  not  instrumented  at  this  time. 

The  azimuth  histogram  for  KLC  Houston  at  346  km  at  2120  hours  is 
shown  in  Figure  8.  As  before,  the  sample  sizes  are  each  150  frames 
taken  in  a  3-5  minute  period.  As  expected,  there  is  more  spread  in 
night  time  data  than  the  earlier  example.  The  mean  deviation  here  is 
-0.25°  with  a  standard  deviation  of  1.25  .  (In  the  first  series,  Figure 
4,  the  mean  was  4-0.18  while  the  standard  deviation  was  0.37c.)  The 
elevation  data  in  Figure  9  are  uneventful  except  that  it  is  interesting 
that  the  elevation  has  lowered  from  the  earlier  value  of  60°  to  53°, 
which  suggests  that  the  afternoon  propagation  was  probably  two  hop  E- . 
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Azimuth  data  on  the  L  u  i  vers  i  t  y  ul  Houston  LarguL--8.0  MHz  and  li  L>  0  km 
range  l liken  at  21  It)  just  10  minules  prioi  Id  Lite  Kl.C  data  j us L  siiown —  are 
shown  in  Figure  it).  there  is  considerably  more  spread  to  these  data, 
showing  a  a./J  standard  deviation.  The  spread  is  interesting  because 
the  two  targets,  K.LC  and  KA2XTZ,  are  separated  by  not  more  than  11)  miles 
and  oo7  kHz  in  trequency.  Transmitting  character isti'-s  have  not  been 
studied  although  it  is  known  that  we  are  working  oil  a  side  lobe  of 
the  University  of  Houston  antenna  which  is  directional  toward  Uroana 
and  chat  kAZXi'Z  is  transmitting  continuous  carrier  while  KLU  is  1CW. 

Before  we  look  at  the  elevation  data,  a  histogram  of  this  same 
signal  just  prior  to  the  path  going  out  an  hour  later  at  2225  hours 
(Figure  11,  demonstrates  that  comparative  statistics  on  a  fixed  sample 
size  are  no  longer  relevant  since  21  points  were  cast  out  on  the  high 
side.  We  do  want  to  emphasize,  however,  that  there  is  no  reason  to  sus¬ 
pect  a  deterioration  in  array  performance  from  those  earlier  examples. 

The  elevation  measurement  performances  for  these  two  azimuth  sets  are 
shown  in  Figures  12  and  13  The  first--well  before  the  path  went  out-- 
has  a  mean  elevation  of  56°  while  the  second  just  prior  to  fade  out 
shows  considerable  deterioration  but  the  mean  elevation  has  gone  up 
to  59°. 

Azimuth  data  for  the  short  range  64  km  target  at  6  3  MHz  are  shown 
in  Figure  14,  taken  around  2210 — just  15  minutes  before  the  8  MHz, 

350  km  path  went  out.  Conditions  are  still  reasonably  stable.  The  mean 
deviation  is  8.2°  while  the  standard  deviation  is  about  4°,  which  has 
been  considered  extremely  good  in  the  past  for  signals  arriving  at  82° 
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elevation  ^Figure  1 1> )  ;  ho  we  vet ,  ah  we  sail!  earlier,  t 
be  misleading  particularly  lot  aLeepiy  downcomiug  signals  and  need  to  be 
judged  with  cons iderat ion  given  Lo  the  lime  interval  over  which  the  data 
were  taken. 

Ihe  last  sets  ol  data  were  obtained  1  roiu  the  8  5  MHz,  4o  km  target 
at  dew  Berlin  on  1j  April  19011.  These  data  demonstrate  the  adverse  ef¬ 
fects  of  building  histograms  or  calculating  statist itaL  performance  data 
on  high  angle  targets  affected  by  ionospheric  tilts  or  other  perturbances, 
in  Figure  Lb,  the  mean  shown  is  about  -7°  while  the  standard  deviation 
is  9.J0.  Compare  the  data  with  those  of  Figure  17.  which  is  a  20  minute 
sample  but  beginning  about  a  half  hour  later  The  mean  deviation  is 
+8.2°  with  a  standard  deviation  of  8°.  As  shown,  the  data  do  not  ap¬ 
pear  to  be  helpful  but  we’ve  shown  these  to  demonstrate  the  need  for 
filtering  and  placing  more  reliance  on  a  time  series  type  plot. 

Filtering  is  shown  in  Figure  18..  These  data,  taken  between  the 
time  intervals  of  Figures  16  and  17.  represent  a  twelve-minute  interval. 
The  azimuth  and  elevation  variations  (the  top  two  plots)  were  obtained 
by  a  50  point  sliding  average  taken  10  frames  at  a  time.  The  total 
number  of  frames  is  about  800*  Note  that  the  overall  standard  devia¬ 
tion  would  be  relatively  large  with  a  mean  near  zero.  However,  if  we 
had  taken  sample  sizes  of  150  frames  as  we  did  in  the  earlier  data, 
the  standard  deviation  would  b-  low  but  the  mean  would  vary  almost  in 
the  manner  of  the  curve  shown. 

One  of  the  objectives  of  the  HF  interferometer  research  on  targets 
of  short  and  medium  range  is  target  location  by  vertical  triangulation. 
The  mechanics  of  the  problem  are  illustrated  in  Figure  19,  In  order  to 
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OBSERVATIONS  OF  8.5  MHZ  TARGET  AT  NEW  BERLIN 
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Figure  18.  Observations  of 


Figure  19.  Ionosphere  Tilt  Effect. 
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perform  vuiLio.ii  t  r  i  anguiaL  ion ,  angle  ol  arrival  data  must  b .  corrected 
lor  ionospheric  tiLt.  The  virLual  height  or  the  path  length  must  be 
measured.  Other  measured  parameters  are  the  angles  oi  arrival  (airimuLli 
and  elevation).  II  the  target  location  is  not  known,  the  Lilt  angle 
must  be  measured  by  Dl-'ing  on  a  vertical  incidence  probe  or  measure¬ 
ments  from  a  known  station  can  be  used  as  a  correction  factor  applied 
to  azimuth  and  elevation  data  of  the  unknown. 

The  12-minute  example  in  Figure  18  demonstrates  the  calculation 
of  the  tilt  parameter  from  the  known  target  at  New  Berlin  —  a  48  km 

range  —  at  8.5  Milz.  A  linear  FM  CW  ionosonde  transmitter  colocated 
with  the  target  enabled  the  path  delay  or  path  length  to  be  measured 
directly.  Data  from  three  ionograms  were  used  to  plot  the  variation 
of  path  length  shown  as  a  function  of  time  in  Figure  20.  Data  were 
just  after  local  sunset. 

By  knowing  the  deviation  from  the  true  azimuth  and  the  elevation 
angle,  along  with  the  path  length  of  the  tilt  angle,  the  inclination 
of  the  normal  to  the  reflecting  surface  with  respect  to  the  vertical 

can  be  calculated! as  shown  in  the  third  plot  of  Figure  18.  The  calcu¬ 

lated  azimuth  of  the  vertical  plane  that  includes  the  tilt  vector  and 
passes  through  the  midpoint  of  the  GCB  line  is  plotted  on  the  bottom. 

As  an  exercise  and  to  gain  insight  into  the  influence  of  the  tilt 

effect,  we  assumed  the  target  as  unknown  and  calculated  target  loca¬ 
tions  by  vertical  triangulation  using  only  the  measured  parameters  of 
azimuth,  elevation  angle  and  the  near  vertical  incidence  virtual  height 
(assumed  plot  ionospheric  parallel  to  the  earth).  The  variations  in 


location  using  data  that  include  the  12-minute  sample  f  Figure  1 y 


WM  H19N31HlVd 


demonstrate  that  a  single  tilt  co  i  rue  L  ion  lor  a  data  sample  as  short 


as  U  minuLes  would  bo  ol  iitLlo  uso  lor  short  range  targets  The 
results  .suggest  that  1 1 1 L  paraueters  are  required  in  real  time  along 
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DIRECTION  OF  ARRIVAL  STUDIES  ON  A  LONG  RATH  (Summary) 

Lowell  Tveten  and  1  will  be  presenting  some  of  the  results  obtained 
in  two  series  of  UF  directional  propagation  experiments  over  long  range 
paths.  Our  paths  were  orthogonal  to  one  another.  I've  called  the  ex¬ 
periment  a  "cooperative"  one,  Lowell  calls  it  a  "coordinated"  one  1 
don't  believe  there  was  that  much  difference.  To  cooperate  means  "to 
act  or  work  together  with  another  or  others  for  a  common  purpose  "  Co¬ 
ordination  implies  harmonious  adjustment  or  functioning.  From  my 
viewpoint  we  have  had  both  ,  The  relationship  is  continuing  because 
there  is  much  left  to  do  with  the  data  that  were  acquired,  We've  barely 
skimmed  a  small  portion  of  the  available  data. 

The  body  of  the  presentation  was  concerned  with  information  found 
in  the  ATLAS  of  DATA  obtained  in  a  cooperative  HF  directional  propaga¬ 
tion  experiment  over  a  1330  Kilometer  path  (March  1967  -  March  1968) . 
Slides  were  shown  which  demonstrated: 

1.  The  propagation  path  geometry 

2.  The  data  vector 

3.  The  format  for  the  reduced  data 

4.  An  example  of  a  typical  page  of  data. 

1.  "An  Atlas  of  Reduced  Data  Obtained  in  a  Cooperative  HF  Directional 
Propagation  Experiment  Over  a  1330  Kilometer  Path  (March  1967  - 
March  1968)"  A.D.  bailey,  EW  Ernst,  L.J  Miller,  W.W,  Wood, 

RRL  publication  No.  344;  Radiolocation  Research  Laboratory,  De¬ 
partment  of  Electrical  Engineering,  University  of  Illinois,  Urbana, 
Ill inc is . 
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Certain  dominant  luaLur-.w.  in  Lhe  d^ta  were  pointed  mil  as  follows: 

1  At  night,  single  (dominant.)  mode  propagation  was  in  evidence  and 
tile  data  represent  true  angle  oi  arrival  as  received.  The  fluctuations 
in  the  data  are  due  to  the  quasi-per ic>dic  lateral  deviation  effects 
(the  "CU-minute"  or  Ross  eifect)  and  the  diurnal  effect.  Twenty-minute 
and  d-liour  "moving  average  windows"  brought  out  the  significance  of  these 
two  deviations. 

in  tlie  dayLime,  considerable  wave  interference  efiects  were  seen 
between  the  two  or  three  modes  of  propagation  when  all  are  of  approxi¬ 
mately  equal  amplitudes.  These  interference  effects  were  removed  by  the 
"averaging  windows." 

J.  The  importance  oi  the  standard  deviation  of  the  sample  as  an  im¬ 
mediate  indicator  of  wave  interference  was  pointed  out. 

4.  Also  discussed  was  Che  dramatic  appearance  of  sidescatter  par¬ 
ticularly  after  the  main  path  fails  because  the  frequency  is  higher 
than  the  MUF  for  the  direct  path.  Of  course,  the  gain  of  the  system 
had  to  be  increased  to  see  the  scatter. 

The  elevation  angle  measuring  properties  of  the  CDAA  were  pointed 
out.  These  properties,  though  restricted  in  elevat ion,  are  neverthe¬ 
less  useful  at  intermediate  ranges 

The  need  for  amplitude  information  in  oblique  incidence  ionograms 
was  stressed. 

The  need  for  more  rapid  sampling  of  the  data  (greater  than  the 
sampling  rate  shown)  was  stressed^  In  this  connection  one  example  of 
some  results  obtained  from  a  tracking-type  system  was  shown.  This 


spoke  to  the  need  for  rapid  sampling  and  indicated  in  a  quantitative 


way  what  one  must  do  it  lie  would  have  his  "druthers"  lor  the  cases  that 
appeared  to  prevail. 

The  ATLAS  shows  the  way  one  time  series  oi  data  may  be  expected  to 
behave  over  a  twenty-four  hour  period.  it  is  a  conversation  with  "Mother 
Nature."  if  one  asks  the  right  questions  and  is  patient,  lie  gets  back 


a  lot  of  answers. 


The  general  plan  and  objectives  of  a  coordinated  cooperative 
experiment  performed  by  the  University  of  Illinois  and  the  Institute 
for  Telecommunication  Sciences,  E8SA  Research  Laboratories,  are  dis¬ 
cussed.  Results  so  far  indicate  that  it  is  feasible  to  use  back- 
scatter  to  assess  DF  reliability  over  simple  one-hop  paths.  Until 
more  is  understood  about  the  complex  interaction  between  radio  waves 
and  the  irregular  ionosphere,  it  will  not  be  generally  possible  to 
correct  DF  measurements  on  the  basis  of  backscatter  irregularity 
patterns.  Larger  scale  tilt  determinations  by  a  backscatter  mapping 
technique  may  ultimately  prove  useful  in  DF  correction.  * 


*  This  paper  has  not  been  released  for  publication  by  ESSA  and 


should  not  be  referenced. 
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1  I  .  SUMMARY 

The  coordinated  HR  propagation  experiment  was  designed  to  obtain 
correlations  between  propagation  parameters  and  d  t  rec L ion-o f-arr ival 
(L)OA)  measurements  by  a  mul L ipl ie i L y  oi  observations  and  observing 
techniques,  the  institute  tor  Telecommunication  Sciences  (ITS)  of 
USSA,  the  University  ol  Illinois,  and  the  University  of  Houston  co- 
opeiated  in  the  experiment.  The  center  ior  the  ITS  observations  was 
tne  high- resolution  antenna  array  at  Table  Mountain  near  boulder,  shown 
in  Figure  i.  The  azimuth  array  forms  a  beam  very  narrow  in  azimuth 
and  wide  in  elevation.  The  elevation  array  forms  a  beam  narrow  in 
elevation  and  wide  in  azimuth.  The  azimuth  array  has  an  aperture  of 
1,392  feet.  The  elevation  array  has  an  effective  aperture  of  1,000 
feet,  since  it  is  mounted  on  a  500-foot  tower  with  the  beams  formed 
in  s  .•.h  a  way  that  both  direct  and  ground-reflected  rays  are  used. 

These  arrays  have  a  nominal  frequency  range  of  12  to  25  MHz  and  equiva¬ 
lent  beam  widths  from  1.5  to  3  degrees  in  the  case  of  the  azimuth 
array,  and  about  2  to  4  degrees  in  the  case  of  the  elevation  array. 

The  effective  scan  sector  for  the  azimuth  array  varies  between  40  and 
90  degrees  and  the  effective  scan  sector  for  the  elevation  beam  varies 
between  22  and  52  degrees,  approximately.  Figure  2  shows  a  view  of 
the  scan  and  beam-forming  elec.ronics  inside  the  building;  not  shown 
are  the  various  recording  interfaces  and  timing  systems  for  the  re¬ 
cording  and  processing  cf  the  data. 

Figure  3  illustrates  some  of  the  types  of  experiments  and  experi¬ 
mental  observations  that  were  conducted.  At  boulder,  high-resolution 
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Figure  1.  Azimuth  and  Elevation  Arrays  at  Table  Mountain  Field  Site. 
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I.T.S.A.  HF  BACKSCATTER  AND 
DIRECTION -OF-ARRIVAL  INVESTIGATION 


Station 

Instrumentation 

BO- Boulder,  Colo. 

V-I  Sounder,  HF  high  resolution 
backscatter  sounder  and 
receiving  •  system,  Granger 
oblique  receiver 

PC-Ponca  City,  Okla. 

V-I  Sounder 

SP-Sprinqfield,  Mo. 

V-I  Sounder  &  C  W  Doppler 

TX-Texarkana,  Ark. 

II  n  H  II  II 

GW-Greenwood ,  Miss. 

n  n  M  n  ii 

UR-Urbana,  III. 

V-I  Sounder,  DF  system, Granger 
receiver,  CW  Doppler 

AL-Alta  Loma,  Tex. 

Granqer  oblique  sounder 

MR-Meridian,  Miss. 

Transponder 

OR-Orlando,  Fla. 

Transponder 

LR-Little  Rock,  Ark. 

Earth  currents 

AN-Annapolis ,  Md. 

Oblique  sounder 

Figure  3.  I.T.S.A.  HF  Backscatter  and  Direction-of-Arrival  Investigation 
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azimuth  and  I'li'V.iLiuii  scan  backs'  a  L  t  c  i  recordings  were  made  on  five 
i requeue ies  .  UUA  and  signal  strength  lor  cji.li  ol  LWu  propagation  modes 
from  each  e> l  two  transponders,  one  iocaLed  aL  Meridian,  Mississippi, 
and  the  other  at  Orlamlo,  Florida,  were  made  Vertical  sounder  stations 
were  located  at  l’onca  City,  Oklahoma;  Springfield,  Missouri;  Texarkana, 
Arkansas;  and  tireenwood,  Mississippi.  in  addition,  steep- incidence 
phase  measuring  transmitters  were  installed  at  lireenwood,  Texarkana, 
and  Springfield  and  received  at  Little  Rock,  Arkansas,  The  University 
of  Illinois  made  00A  measurements  between  AlLa  Lo.na,  Texas,  and  Irbana, 
Illinois,  ft  is  obvious  the  two  propagation  paths  crossed  at  very 
nearly  right  angles.  Oblique  sounder  transmissions  from  Texas  were  re¬ 
ceived  at  both  Urbana  and  Boulder,  Measurements  were  made  24  hours  a 
day  on  5  consecutive  days  each  month  from  August  1966  to  June  1968. 

The  ITS  objectives  were:  (1)  to  determine  the  feasibility  of  using 
backscatter  for  (a)  assessing  DUA  measurement  quality  or  (b)  correcting 
specific  OF  readings;  (2)  to  check  the  ITS  3-D  ray-tracing  program; 

(i)  to  identify  and  measure  propagation  modes  and  to  obtain  experimental 
data  for  testing  particular  processing  and  analytical  techniques;  and 
(4)  to  study  anomalous  propagation,  specifically  sidescattexA 

Figure  4  illustrates  simultaneous  azimuth  and  elevation  scan  back¬ 
scatter  records.  The  horizontal  axes  are  labeled  in  degrees  for  both 
azimuth  and  elevation;  the  vertical  axis  is  labeled  in  milliseconds  As 
one  might  expect,  the  backscatter  shows  high  elevation  angles  at  .".he 
smaller  delays  and  lower  angles  as  the  delay  increases  -  the  classic 
elevation  angle-backscatter  delay  relationship.  The  azimuth-versus-cielay 
presentation  on  the  left  shows  a  number  of  intensifications  caused  by 


ionospheric  disturbances. 


Figure  4.  Azimuth  Elevation  Scan  Backseatter  Records  at  15.67  MHz  Taken 

at  1215  and  1630  MST  November  19,  1966.  Skip  Echoes  are  Saturated. 


I'o  chock  i  ho  lejii  iiliLy  ol  using  backscatter  to  assess  OF  qualiLy, 


0-1 


the  backseat  tor  data  lor  each  hour  ol  t ho  day  wore  categorized  according 
to  live  cl  ass i t  icaL ions  illustrated  in  Figure  j  The  precise  charac¬ 
teristics  ol  those  categories  are  not  especially  important  They  ap¬ 
peared  consistently  throughout  the  live  months  (distributed  over  a  year's 
time)  used  in  the  analysis  UoX  data  were  class  Hied  as  "bad"  if  the 
deviations  exceeded  one  degree  during  a  given  hour  and  as  "good"  other¬ 
wise.  bach  hour's  data  were  compared  with  respect  to  L)F  quality  back- 
scatter  category  .  Results  of  this  comparison  are  shown  in  Table  1. 

Table  1,  Meridian  One- Hop  F  UUA  Compared  with 
Backseat  ter  Characteristics 


A 

Back 

B 

scatter 

C 

Category 

D 

E 

Good  DOA,  per  cent 

19 

27 

37 

100 

100 

Relative  occurrence,  per  cent 

19.3 

28c  9 

VA 

20,2 

14  2 

Order  of  preferred 
distributions 

4 

3 

2 

lb 

la 

The  traveling  ionospheric  disturbance  category  has  only  19  per  cent 
good  DF  data,  Backscatter  category  E  has  100  per  cent  good  data,  with  a 
relatively  smooth  ascending  percentage  between  categories  A  and  E  The 
next  row  shows  the  relative  occurrence  of  the  various  categories,  and 
we  see  that  the  distribution  is  relatively  uniform  over  the  time  period. 

The  last  row  shows  the  order  of  preference,  if  one  were  to  assess  the 
quality  of  the  DOA  data  with  respect  to  the  various  backscatter  categories. 
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Wo.  therefore  conclude  thai  it  is  indeed  ieasitilo  lo  use  buckscatter  to 
assess  LH)A  measurement  quality  or  accuracy  at  the  same  time  and  in  the 
same  general  region  over  simple  one-hop  I'  patlis 

To  actually  correct  DUA  measurements  Is  a  bit  more  complex  Une 
would  expect  that  since  Lhe  toe using  irregularities  for  backseat  ter  are 
the  same  irregularities  that  cause  deviations  in  L)OA,  there  would  be  a 
rather  specific  relationship  between  the  two.  To  check  this  relation¬ 
ship,  we  examined  several  months' scan  backscatter  data  and  defined  the 
times  at  which  a  specific  irregularity  appeared  at  the  apparent  midpoinL 
of  the  path  between  boulder  and  Meridian,  Mississippi  The  DOA  data 
were  examined  for  the  same  time.  In  general,  little  correlation  was 
noted  except  in  the  December  1967  data,  in  which  the  DOA  deviated  to 
the  south  within  two  minutes  of  an  irregularity  over  path  midpoint  in 
50  per  cent  of  the  cases.  The  probability  of  this  occurring  on  the 
basis  of  random  distributions  is  extremely  small. 

To  gain  insight  into  the  problem,  we  did  some  ray  tracing  with 
ionospheric  irregularity  models.  Figure  6  is  an  example  of  simultaneous 
azimuth  and  elevation  DOA  determined  from  ray  tracing  a  moving  irregu¬ 
larity  model.  It  can  be  compared  directly  in  Figure  7  with  an  experi¬ 
mental  record  taken  during  December  1967  in  which  a  number  of  TID’s 
are  moving  through  the  area.  The  two  are  quite  similar;  however,  one 
can  obtain  different  relationships  in  the  ray-tracing  outputs  simply  bv 
using  the  different  heights  of  travel  for  the  disturbance  with  respect 
to  the  ray  path  geometry,  by  changing  the  direction  of  travel  of  the 
disturbance  somewhat  with  respect  to  the  ray  geometry,  or  by  changing 
the  type  of  irregularity  from  an  ionization  intensification  model  to  a 
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Figure  5.  Kx.imples  of  b  ickscat Ler  Categcr  les . 
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raretaet  ion  moil  cl.  WiLli  all  Lhese  possibil i L i us ,  iL  Is  not  difficult 
to  understand  why  many  diliurent  re  Lit  ionsliips  can  be  seen  in  the  actual 
records  when  one  compares  backscaLLer  irregular i  Ly  structures  aiu  ■  A 
behavior  and  even  behavior  between  the  azimuth  and  elevation  angle  de¬ 
viations.  On  the  basis  oi  what  we  know  so  far,  we  conclude  that  it  is 
not  usually  possible  to  make  specific  DOA  corrections  on  the  basis  of 
irregularity  positions  on  backscatter  records.  Ultimately  however,  witli 
more  knowledge  oi  t he  mechanisms  of  interaction,  one  may  be  able  to  go 
from  the  backscatter  pattern  back  to  the  generating  ionospheric  irregu¬ 
larity  and  so  deduce  the  correction  necessary,  but  much  more  work  in 
modeling  and  ray  tracing  is  required  for  knowledge  of  this  kind 

The  ionosphere  may  be  thought  of  as  being  made  up  of  a  spectrum 
of  ionospheric  irregularities  with  periods  of  less  than  a  second  to 
several  hours.  Another  application  or  possibility  for  DF  correction 
is  to  correct  for  the  longer  period  or  larger  area  tilts..  This  could 
perhaps  be  done  by  mapping  the  generalized  ionization  densities  by  the 
backscatter  technique.  Figure  8  is  another  azimuth  and  elevation  back¬ 
scatter  record  that  shows  a  very  classic  elevation  angle  versus  delay 
relationship.  From  this  relationship,  in  conjunction  with  transmission 
curves  for  an  assumed  ionosphere  model,  one  can  deduce  plasma  frequencies 
versus  h'  at  various  locations  along  the  path  corresponding  to  the 
different  delays.  The  use  of  multiple  frequencies  allows  these  plasma 
frequencies  to  be  determined  at  different  height  levels  in  the  ionosphere. 
With  a  generalized  mapping  program,  one  can  map  the  virtual  height  versus 
plasma  frequency  over  a  large  geographic  area  to  be  subsequently  used 
for  defining  the  broad  general  overall  tilts  over  a  rather  large  area 


for  correction  purposes. 
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Figure  8.  Simultaneous  Azimuth  and  Elevation  Scan  Display  of  Backscatter 

at  about  15.5  MHz  Taken  during  a  November  Day  in  19bb  at  1405  MSI. 
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For  such  mapping,  one  needs  azimuthal  discrimination  for  the  ele¬ 
vation  ungLe  versus  delay  time  response  lor  backseat  Let’ ,  wliicii  can  be 
achieved  by  cross-correlat  tug  Lise  uximuLh  and  elevation  beam  outputs  to 
form  an  eilective  pencil  beam,  a  method  LiiaL  lias  been  attempted  at  lid 
It  has  not  proved  wholly  successful  because  ot  our  array  configuration 
and  spacing,  but  the  addition  ot  a  lew  more  elements  to  Lhe  azimuth  ar¬ 
ray  on  the  unoccupied  side  oi  Lhe  300-Loot  tower  would  we  think,  bring 

SUCCESS „ 

One  could  also  apply  the  mapping  technique  if  the  transmitting 
beam  were  sufficiently  narrow.  Using  a  rhombic  antenna  as  an  illumi¬ 
nator  for  transmission,  we  have  made  estimates  of  plasma  frequencies 
versus  height  from  backscatter  records  and  compared  them  with  actual 
plasma  frequencies  versus  height  at  particular  ionosonde  locations 
The  errors  have  been  on  the  order  of  10  to  15  percent  and  we  believe 
that  with  refinement  of  the  technique  they  can  be  reduced  considerably. 

Another  objective  was  to  check  the  3-D  ray-tracing  program  This 
we  did  by  reducing  vertical  ionosonde  data  to  true  height  profiles  at 
several  locations,  connecting  them  smoothly  over  the  area  in  question 
with  a  generalized  mapping  program,  and  ray- tracing  through  the  re¬ 
sulting  ionosphere.  The  directions  predicted  from  the  ray-tracing 
program  were  compared  with  the  actual  measurements  made  at  the  same 
time  over  the  path,.  Figure  9  is  an  example  of  this  type  of  comparison 
The  continuous  lines  show  the  measured  directional  variations  of  azimuth 
and  elevation,  and  the  dots  are  the  calculated  points  at  the  specific 
times.  The  agreement  is  not  exact,  but,  under  the  circumstances,  is 
considered  close  enough  to  validate  the  ray-tracing  program  Other 
examples  have  been  checked  and  found  to  be  in  reasonable  agreement  also. 


Universal  Time,  September  21, 1967 


Comparison  of  actual  c!irection-of-ar rival  (continuous  curve) 
and  di  rec  tio  n  -  of -a  r  r  i  val  determined  by  ray  tracing  true  height 
mapped  ionosphere  from  vertical  ionograms  (dots). 
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Another  area  a!  cl  Lai  l  was  i  lie  dc  Lc  rm  i  na  L  ian  ol  propagation  moUvs 
by  elevation  -ingle  measui  eiiients .  Ihe  Lu.  Ini  iipie  lor  doing  this  is  to 
generate  a  large  number  o 1  ground  rangc-vci sus-elcvnl ion  angle  curves 
lor  dittet'ent  possibLe  propagation  modes  and  i  oiiibinat  ions  of  prupagutiui 
modes.  Since  Lhe  source  el  L  ansmission  must  be  at  one  point,  a  clustei 
ol  these  modes  must  indicate  a  possible  source.  1 i ,  in  addiLiun.  back¬ 
seat  ter  data  are  used  to  assess  the  propagation  conditions  of  the  iono¬ 
sphere,  a  particular  range  can  be  chosen.  This  choice  determines  the 
range  to  the  site  as  well  as  the  identification  of  the  pertinent  propa¬ 
gation  modes.  We  compiled  a  set  of  data  from  our  tran.s-  nder  locations 
using  elevation  angle  and  backseat  ter  measurements  and  derived  range 
estimation  errors  with  a  standard  deviation  of  about  d  percent.  We 
believe  this  can  possibly  be  reduced  even  further  and  are  currently 
working  on  an  automated  computer  technique  for  making  these  estimates 
to  determine  whether  we  can  duplicate  the  original  accuracy. 

figure  10  shows  another  azimuth-elevation-versus-delay  scatter 
record  that,  in  this  case,  is  sidescatter  from  the  region  of  the  Gulf 
of  Mexico  from  transmissions  originating  near  Havana.  Illinois.  At  the 
present  time,  we  are  endeavoring  to  make  quantitative  measurements  of 
sidescatter  parameters. 

As  useful  fallout  from  this  experiment,  we  are  planning  to  extract 
DF  statistics  concerned  with  mode  of  propagation  and  periods  of  the 
variations  and  to  perform  correlations  with  various  geophysical  para¬ 
meters,  etc.  We  hope  to  obtain  statistics  on  Till  occurrence  and  to 
learn  much  more  about  the  mechanisms  oi  wave  interaction  with  the 
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Figure  10.  Azimuth  and  Elevation  Record  of  Direct  and  Sidescatter  Modes 

Received  from  Havana,  Illinois,  at  a  Frequency  of  19.36  MHz  at 
2250  GMT  on  March  30,  1967.  Left  Side  of  Azimuth  Scale  is  Labeled 
in  Azimuth  of  First  Grating  Lobe  on  the  Right  as  Well  as  Main 
Lobe . 
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Other  jppi  l  ev  l  ions  I  ji  L  lie  data  wiJl  duulalcbb 


arise.  Sim  c  lias  expel  imenL  .ovets  nearly  a  J-ycar  |j<_-  i  i  >  >  d  ,  tin.  slo- 
lL->lics  slum  I  •  he  tail  tv  n.jelul. 

La  summa  i  y  t hen ,  we  may  b>ay  that  backseat  Li  t  can  he  Ubcd  Lu  ab.be 
reliability  ot  aa'iii'jiy  ol  t'OA  mcasuied  neat  U  ie  location  or  j’va  ob¬ 
served  by  the  backseat ler  itowevei .  s pec L l i c  backseat ie r  irrogularit 
coni  igurat  ions  camu>L  yet  be  used  to  col  lect  IJUA  data.  We  are  opti¬ 
mistic  t ha L  a  general  mapping  program  based  on  backseat Ler  elevation 
angle  information  may  be  used  to  correct  lor  large  generalized  Lilts. 
However,  the  necessary  implementation  in  the  form  oi  a  Mills  T  pro¬ 
cessing  aiiay  still  requires  some  Lurther  development.  I  lie  IIS  3-D 
ray- tracing  program  appears  to  be  an  accurate  and  useful  tool  for 
ionospheric  study  and  it  does  check  well  with  experimental  data  We 
are  hopeful  that  backseatter  and  elevation  angle  information  can  be 
csed  to  estimate  range  rather  accurately  Sidescatter  lias  been  ob¬ 
served  rather  extensively,  and  we  expect  to  define  its  various  char¬ 
acteristics  to  make  it  more  readily  identifiable  to  DF  systems 
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KAY  TRACING  blMLLATTON  OF  IONOSPHERIC  EFFECTS  IN  1JICII  FREQUENCY 
RAD 10  LOCATION 

L.  SUMMARY 

Deviation  of  high  1 requuncy  radio  waves  l rum  the  great  circle  plane 
containing  the  transmitter  and  tl;e  receiver  due  to  horizontal  gradients, 
irregularities  and  other  effects  in  the  ionosphere  is  a  dominant  factor 
in  HF  radiolocation.  In  order  to  correct  for  the  deviations  in  the 
experimentally  measured  angles  of  arrival  and  determine  the  unknown 
location  of  the  source,  it  is  first  necessary  to  fully  understand  the 
role  played  by  the  ionosphere.  This  is  done  by  the  application  of 
radio  ray-tracing  techniques  to  signals  from  sources  of  known  location 
The  Radiolocation  Research  Laboratory,  in  cooperation  with  the 
ESSA  Research  Laboratories,  conducted  a  series  of  experiments  in  which 
the  direction-of-arrival  observations  of  high  frequency  signals  trans¬ 
mitted  from  a  site  near  Houston,  Texas,  were  made  continuously  for  a 
period  of  one  week  per  month  using  the  Wullenweber  array  at  the  Bond- 
ville  Road  Field  Station  during  the  period  March,  1967  -  March,  1968. 
This  paper  is  concerned  with  the  ray  tracing  simulation  of  certain 
ionospheric  effects  observed  during  this  experiment.  The  bearing 
deviations  obtained  from  the  ray  tracing  simulation  are  compared  with 
the  measured  values  during  The  same  period. 

Ionospheric  effects  on  the  direction-of-arrival  of  hf  signals  due 
to  the  following  phenomena  are  discussed: 
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SHORL  lUNlih  liu;ii  ANULh  Si'L’UlhS 

La  attempting  to  perform  short-range  radiolocation,  it  is  neces¬ 
sary  to  work  with  a  signal  which  arrives  at  the  radiolocation  site  at 
angles  near  the  zenith  (directly  overhead),  bignals  i  rom  less  than 
100  kiu  distant  retlectiug  from  an  ionosphere  aL  a  virtual  height  of 
more  Lhan  lot)  km, lor  example,  will  have  this  property.  The  ray  path 
geometry  is  shown  in  Figure  L.  This  signal  will  also  consist  of  two 
independent  components,  tne  ordinary  component  and  the  extraordinary 
component . 

Also,  since  a  radio  signal  is  polarised  in  the  plane  normal  to  its 
propagation  direction,  it  will  appear  horizontally  polarized  It  is 
this  kind  of  signal  that  must  be  used  in  short  range  RDF. 

This  set  of  conditions  makes  several  demands  on  the  antennas  and 
indicating  equipment.  The  antennas  must  be  horizontally  polarized  so 
that  sensitivity  is  not  lost.  They  should  also  be  circularly  polarized 
so  that  the  ordinary  mode  (left-hand  circularly  polarized)  and  the 
extraordinary  mode  (right-hand  circularly  polarized)  do  not  interfere 
with  each  other.  Finally,  they  should  be  omnidirectional  so  that  sig¬ 
nals  from  all  directions  will  be  acceptable, 

The  expected  direction  of  arrival  of  the  signal  also  places  some 
special  requirements  on  the  indicating  equipment,  It  should  display 
the  incidence  angle  as  well  as  the  bearing  of  the  received  signal 
This  incidence  angle  information  can  be  used  to  make  an  estimate  of  the 
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range  ot  L  lie  transmitting  station  when  coupled  with  suiue  available 
ionospheric  inturmat ion. 

an  example  o  1  a  sysLem  ot  Lliis  kind  is  shown  .11  Figure  2  'lliis 
l Lgure  ts  a  block  diagram  oi  Lhe  University  ul  Illinois  triple  radio 
inter t urometer  system.  A  basic  set  ot  Llit'ee  antennas  is  used  to  give 
three  interterometer  pairs.  iliese  Lhree  pairs  ol  antennas  are  Lime- 
shared  by  a  dual  channel  receiver  and  phase  display  system.  The  direc¬ 
tion  ot  arrival  is  computed  with  the  available  analog  computer  or  by  a 
digital  computer  using  the  phase  information  stored  on  punched  paper 
tape. 

Some  more  information  should  be  given  on  the  antennas  which  are 
used  in  this  system.  Each  antenna  is  a  pair  of  crossed  loops  connected 
to  a  quadrature  hybrid  coupler.  The  result  is  a  horizontally  polarized 
antenna  with  two  outputs,  a  right-hand  circularly  polarized  output  and 
a  left-hand  circularly  polarized  output.  We  may  thus  choose  the  polari¬ 
zation  we  want  at  any  time.  These  antennas  are  also  omnidirectional  in 
bearing  angle  and  nearly  omnidirectional  in  zenith  angle.  Thus,  all  of 
the  antenna  requirements  for  radio  direction  finding  at  short  range  and 
high  angles  have  been  met.  The  physical  separation  between  adjacent 
antennas  is  100  meters. 

The  signals  from  these  antennas  are  received  as  three  pairs  of 
signals  using  time  sharing  of  a  dual  channel  receiver.  The  three  dif¬ 
ferencial  phase  angles  are  presented  as  shown  in  Figure  2  to  an  analog 
system  and  (via  punched  paper  tape)  to  a  digital  computer.  These  com¬ 
puters  solve  equations  such  as  shown  in  Figure  3  and  display  the  angles 


of  arrival  of  the  signal. 


To  be  used  subject 
to  the  criterion: 


Figure  2.  Functional  Block  Diagram.  Triple  Interferometer 
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Angle  vjr  to  be  measured 


Figure  3.  Interferometer  Pair. 


One  other  piece  ol  eijuipmenL  is  used  in  the  short  range  studies. 
Located  at  the  RDF  siLe  is  a  vertical  incidence  ionosonde.  ionograius 
are  instantly  available  and  also  recorded  on  lilm  ior  further  analysis. 

these  data  can  now  be  used  Lo  estimate  the  range  Lo  the  transmitter, 
the  first  approximal ion  is  the  so-called  "flat  ionosphere"  model  shown 
in  Figure  a.  l'he  angle  ot  incidence  is  used  in  conjunction  with  the 
virtual  height  for  that  frequency  to  find  Lite  range  using  the  geometry 
ot  Figure  a.  lilt',  bearing  to  the  transmitter  is  just  t he  bearing  of  the 
signal.  We  have  thus  "located"  the  transmitter  using  a  single  RDF  site. 

More  sophisticated  techniques  may  be  used,  of  course,  beyond  this 
first  approximation.  We  may  go  one  step  beyond  the  assumption  of  a 
horizontal  flat  mirror-like  ionosphere  and  use  a  "tilted"  mirror.  This 
tilted  mirror  model  can  be  used,  as  a  first  approximation,  to  absorb 
all  sources  of  deviation  of  the  signal  from  the  flat  ionosphere  model 
if  some  method  of  determining  the  "apparent  tilt"  parameters  of  the 
ionosphere  can  be  found. 

Even  more  sophisticated  models  are  possible.  Three-dimensional 
ray’  tracing  using  some  method  of  estimating  the  lateral  as  well  as 
vertical  electron  density  distribution  would  perhaps  yield  a  better 
estimate  of  the  position  of  the  transmitter. 

A  typical  short  range  radiolocation  experiment  was  performed  on 
o  March  1969.  Two  transmitters  were  set  up  as  shown  in  Figure  5.  Data 
were  taken  at  alternating  times  on  the  target  transmitter  and  on  the 
beacon  transmitter.  The  flat  ionosphere,  model  used  on  the  target  trans¬ 
mitter  produced  the  results  shown  in  Figure  6.  Note  that  the  average 
fix  is  about  18  km  from  the  true  position.  Next,  data  from  the  beacon 
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i TARGET  (MATTOON) 


Figure  5.  Geometry  of  Experiment. 
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Figure  6.  Fix  on  Mattoon  Lake  before  Correction. 


transmit  let  were  used  to  compute  '.lie  ionosphc  r  it  i  i  1  L  parameters.  ihe 
tilted  ionosphere  model  is  shown  in  Figure  /.  This  tilted  model  was 
Lhen  used  to  recalculate  the  position  ot  the  targeL  transmitter.  'lhe 
results  are  shown  in  Figure  S  Mote  that  the  average  I  ix  is  now  only 
about  j  ku.  I  rom  the  true  position  thus,  we  were  able  to  improve  on 
the  tLat  ionosphere  model  in  'his  example. 

More  work  oi  this  kind  is  being  done  The  ti  l  Led  ionosphere  model 
and  some  of  the  more  sophisticated  models  need  to  be  fully  studied  if 
thev  are  to  be  used  in  the  future. 
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Figure  8.  Fix  on  Mattoon  Lake  after  Corrections. 
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i.  1M  KOULX'  1  ION' 

A  dual  i nte rf e remote r  radio  direction  finding  (RDF)  system  utilizc-.s 
both  a  small  aperture  and  a  large  aperture  antenna  array.  The  distance 
between  the  antennas  of  a  small  aperture  array  is  less  than  one-half 
wavelength  over  the  range  of  operating  frequencies.  The  opposite  is 
true  for  a  large  aperture  array.  in  the  past  Lhe  small  aperture  (or 
large  aperture)  distance  was  usually  dictated  by  considerations  other 
than  those  dealing  with  the  relationships  between  the  antenna  aperture 
and  the  error  in  the  bearing  information .  The  need,  therefore,  exists 
to  determine  this  relationship. 

This  paper  is  a  study  of  the  error  generated  in  the  angles  of 
arrival  and  how  it  varies  with  the  aperture  size  of  the  antenna  array. 
This  relationship  will  be  examined  under  the  effects  of  signal  noise, 
equipment  noise,  and  interference  by  a  secondary  wave.  A  three  antenna 
interferometer  system  will  be  simulated  using  a  digital  computer.  The 
root  mean  square  (RMS)  error  of  the  angles  of  arrival  will  be  calculated 
at  several  aperture  sizes  ranging  from  1/20  wavelength  to  5  wavelengths 
The  relationship  between  the  (RMS)  error  and  aperture  size  will  be  com¬ 
pared  with  the  relationship  obtained  from  a  statistical  approach  to  the 
problem  in  order  to  corroborate  the  results  of  both  approaches- 

The  simulation  model  used  in  this  paper  is  based  on  the  work  by 

Church  which  unified  and  expanded  part  of  the  previous  works  of 

„  (2)  _  .  (3)  ,  ...  (4) 

Creasy  ,  Grush  ,  and  Allen 

*Submitted  in  partial  fulfillment  of  the  requirements  for  the  degree  of 
Master  of  Science  in  Flectrical  Fngineering  .  Mr.  Talbott  is  currently 
with  TRW,  Redondo  Beach,  California. 
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The  Interferometer  RDF  system  studied  in  this  paper  uses  three 
monopole  antennas  located  as  shown  in  Figure  1.  It  will  be  assumed 
thac  two  plane  waves  are  incident  upon  the  antenna  array.  The  stronger 
amplitude  wave  is  designated  as  the  primary  wave  and  the  other  is  called 
Che  secondary  wave. 

The  ray  path  of  either  wave  is  defined  to  be  any  line  which  is 
perpendicular  to  the  phase  front  of  that  wave.  The  azimuthal  angle 
(:■)  of  an  incident  wave  is  taken  to  be  the  angle  between  the  North 
reference  line  and  the  ground  plane  projection  of  the  ray  path  as 
shown  in  Figure  1.  The  azimuthal  angle  can  vary  between  +  n  and  is 
measured  positive  in  the  clockwise  direction.  The  incidence  angle  (6) 
is  the  angle  between  a  normal  to  the  ground  plane  and  the  ray  path. 

The  incidence  angle  can  vary  between  0  and  n/2  and  is  measured  positive 
from  the  normal. 

The  distance  (D)  measured  between  antennas  1  and  2  or  antennas  1 
and  3  is  called  the  aperture  size  of  the  interferometer.  Although 
the  angle  (y)  between  antenna  pairs  (1,2)  and  (1,3)  is  assumed  to  be 
60°  in  this  paper,  it  can  have  other  values. 

In  general,  the  directions  of  arrival  of  the  primary  wave  and  the 
secondary  wave  will  differ.  The  direction  of  arrival  is  determined  by 
the  azimuthal  angle  and  the  incidence  angle  of  the  wave.  The  difference 
in  the  angles  of  arrival  of  the  two  waves  is  caused  by  ionospheric 
disturbances.  Thus,  the  system  sees  a  secondary  wave  which  traveleu 
along  a  different  propagation  path  than  the  primary  wave.  The 


X  =  it  cos  a  sin  0 
p  P  P 


X  =  it  cos  a  sin  0 
s  ss 


Y  =  rr  cos  (a  -  y)  sin  0 
p  P  P 


Y  =  rr  cos  (a  -  y)  sin  0 
s  s  s 


and  Re  is  the  real  part  operator.  The  other  symbols  not  previously  de¬ 


fined  are  defined  as  follows: 


w  =  signal  frequency,  in  radians/sec 


t  =  time,  in  seconds 
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\  =  wavelength,  In  meters 

h  =  relative  amplitude  factor,  where  0  Hi  J;  1. 

The  antenna  voltages  can  be  expressed  in  phasor  notation  as 


follows : 


=  E0  +  h‘E0 

-  -  1NX  ~  1  (^  +  NX  ) 

E2  =  E0  eJ  p  +  h-E0 

-  pr  jNY  ,  ,  p:  J  ('t>  +  nys) 

i?  _  v  -•  o  p  +  h  •  E0  e 


E3=  E 


For  simplicity  let  us  normalize  those  phasor  voltages  by  setting 
E0  =  1.  We  obtain: 


El  a  [1  +  h-cos  4]  +  j  [h.sir  <*>] 


(1-1) 


E  =  [cos  NX  +  h*  cos (4)  +  NX  )]  +  j  [sin  NX  +  h-sin(<j>  +  NXg)  ]  d  2) 
J2  P  &  v 

E  =  [cos  NY  +  h-cos(«|>  +  NY  )]  +  j[sin  NY  +  h-sin(d)  +  NYg)]  (1-3) 

3  P  s  H 

Creasy(2)  determined  (using  the  magnitudes  of  the  phasor  voltages 
and  the  phasor  difference  of  the  voltages  at  antennas  i  and  j)  the 
phase  difference  (B^)  between  the  induced  voltages  of  the  two  antennas. 

The  phase  difference  is: 

«  “  3  lEj2  +  l^l2  -  IE.  -  EiP 


core  E'  =  Ej  eJ 


90' 


Equations  (1-1)  through  (1-4)  indicate  that  the  phase  difference 
between  antenna  voltages  (B^)  is  a  function  of  the  phase  difference 


(.$)  between  the  two  Incident  waves. 


In  his  investigations  Click 


assumed  that  the  relative  phase  could  be  approximated  to  be  uniformly 
distributed  if  the  relative  phase  of  the  Incoming  wave  is  sampled  at 
uniform  intervals  of  time.  Applying  this  assumption  Glick^^  showed 


that  the  average  (with  respect  to  j)  phase  difference  (AVB^)  was  a 


function  only  of  N  and  the  primary  angles  of  arrival.  Using  Glick's 


.-<5) 


results,  Church'  showed  that  the  average  phase  difference  between 


antennas  is : 


AVB.„  =  NX  =  N  n  cos  a  sin  d 

12  p  p  p 


(1-5) 


AVB. .  =  NY  =  N  n  cos (a  -  y)  sin  0 
13  p  p  p 


(1-6) 


Solving  equations  (1-5)  and  (1-6)  for  a  and  0^  yields: 

-1  ^B13  ~  ^B12  C0S  y 

a  =  tan  [ - -rrr - : - ] 

p  AVB  „  sm  y 


(1-7) 


i  =  tan  - 

[i1"1  c"  V  .  !  ,1/2 
1  avb12  1 


(1-8) 


where  in  equation  (1-8)  ij  that  determined  in  equation  (1-7). 


CHAPTER  2 


SIMULATION  ANALYSIS  OF  APERTURE  SIZE 


')  j 


The  interferometer  RDF  system  lias  been  simulated  with  a  digital 
computer.  The  basic  structure  of  the  simulation  parallels  the  more 
rigorous  simulation  scheme  used  by  Church  The  simulation  done  in 

this  paper  was  considerably  simplified  through  the  use  of  a  priori 
information.  The  simulation  has  a  closed  loop  form  in  that  angles  of 
arrival  are  used  to  calculate  the  same  angles  of  arrival.  Any  devi¬ 
ations  of  the  angles  after  simulation  from  the  values  before  simulation 
are  considered  to  be  errors  due  to  the  interferometer  system.  The 
simulation  of  an  interferometer  system  has  many  parameters  all  of  which 
have  an  effect  on  the  final  results.  Consideration  of  all  of  them  and 
their  interactions  with  one  another  would  be  difficult.  In  this  paper 
the  primary  angles  of  arrival,  aperture  size,  signal  noise,  and  wave 
interference  are  considered  to  be  of  major  interest,  while  the  secondary 
angles  of  arrival,  phase  meter  and  system  noise,  signal  frequency,  and 
perturbations  of  the  relative  phase  between  the  primary  and  secondary 
signals  are  considered  to  be  of  lesser  interest. 

In  the  interferometer  simulation  by  Church ^  two  problems  of 
considerable  importance  have  to  be  resolved.  The  first  problem  has  to 
do  with  a  "2tt"  ambiguity  which  exists  between  the  measured  phase  dif¬ 


ference  (BM, .)  from  the  actual  antenna  voltages  and  the  theoretical 
phase  difference  (B^)  obtained  from  equation  (1-4).  The  relation  which 
exists  between  these  two  phase  differences  can  be  expressed  as: 


B  =  BM  +  k  ’2TT,  k  =  0,1,2,... 
ij  ij  -  ij 


(2-1) 


UK) 


where  k  has  to  be  determined  and  the  pair  (i))  can  be  either  (12)  or 
(L3).  ihe  actual  measured  phase  difference  can  only  vary  from  - to 
+  1  ;  whereas,  B  can  vary  over  a  much  larger  range  of  values.  Church^ 
described  in  his  paper  a  method  whereby  k  could  be  evaluated  using 
a  dual  interferometer  system.  Since  this  paper  does  not  focus  on 
this  problem,  a  way  is  found  to  avoid  it.  This  method  incorporates 
equations  (1-5)  and  (1-6)  which  tell  exactly  what  B  will  yield  when 
it  is  correctly  averaged  with  respect  to  ■! .  k  is  found  with  a  know¬ 
ledge  of  AVB^  from  equation  (1-5)  or  (1-6)  since  it  represents  the 
multiple  of  2r  needed  to  normalize  AVB_  to  within  the  range  from 
--  to  +-1.  The  "2^"  ambiguity  is,  thus  avoided  for  the  present  simula¬ 
tion. 

The  second  problem  is  that  of  finding  the  proper  branch  cut  of 
the  inverse  tangent  function  in  equation  (1-4)  so  that  the  instan¬ 
taneous  B  will  form  a  continuous  distribution.  If  a  branch  cut  is 
chosen  that  makes  the  B^  distribution  discontinuous  the  correct  value 

of  AVB_  will  not  be  obtaihed  when  the  B^  ' s  are  averaged  with  respect 
(4) 

to  ?.  Allen  developed  techniques  to  properly  determine  the  branch 
cut  when  real  interferometer  data  are  used;  however,  this. problem  can 
also  be  avoided  in  the  present  simulation  by  using  the  calculated 


values  of  AVB^  from  equations  (1-5)  and  (1-6).  The  correct  branch 
cut  can  always  be  found  when  one  knows  AVB .  The  relationship  between 


the  two  can  be  expressed  as: 


=  AVB  -  180' 


(2-2) 


where  n  is  the  branch  cut.  It  should  be  emphasized  that  the  AVB^  used 
in  the  present  simulation  to  avoid  the  "2n"  ambiguity  and  branch  cut 


io : 


problems  Is  not  available  when  one  uses  real  1  n t er t erorne ter  data. 

In  order  that  the  simulation  more  closely  approximate  the  real 
interferometer  system,  random  laciors  were  incorporated  into  the 
simulation  to  account  for  the  perturbations  and  noise.  The  random  factors 
were  all  of  the  form  X  •  RNI)  ( A)  ,  where  RND(A)  is  a  number  taken  from  a 
pseudo  random  number  generator  that  generates  values  from  a  Gaussian 
distribution  with  a  mean  of  zero  and  standard  deviation  of  one.  X  is 
a  scaling  factor.  The  voltage  phasor  of  antenna  1  given  by  equation 
(1-1)  can  then  be  represented  with  a  perturbation  as: 

E  =  [1  +  h . cos  ;  +  T-RND(A)]  +  j  [sin  1  +  T-RND(a)]  (2-3) 

where  the  scaling  factor,  T,  can  be  seen  to  be  the  standard  deviation 
of  the  signal  noise.  The  voltages  for  the  other  antennas  are  modified 
in  a  similar  manner,  Phase  meter  and  system  noise  are  Included  by 
adding  S‘RND(A)  to  B  .  Also  a  perturbation,  PER-RND(A),  is  added  to 


every  o. 


Once  the  B^'s  have  been  averaged  to  yield  AVB^ ,  equations  (1-7) 
and  (1-8)  can  be  used  to  find  the  primary  azimuthal  angle  and  the 
primary  incidence  angle,  The  calculated  angles  of  arrival  are  deter¬ 
mined  sixteen  times  for  each  aperture  size  (D)  in  order  to  obtain  the 
root  mean  square  (RMS)  error  of  the  arrival  angles.  Each  of  the  sixteen 
calculations  would  be  the  same  if  the  same  value  of  RND(A)  were  used. 
However,  different  values  of  RND(A)  are  used  for  each  of  the  calcula- 

1/2 


tions.  The  RMS  error  can  be  expressed  as: 

""16 


RMS  error  = 


^  (Calculated  angle  -  Actual  angle)' 


(2-4) 


16 


The  results  of  this  analysis  are  described  in  Chapter 


programs 


which  were  used  are  listed  in  Appendix  C. 


CHAPTER  3  —  STATISTICAL  ANALYSIS  OF  APERTURE  SIZE 


10  i 


Equation  (2-4)  was  used  to  find  the  RMS  error  of  the  angles  of 
amvai.  However,  the  value  which  was  obtained  can  be  considered  to 
be  only  an  estimate  of  the  RMS  error  based  on  sixteen  samples.  The 
confidence  that  we  have  in  the  estimate  is  related  to  the  size  of  the 
group  on  which  the  estimate  is  based.  Ideally  we  prefer  to  base  our 
estimate  on  a  group  with  a  very  large  number  of  samples. 

Sixteen  samples  were  used  in  the  simulation  analysis  because 
larger  group  sizes  necessitated  the  use  of  an  impractical  amount  of 
digital  computer  time.  This  group  size  was  still  able  to  yield  mean¬ 
ingful  tendencies  in  the  results  of  the  simulation  analysis.  Apply 
ing  statistical  techniques  we  can  determine  the  RMS  error  based  on  a 
very  large  group  size  using  only  a  nominal  amount  of  computer  time.  The 
tools  to  be  used  for  the  statistical  analysis  are  the  standard  devia¬ 
tions  of  both  the  primary  azimuthal  angle  (a  )  and  the  primary  inci- 

“p 

dence  angle  (o^.  ). 

P 

Equation  (1-7)  indicates  that  a  is  a  function  of  AVB^j  AVB^> 
and  y.  The  standard  deviation  of  a  can  then  be  expressed  as: 


c)a 


[(- 


9AVB 


12 


.2  2 

Q  AVB12 


where  a  represents  variance  and 
is  evaluated  with  a  given  set  of 
equations  (1-7)  and  (1-8)  yield: 


)2-%V/2 

(3-1) 

Q  implies  that  the  partial  derivative 
parameters.  In  a  similar  manner 


9a 


+  (■ 


3AVB 


13 


2  2  9q 

)■  +  ( 

Q  AVB13  9 ' 


10-< 


Jo 


,  =  ( — 

0  1  V  3AVB 


12 


,2  2 


Q 


+  ( - E - 

AVB^  ^  jAVB^ 


3  0 

2  2 

Q’  +  (j' 


.  2  2,1/2 


(3-21 

Since  >  is  a  constant,  e  2  is  zero. 

The  partial  derivatives  needed  to  evaluate  equations  (3-1)  and 
t3-2)  can  be  expressed  as: 


3a 


AVB 


13 


3  AVB 


12 


(AVBl2r  *  sin  y  * 


(3-3) 


3a 


3  AVB 


13 


AVB^  '  sin  y  ’  X^ 


(3-4) 


where 


AVB  -  AVB  .  •  cos  Y  j 

\  -  1  +  <  AVB12  T  8in-Y— } 
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AVB ^2  -  cos  Y  *  ^^13 


3AVB12  (N2-TT2-sin2Y“X2  -  X2)1/2 


(3-5) 


3  o 


AVB 


13 


cos  Y  ’  AVB 


12 


3  AVB 


13 


,..2  2  2  v  v2. 1/2 

(N  -it  -sin  y*X2  -  X2) 


(3-6) 


where  X2  =  (AVB^)^  +  (AVB^)2  -  2* cos  y 


AVBi2  •  AVB13 


The  variances  of  AVB.^  (°^yB  )  an^  AVB^  ^°AVB  ^  Can  ex_ 

pressed  in  similar  forms.  Therefore,  the  use  of  AVB^,  where  j  equals 
2  or  3,  will  considerably  simplify  the  work  to  follow.  Since  AVB^  can 
be  represented  as: 


«*"RI«wpps 9PW|pV«|||i9Pi| 


10'. 


AVB, 


B11  +  Ul^  +  •  •  ■  +  1S1  j 

180 


180 


its  variance  can  be  written  as: 
180 


•i 

AVB 


\ 


3  AVB 


,  ( 


il 


ij  i=i 


sb 


ij 


>2  V 


ij 


When  the  partial  derivatives  are  evaluated  the  result  is: 

180 

1  ^  2 

?,  o„i 


a  2 

AVB 


lj  (180) 2  i^l  Blj 


(3-7) 


In  order  to  find  o  i  we  write  equation  (1-4)  in  a  simplified 

Bij 


form: 


Ri  _  -1  rElR  E1I  E1R  EHi 

j  ^  1 p  p  4-  F  F 

lj  E1R  EjR  +  E1I  Ejl 


(3-8) 


Using  equation  (3-8)  a  i  can  be  written  as: 

lj 


3b: 


°Bi  =  (3E 

BU  3eir 


>2  4  +  (• 

Q  IF 


SB 


3E 


ii 


II 


2  2  >*1, 

Q  E1I  9EjR 


)2  o\  + 

Q  jR 


3  B 


ii 


3E 


jl 


v 2  2 

}  °E 

Q  ji 


(3-9) 


Equation  (2-3)  shows  that  the  variance  of  the  real  or  imaginary 

2 

part  of  the  phasor  voltage  is  equal  to  T  .  Since  this  is  valid  for 
all  three  antenna  voltages,  we  have: 


2  2  2 

aE'  =  °E  =  °E 

1R  II  2R 


l  2  2  t2 

?  =  Op.  =  oF  =  T 

J2I  n3R  31 


(3-10) 


Substitution  of  the  partial  derivatives  and  equation  (3-10)  into 


equation  (3-9)  gives 


o  i  =  T  ' 

BU 


<E?i  +  e?r)(4+eii)2  + 

(eir  ejr  +  EU  Ej!  +  E1R  +  EJr  Eii) 


(3-11) 


Equations  (3-1)  through  (3-7)  with  equation  (3-11)  can  be  combined 
to  give  the  standard  deviations  of  the  primary  angles  of  arrival  as 
a  function  of  T,  the  signal  noise.  The  computer  programs  which  applied 
this  analysis  are  listed  in  Appendix  D.  The  results  of  these  programs 


are  described  in  Chapter  4. 
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CllAPTKR  4  —  RESULTS  AND  CONCLUSION.-. 

The  number  of  parameters  used  in  the  present  analysis  necessitated 
that  some  should  be  assigned  a  constant  value,  which  is  typical  of  all 
the  values  in  its  range,  and  others  should  be  varied  over  a  range  of 
values.  The  primary  incidence  angle  ('•>),  the  secondary  angles  of 
arrival  (a  ,  0  )  f  gamma,  the  frequency  oi  the  incident  waves  (f),  and 
Che  perturbation  (PER)  of  3>  are  all  assigned  typical  values.  Although 
is  varied,  its  effect  on  the  bi-ari  'g  error  -  aperture  size  relation¬ 
ship  is  not  investigated  in  depth. 

The  dependent  variables  for  the  pr^ent  analysis  are  the  bearing 
errors  (the  errors  in  a  and  9  ).  The  two  major  independent  variables 
are  the  aperture  size  (D)  and  the  standard  deviation  of  the  signal 
(!)  .  The  units  of  T  are  volts.  The  minor  independent  variables  are 
the  wave  interference  factor  (h)  and  the  standard  deviation  of  the  phase 
neter  and  system  noise  (S)  .  The  parameter  h  can  vary  from  zero  to  one 
and  represents  the  ratio  of  the  amplitudes  of  the  secondary  incident 
wave  to  the  primary  incident  wave.  The  parameter  S  has  been  described 
as  the  phase  meter  noise;  however,  in  the  simulation  analysis  it  was 
chosen  to  account  for  all  of  the  noise  produced  by  the  electronic 
equipment  which  processed  the  signals.  The  units  of  S  are  radians - 

The  data  obtained  fVom  the  statistical  analysis  and  the  simulation 
analysis  are  listed  in  Appendix  A  and  Appendix  B,  respectively .  Figures 
2  through  13  represent  these  data  in  pictcri  1  form;  however,  tne  repre¬ 
sentation  is  somewhat  distorted.  Figures  2  through  13  are  traces  of 
actual  digital  computer  ou-liue  printer  output.  The  distortion  arises 
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Figure  9  RMS  Incidence  Angle  Error  vs  T  —  Param 
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Figure  13  RMS  Incidence  Angle  Error  vs  D  —  Parameter  h 
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because  an  on-line  prinCer  can  only  plot  discrete  values  oi  a  variable. 
The  distortion  Introduced  by  this  discrete  plotting  process,  however, 
does  not  affect  our  ability  to  recognize  general  tendencies  in  the  data. 
If  a  more  critical  analysis  is  required  reference  should  be  made  to  the 
precise  tables  in  Appendix  A  and  B. 

For  a  given  set  of  parameters  it  is  evident  from  Figures  2  through 
13  that  the  azimuthal  angle  error  and  the  incidence  angle  error  have 
the  same  relative  magnitude.  It  should  be  pointed  out  that  this 
condition  results  from  the  present  choice  of  parameters.  It  could  have 
been  easily  shown,  for  example,  that  if  the  incidence  angle  was  near 
90°  the  magnitude  of  the  incidence  angle  error  would  have  been  larger 
than  the  azimuthal  angle  error. 

Figures  2  through  5  show  the  bearing  errors  versus  D  at  five  values 
of  T.  Figures  2  and  3  represent  the  statistical  results  and  Figures 
4  and  5  represent  the  simulation  results.  It  is  seen  from  these  figures 
that  the  simulation  results  approximate  the  statistical  results  as  they 
should.  At  h  =  .7  and  T  =  .5  volts  the  simulation  analysis  shows  an 
azimuthal  error  of  approximately  1°  for  an  aperture  size  of  100  meters 
(1A  at  3  Me.).  For  the  same  set  of  parameters  this  is  exactly  what  the 
statistical  approach  predicts. 

Let  us  assume  that  the  maximum  bearing  error  that  we  will  tolerate 
is  1°.  At  T  =  .5  volts  Figures  2  and  3  indicate  that  the  aperture  size 
must  always  exceed  100  meters  (1 A ) .  This  result  is  general  though  for 
it  implies  that  for  any  frequency  the  aperture  size  must  exceed  IX. 

If  the  signal  noise  is  now  decreased  to  T  =  .1  volts  the  aperture  must 
now  exceed  about  33  meters  or  1/3  A.  The  previous  description  is  only 


1 . 1 


one  Interpretation  of  the  data.  The  Interpretation  of  the  results 
will  depend  on  particular  needs. 

Figures  6  through  9  display  the  same  data  shown  in  Figures  2 
through  5  in  a  different  form.  These  figures  present  the  bearing 
errors  as  a  function  of  T  for  different  values  of  D.  It  is  easily 
seen  from  these  figures  that  the  bearing  errors  vary  linearly  with  T. 

As  D  is  increased  the  bearing  errors  approach  zero  as  is  demonstrated 
by  the  decreasing  slopes  of  the  plotted  lines. 

Figures  10  through  13  show  the  effects  of  h  and  S  on  the  bearing 
errors.  Figures  10  and  11  show  that  at  h  =  0  and  D  =  5  meters  (1/20 
the  bearing  errors  are  about  2.5°  with  S  =  .1  radian.  This  represents 
a  relatively  small  effect  for  such  a  substantial  noise.  With  S  =  .1 
radian  Figures  12  and  13  show  there  is  no  effect  on  the  bearing  error 
when  h  is  changed  from  0  to  .7,  When  h  is  increased  to  «_ b'-oct  maximum, 
h  =  .99,  the  error  only  changes  by  1°.  Thus,  wave  interference  is 
seen  to  have  a  relatively  small  effect  on  the  bearing  errorB  if  the 
differential  phase  is  averaged. 

a  was  chosen  to  be  0°  for  most  of  the  calculations  in  this 
P 

chapter.  This  value  was  arbitrarily  chosen  to  represent  the  complete 
range  of  a^’s.  Tables  5  and  6  show  the  variation  of  bearing  errors 
with  D  for  different  values  of  a  .  The  small  variation  of  bearing 
error  with  ct^  suggests  that  the  results  for  =  0°  may  be  extended  to 

other  values  of  a  . 

P 

In  conclusion  the  resultB  obtained  in  the  present  analysis  sub¬ 
stantiate  our  intuitive  feelings  about  the  response  of  the  bearing  errors 
to  changes  in  D  and  T.  As  D  is  decreased  the  errors  increase  and  as  T 


u. 

increases  the  errors  increase.  The  precise  variation  is  shown  in  Figures 
2  through  13.  The  agreement  between  the  simulation  and  statistical 
approaches  to  the  problem  strengthens  our  belief  that  this  analysis  is 
valid.  One  area  which  could  he  considered  ior  further  studies  by 
simulation  would  be  to  find  the  effects  of  those  parameters,  which  were 
left  constant  in  this  report,  on  die  bearing  errors. 

It  is  important  to  note  that  the  results  shown  here  are  derived 
from  a  particular  model  of  the  interferometer  RDF  system.  Thus,  a 
further  experimental  study  is  needed  to  validate  the  model  used.  An 
investigation  of  errors  in  determining  angles  of  arrival  for  a  small 
aperture  size  is  suggested. 
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INI*  ERF  EKUMET  EK  1 N  ST RUMEN  TAT  I  UN  FOR  PROPAGATION  RESEARCH 

This  paper  outlines  some  oi  the  111’  propagation  capabilities  of  the 
L)F  interferometer  instrumentation  now  in  development  or  in  use  at  South¬ 
west  Research. 

this  review  emphasizes  the  versatility  for  propagation  research 
provided  by  the  interferometer  class  of  instruments  both  in  amplitude 
and  phase  measuring  versions.  We  will  briefly  describe  some  propagation 
work  in  progress  (ur  soon  to  be  in  progress)  which  supports  the  Army- 
sponsored  single  site  location  program. 

The  amplitude  interferometer  techniques  in  use  today  have  their 
origins  in  the  Adcocks  and  spaced  loops  of  the  period  before  World  War 
II.  These  early  amplitude  difference  measurements  provided  an  indirect 
means  to  obtain  phase  measurements.  The  Adcocks , be ing  a  relatively 
easy  development , came  at  an  early  date;  the  spaced  loop,  being  a  more 
difficult  development,  came  much  later. 

There  are  several  references  to  interferometers  in  the  1930’s. 
Eckersley  in  1938  was  probably  the  first  to  use  spaced  loops  in  research 
problems  which  required  an  interferometric  approach.  Eckersley  measured 
direction  of  arrival  and  estimated  visibility  functions  using  the  spaced 
loop  as  a  Michelson  interferometer, 

After  World  War  II,  phase  measuring  interferometer  developments 
moved  toward  multiple  wavelength  baselines  as  are  now  in  use.  The  de¬ 
velopment  of  rapid  phase  measuring  techniques  in  the  1930’s  and  small 
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general.  purpose  ilii’iiaL  coiii|miutb  in  C lDUO's  lias  produce^  rapid  ad¬ 
vancement  in  111-’  interlcromeLer  applications  lor  direction  finding  and 
voiy  recently  lor  single,  sice  location. 

Figure  1  aL  the  top  summarizes  ellects  imposed  by  the  ionospheric 
medium  on  a  target  signal  which  is  illustrated  by  the  box  on  the  left, 
ihe  Larget  signal  passes  through  a  propagating  medium  which  establishes 
a  mode  structure,  alters  the  polarization, and  causes  lateral  deviation, 
as  illustrated  by  the  box  in  Lhe  center.  At  the  rlgliL,  we  state  some 
major  measurements  on  the  emerging  signal  which  are  important  fur  propa¬ 
gation  research  and  radiolocation  technology. 

From  the  point  of  view  of  propagation  research,  measurements  pro¬ 
viding  directional  mode  resolution,  angular  width,  and  specular  to  dif¬ 
fracted  ratio  and  others  can  be  interpreted  in  terms  of  the  ionospheric 
irregularities  which  cause  them.  These  measurements  yield  data  on  the 
details  of  irregularities,  effective  ionospheric  tilts,  etc.,  —  all 
of  importance  for  understanding  die  propagation  effects. 

At  the  same  time,  these  measurements  are  significant  for  radio¬ 
location  technology.  Directional  mode  resoluLic'n  permits  reduction -of 
wave  interference  UF  error  while  knowledge  of  source  angular  width  anu 
specular  to  diffracted  ratio  determines  optimum  system  design  including 
aperture 

The  interferometer  as  a  direction  finder  provides  a  uniquely  ver¬ 
satile  system  for  measuring  these  cliarac  teristics .  Specular  to  dif¬ 
fracted  ratios  are  obtained  from  phase  angle  distribution  measurements. 
Angular  width  is  obtained  by  fringe  visibility  measurements.  Directional 


mode  resolution  can  be  obtained  by  either  F  and  A  polygon  displays  or 


Four  ici'  L runs  term  ot  measured  umpi  i  t  ude  and  phase  distributions  provided 
suliiciont  daLa  samp  Los  arc  available.  At  the  same  Lime,  a  combination 
ol  long  and  short  baseline  phase  and  amplitude  measuring  interferometers 
permits  separation  ot  polurixuL  ion  1  rom  wave  interference  effects  while- 
observing  boLli.  We  will  specifically  indicate  this  capability  by  ex¬ 
amples,  but  first  we  will  review  the  instrumentation. 

the  wide  aperture  phase  measuring  interferometer  (Figure  2)  at 
SwRi  has  been  previously  described .  the  equipment  includes  two  data 
readout  systems,  the  digital  system  shown  here  and  an  analog  system 
based  on  a  University  of  Illinois  development  and  shown  in  the  next 
figure . 

The  analog  arrangement  (Figure  3)  permits  directional  resolution 
of  multiple  modes  by  providing  sum  and  difference  phase  meter  envelope 
displays  in  the  form  of  polygons  obtained  with  a  storage  oscilloscope. 
Typically  the  polygon  displays  are  developed  when  observations  are  made 
over  intervals  on  the  order  of  15  to  30  seconds.  Two  component  signals 
provide  parallelograms,  three  provide  hexagons,  etc.  The  phase  angles 
are  measured  by  the  slopes  of  the  sides  of  the  parallelograms.  These 
phase  angles  may  be  converted  into  bearing  computations  either  by  feed¬ 
ing  the  data  directly  into  the  computer  in  a  semiautomatic  manner  or 
by  the  use  o f  an  analog  calibrator  (piiase  angles  are  slopes  of  displays). 

The  four  baselines  of  the  interferometer  are  displayed  in  four 
quadrants  on  the  CRT,  two  short  baselines  at  the  top,  two  long 
at  the  bottom.  A  grating  alidade  is  provided  to  transfer  envelope 
slopes  to  the  digital  computer  where  bearings  are  computed  after  four 
phase  angles  have  been  measured  in  the  proper  sequence.  Alternatively, 
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an  analog  calibrator  which  generates  electronic  sLrobe  lines  superim¬ 
posed  over  (.he  stored  polygons  can  be  used  Lu  direclly  obtain  Lhe  azi- 
rnuth  and  elevation  oi  any  one  mode.  The  calibration  is  shown  in  the 
lower  part  oi  Llie  iigure  where  separate  azimuLh  and  elevation  controls 
are  indicated.  The  calibrator  shown  is  based  on  a  development  at  the 
University  oi  Illinois  with  a  number  of  changes  . 

in  use  the  procedure  is  to  apply  the  strobes  first  to  the  short 
baseline  displays  to  obtain  rougii  cuts  and  then  to  refine  the  controls 
with  a  second  adjustment  witli  the  strobes  superimposed  on  the  long  base¬ 
line  displays.  This  procedure  eliminates  the  phase  ambiguities  and  per¬ 
mits  the  bearing  to  be  obtained  without  trial  and  error  searching  for 
the  proper  values  since,  in  general,  the  short  baselines  will  provide 
only  one  answer.  This  method  of  obtaining  bearings  is  completely  apart 
from  the  fully  automatic  digital  system  which  uses  a  conventional  coin¬ 
cidence  type  phase  meter  and  computes  azimuth  and  elevation  of  the  strong¬ 
est  mode  from  phase  angles. 

Figure  4  shows  a  two  component  field  near  8,5  MHz.  The  amplitudes 
of  the  two  components  are  evident  from  the  length  of  the  sides  of  the 
display . 

Figure  5  shows  the  result  of  applying  plane  wave  editing  to  an 
observed  two  component  signal.  It  will  be  noted  that  the  parallelogram 
display  for  at  least  one  baseline  is  replaced  with  an  hour  glass  figure 
or  so-called  bundle  of  sticks  figure  which  corresponds  to  a  number  of 
instances  where  equal  amplitudes  were  recognized  at  the  ends  of  the 
baseline.  Some  early  tesc  results  with  plane  wave  editing  have  been 
provided  in  one  of  the  ocher  papers. 
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spaced  loi-p  aiitL'iina  may  be  unduLSLood  l  rum  many  points  a)  view 
including  that  ol  a  roLaling  baseline  air.pl  i  Lude  interferometer  Figure 
t>  .shows  ilia  external  view  ol  a  rotating  coaxial  spaced  loop  anLenna  da- 
signed  tot  -tUU  Kl'M  rotation.  this  is  an  Army  development  for  tactical 
app 1 icut ions. 

Figure  7  shows  Lite  spaced  loop  indicator  and  receiver  used  in  Lite 
tilm  which  follows  in  a  tew  minutes. 

Figure  8  shows  calculated  spaced  loop  response  patterns  fur  a 
single1  signal  as  a  function  of  polarization  at  an  elevation  angle  of 
.  The  characteristic  spaced  loop  response  patterns  permit  continu¬ 
ous  qualitative  information  on  the*  incident  polarization.  Thinking  of 
the  antenna  as  a  spaced  array  with  loop  elements,  the  array  null  (or 
interferometer  null)  remains  fixed  on  the  bearing  independent  of  polari¬ 
zation  while  the  element  null  (or  loop  null.'  rotates  and  blurs  accord¬ 
ing  to  the  incident  polarization,  A  blurred  loop  null  means  elliptical 
polarization.  A  sharp  loop  null  indicates  linear  polarization  and  ro¬ 
tates  with  the  orientation  angle  of  the  incident  polarization  ellipse. 
Vertical  polarization  is  on  the  left,  horizontal  on  the  right  and  cir¬ 
cular  at  bottom  center.  Equivalent  simple  loop  patterns  are  shown  in 
some  of  the  patterns  as  dotted  lines.  In  the  film  which  follows,  all 
of  these  patterns  can  be  found  in  varying  degrees  in  combination  with 
wave  interference  effects.  We  have  not  attempted  to  show  here  theoreti¬ 
cal  patterns  for  the  wave  interference  effects. 

As  an  introduction  to  rotating  spaced  loop  data,  we  want  to  show 
a  sequence  of  film  taken  on  two  typical  skywave  signals,  Varying  polari¬ 
zation  will  be  seen  on  both  signals  as  the  film  progresses,  while  the 
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observed  bearing  remains  nearly  iixed  depending  on  yiaw.  interference- 
cond i t Lons 

I'ho  first  sequent  u  shows  WWV  uL  lU  Mliz  wiinh  lias  a  range  ol  up- 
proximnteiv  8\>U  miles.  this  sequence  ol  lilm  shows  relatively  stable 
elliptical  po  la  r  i  ;:u  L  ions  as  indicated  by  the  depth  of  Liu  simple  loop  or 
element  nulls,  The  photographs  were  taken  wiLli  a  sequence  camera  at 
two  frames  per  sec  >nd  and  are  being  shown  at  eigul  1  ramus  per  second  so 
that  there  is  a  factor  ol  four  increase  in  the  rate  ol  display  At  a 
number  of  points  in  the  film,  a  sudden  change  in  Liu-  shape  of  the  pat¬ 
tern  accompanied  by  a  loss  in  symmetry  will  occur-  These  patten,  ro- 
vide  sense  displays  obtained  by  unbalancing  the  spaced  loop  (wave  inter¬ 
ference  effects  present). 

The  next  sequence  will  shew  KLC,  an  interrupted  carrier  sta¬ 
tion  at  8.5  MHz,  located  near  Uojstoo,  Texas  This  station  has  a  range 
of  approximately  215  miles.  In  contrast  to  the  previous  sequence  of 
elliptical  polarization,  this  target  shows  rotating  linear  polarization 
from  vertical  through  horizontal  sometimes  in  the  clockwise,  other  times 
reversing  to  the  counterclockwise  direction.  It  is  important  to  note 
that  a  qualitative  measurement  of  the  state  of  polarization  is  available 
at  ail  times  making  possible  a  comparison  between  the  polarization  in¬ 
formation  and  any  wave  interference  effects  which  might  be  present. 

One  could  also  compare  other  direction  finders  simultaneously  to 
obtain  a  measure  of  the  degree  to  which  either  polarization  error  or 
wave  interference  error  are  present.  It  is  also  important  to  note  that 
the  wave  in r e r f e rence  or  scattering  error  reduction  effect  of  the  spaced 
loop  would  provide  another  means  for  separating  polarization  error  from 
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wave  intericronce  error  .such  as  in  the  comparison  oJ  an  /.dcock  and  a 
spa  'ed  loop  or  a  simple  loop  and  a  spaced  loop. 

We  now  wish  to  show  a  comparison  oi  spaced  loop  patterns  with  simul¬ 
taneous  sum  and  difference  phase  meter  displays  from  the  wide  aperture 
phase  measuring  interferometer.  During  this  sequence,  the  film  will  be 
stopped  to  show  displays  frame  by  frame.  One  can  observe  the  phase 
angles  and  consequent  directions  of  arrival  of  the  components  in  the 
incident  field,  the  polarization,  and  direction  of  arrival  of  the  re¬ 
sultant.  One  can  also  see  the  number  of  components  present  and  the 
signal  amplitude  of  each. 

This  is  a  sequence  of  film  obtained  on  WWV  at  10  MHz,  The  sequence 
shows  the  rotating  spaced  loop  pattern  on  the  right  and  sum  and  differ¬ 
ence  phase  meter  displays  from  the  crossed  baseline  phase  interferometer 
on  the  left.  The  parallelograms  on  the  far  left  are  developed  from  the 
phase  angles  from  the  long  baselines.  The  parallelograms  (which  are 
nearly  straight  lines)  from  the  short  baselines  are  just  to  the  right 
of  the  long  baselines.  This  sequence  of  data  shows  elliptical  polari¬ 
zations  with  rotating  orientation  angle  as  indicated  by  the  spaced  loop 
response  patterns.  The  phase  parallelograms  for  each  baseline  indicate 
one  mode  much  stronger  than  the  other  as  determined  by  the  length  of 
the  parallelogram  sides.  The  long  baseline  parallelograms  at  far  left 
are  stable  and  show  greater  resolution  than  the  much  narrower  short 
baseline  parallelograms  on  the  right  as  expected  by  the  increased  reso¬ 
lution  of  the  long  baselines.  Tue  polarization  sequences  rapidly  change 
from  the  elliptical  polarizations  towards  linear.  However,  similar  to 
the  previous  WWV  sequences,  one  does  not  observe  much  rotation  of  the 
orientation  angle  (simple  loop  null  position). 


Finally,  il  is  impo'.Lanl  to  note  iliaL  in  the  iuily  automatic  mode 
the  digital  compuLor  output  is  displayed  similar  to  a  manner  illustrated 
by  i he  examples  shown  in  Figure  1.  This  output  provides  azimuth  and 
elevation  angle  automatically  and  repeLiLively  and  optionally  a  variety 
of  oilier  parameters  which  can  be  derived  from  Lhese  including  the  SSL 
position,  ionospheric  tilt  angle  and  direction,  and  all  kinds  of  av¬ 
eraging  statistics.  in  the  present  display,  Lhe  azimuth  is  on  the  left, 
elevation  next,  followed  by  number  of  baseline  phase  sequences  and  Lhe 
number  stored.  The  plots  to  Lhe  rigiit  were  put  in  by  band  to  show  typical 
one-minute  sequence  tor  uhree  target  ranges.. 

in  this  example,  three  time  series  displays,  each  about  one-minute 
long,  are  provided  for  three  targets  at  different  ranges:  at  the  top, 
WWV;  in  the  center.  K.LC,  the  interrupted  carrier  signal  215  miles  to 
the  east ; and  at  the  bottom,  our  own  target  at  New  Berlin,  Texas,  34 
miles  east.  Relatively  stable  azimuth  data  are  obtained  in  the  low 
elevation  angle  case  whereas  relatively  stable  elevation  angle  data 
are  obtained  for  the  very  high  elevation  angle  case.  The  scale  on  the 
azimuth  plot  is  10°;  the  scale  on  the  elevation  plot  is  2° 

Figure  10  shows  one  more  example  —  a  measured  phase  angle  dis¬ 
tribution  along  both  interferometer  baselines.  These  distributions 
can  be  interpreted  in  terms  of  specular  to  diffracted  ratios.  Our 
data  show  variations  from  .1  to  1. 

We  emphasize  in  closing  that  the  capability  of  the  com-ptrtje-f  in¬ 
strumented  and  analog  readout  interferometer  systems  can  now  be  fully 
exploited  for  propagation  research.  Basic  rigorous  interference  tech¬ 
niques  are  available  for  obtaining  new  quantitative  data  for  radio¬ 


location  technology  and  propagation  research. 
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Long  Range  Target,  WWV.10  MHz,  Range  858  miles 


95 


OS 


Azimuth 


Medium  Range  Target,  KLC,  Houston,  Texas 
Range  215  miles 


Short  Range  Target,  New  Berlin,  Texas 
Range  34  miles 

FIGURE  9 
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lONOCKAMS* 

A  study  has  recently  been  conducted  ij y  liainum  ol  the  Stanford 
University  Fleet  conics  Laboratory  (SFL)  wiLb  FM/CW  ionograms  taken  over 
a  1900  km  path  l rum  Texas  to  Staniord,  for  which  an  absolute  time  delay 
scale  was  available  for  each  ic'nogram.  ^  ^  As  a  test  of  the  procedures 
developed  by  A vco  for  calculating  absolute  time  delay  scales  for  oblique 
incidence  ionograms  from  the  inter -mode  time  delays  evident  in  the  data, 
the  method  was  applied  to  the  SFL  Texas  to  Stanford  ionograms,  and  the 
resulting  time  scale  compared  with  the  known  time  scale.  The  Avco  com¬ 
puter  generated  time  delay  versus  virtual  height  plots  were  used  in  this 
analysis;  the  Texas  to  Stanford  ploL  is  shown  in  Figure  1  The  lowest 
order  pure  F  mode  for  this  path  should  be  the  IF  mode.  For  the  iono¬ 
grams  where  no  E  mode  was  present,  the  time  scale  was  determined  by  using 
the  relative  time  delays  between  the  pure  F  modes. 

The  time  delay  scales  determined  by  using  the  computer  plots  without 
regard  to  the  time  delay  axis  of  Barnum’s  records  were  usually  within 
0.03  milliseconds  of  the  absolute  delay  time  scales  given  on  the  data, 
which  had  been  determined  by  two  calibrated  cesium  clocks  For  one  case 
the  time  scales  ditfered  by  0,06  milliseconds,  and  this  was  the  maximum 
deviation  obtained  for  the  ionograms  scaled. 


*Work  Sponsored  by  RADC  under  Contract  F30602-69-C-0096 


1 1'.  Barn  uni's  work  the  If  mode  Iran1  was  used  to  estimate  a  modi- 1 
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lc  was  our  purpose  to  see  ii  the  reflectrix  method  ^  could  be 
used  directly  to  generate  mode  traces,  i.e.,  to  go  from  the  measured 
IF  trace  to  a  reflectrix,  and  then  to  whatever  other  mode  trace  is 
desired,  such  as  the  2F,  without  any  iterative  procedure  or  costly  ray 
tracing.  The  Avco  reflectrix  program  converts  an  oblique  incidence 
ionogram  into  an  equivalent  midpath  vertical  incidence  ionogram  and 
then  plots  virtual  height  versus  ground  range,  parametric  in  frequency, 
for  various  takeoff  angles.  The  program  also  allows  one  to  superimpose 
a  function  of  virtual  height  and  ground  range  onto  the  reflectrix  plots 
for  the  purpose  of  calculating  ranges  to  critical  reflection  regions. 

The  points  of  intersection  of  the  reflectrix  curve  and  this  function 
are  printed  out  for  each  frequency. 

One  of  the  oblique  incidence  ionograms  from  Barnum’s  report  (see 
Figure  2)  was  used  to  generate  a  2F  mode  from  the  actual  IF  mode.  First 
the  cime  delays  versus  frequency  for  the  IF  mode  were  read  off  the  iono¬ 
gram.  These  were  inverted  into  the  equivalent  midpath  vertical  inci¬ 
dence  ionogram  shown  in  Figure  3  and  then  transformed  into  the  re¬ 
flectrix  curves  shown  in  Figure  A.  The  total  path  range  and  the  ground 
range  to  the  first  reflection  point  of  a  2F  mode  were  then  read  in  as 
data.  The  program  determines  the  points  of  intersection  of  this  vertical 
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Figure  3.  Equivalent  Midpath  Vertical  Incidence  Ionogram  Inverted  from 
Mode  Identified  Ionogram. 
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line,  with  the  lu 1  feet i ices ,  and  calculates  the  ray  path  distance  ii'ur. 
the  transmitter  to  the  point  ol  reileotion.  From  this,  by  assuming 
equal  and  symmetrical  hop  cunt  igurut  ions ,  1l  is  possible  U  obtain  the 
time  delay  tor  a  given  frequency  lor  the  21  mode  simply  by  multiplying 
the  distance  computed  above  by  4/c  tor  each  frequency.  figure  b  shows 
a  tracing  ol  the  original  ionogram  with  Lhe  21'  trace  points  calculated 
in  this  fashion  from  the  11'  data.  Whiie  tlie  calculated  times  are  not 
precisely  on  the  Line  drawn  along  the  bottom  of  the  2 F  trace  of  Lne 
actual  ionogram,  they  easily  fall  within  the  spread  of  the  experimental 
trace  at  the  lower  frequencies  (i.e.,  they  are  all  within  0  02  milli¬ 
seconds  of  this  line).  The  biggest  discrepancy  (up  to  0,06  milliseconds) 
occurs  at  the  upper  frequency  end  of  the  2F  mode,  where  a  discontinuity 
in  the  original  data  trace  seems  to  indicate  the  presence  of  a  more 
complicated  layer  shape  or  irregular  region. 
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1  .  HALKOKULNU  ALL  IN'  IKUULL'l  lUN 

l'he  liK'DK  work  at  Clw rgia  l'ecli,  which  has  been  .sponsored  by  the 
L.  S.  Army  Licet  routes  Command,  is  directed  toward  reducing  bearing 
errors  in  small-aperture  simple  loop  UK  direction-finders,  especially 
when  used  for  Lonospherieally- propagated  signals  Also,  as  a  result  of 
the  error  reduction  work,  techniques  have  been  developed  to  improve  sen¬ 
sitivity  and  display  characteristics.  l'he  capability  of  obtaining  a 
measure  of  elevation  angle-of-arrivai  with  a  simple  loop  has  been 
developed . 

A  major  consideration  is  the  tactical  nature  of  simple  loop  DF 
systems  requiring  that  all  the  techniques  be  amenable  to  tactical  usage. 
Techniques  must  be  simple  and  operate  in  real-time  without  the  benefit 
of  ancillary  data  processing  or  computer  devices.  Weight,  volume,  and 
power  consumption  must  be  minimized;  hence,  all  circuitry  is  solid 
state  except  the  receiver  and  displays,  which  are  GFE, 

The  techniques  were  developed,  implemented  in  a  complete  exploratory 
HF/DF  system,  and  operationally  evaluated  in  a  realistic  tactical  en¬ 
vironment  using  both  cooperative  and  uncooperative  transmissions  from 
a  wide  variety  of  locations  and  distances. 

The  following  will  describe:  (1)  the  sources  of  bearing  errors  in 
loop  L)F  systems;  (2)  the  fundamental  approaches  lo  error  reduction; 

(3)  a  general  description  of  the  specific  techniques  for  error  reduction; 
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1  i  .  SOURCKS  OF  Hl.AK INC  JiRROKS 

Large  bearing  errors  are  usual  when  conventional  111-'  loop  direction- 
finding  (.DF)  systems  are  used  for  ionospherically-propagated  signals, 
lho  inaccuracies  are  primarily  determined  by  the  characteristics  of  the 
incident  signal  which,  unfortunately,  change  constantly  as  a  function  of 
time  LI  i  —  L  t>  J  •  The  ionosphere  not  only  affects  signal  parameters  such 
as  amplitude,  polarization,  and  phase,  but  also  introduces  multipath 
propagation  modes  L7]--[10j.  Real-time,  instantaneous  bearing  error  is 
a  function  of  Che  real-time,  instantaneous  state  of  the  incident  signal; 
hence,  as  the  signal  characteristics  vary,  bearing  errors  vary  [  12  J- [ 1 A  J  > 
However,  at  certain  times  the  signal  state  is  such  that  minimum  bearing 
errors  can  be  expected. 

There  are  two  major  error-producing  effects  in  loop  DF  systems: 
signal  depolarization  caused  by  Faraday  rotation  and  wave  interference 
caused  by  multipath  propagation  modes. 

Polarization  Lrrors 

Loop  DF  accuracy  depends  upon  the  existence  of  a  correct  and  con¬ 
stant  relationship  between  the  direction  of  propagation  and  the  polari¬ 
zation  of  the  b-field  vector.  The  correct  polarization  occurs  when  the 


incident  If—  field  contains  no  horizontal  components  perpendicular  to  the 
direction  of  propagation.  in  this  case,  which  generally  occurs  with 


IV) 


ground-wave  signals,  no  polarization  error  exisLs,  and  the  measured  bear¬ 
ing  will  be  Llie  true  great-circle  bearing  (GDI)).  However,  any  E-field 
polarization  Lilt  creating  horizontal  components  perpendicular  to  tiie 
direction  ol  propagation  produces  bearing  error  lor  sky-wave  signals  [2J. 

A  signal  propagated  via  the  ionosphere  is  depolarized  into  a  time- 
variant,  elliptically  polarized  lorm  [11)]  1  oiarination  tilt  introduces 

voltages  in  the  Loop  that  are  not  zero  when  the  plane  of  the  loop  is  per¬ 
pendicular  to  the  GCli.  As  the  loop  is  rotated,  the  resulting  null  in¬ 
dicates  a  bearing  which  is  incorrect  ,  Such  bearing  errors  are  a  function 
of  the  elevation  augle-o£-arrival,  1,  and  polarization  tilt  angle,  1. 

The  general  expression  for  the  total  induced  voltage,  V,  in  a  small, 
vertical  loop  (diameter  less  than  about  0.1>.)  located  in  freespace  is 
given  by 

V  =  KE  [sin  6  cos  1  -  cos  c  sin  4>  sin  1  ]  (1) 

where 

2  7T 

K  =  - —  x  (area  of  loop)  x  (number  of  turns), 

0  =  azimuth  angle  between  direction  of  propagation  and  the 
normal  to  the  plane  of  the  loop, 

A  =  wavelength, 

T  =  polarization  tilt  angle  or  the  angle  between  the  E-field 
vector  and  the  vertical  plane  normal  to  the  plane  of  the 
loop, 

$  =  elevation  angle  between  direction  of  propagation  and  the 
normal  to  the  plane  of  the  loop, 

and  E  =  the  electric  field  strength  in  volts/meter. 

If  0  is  not  zero  when  a  minimum  value  of  V  is  obtained,  a  bearing 
error  results.  It  can  be  shown  from  (1)  that  the  bearing  error,  0,  is 
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given  by 

0  -  Lan  (s  i  11  ;  Lan  :  )  (2) 

when  the  loop  is  rotated  unLil  (he  receiver  output  is  zero.  figure  1 

depicts  b  .is  i  [unction  o!  :  lor  various  selected  values  oi  1  and  clearly 

0 

shows  that  any  nonzero  value  oi  i  produces  a  bearing  error  with  !  •  0  . 

c  o 

A  maximum  error  ol  10  occurs  wiLh  >  =  10  lor  ail  elevation  angles-of- 
arrival  o  is  a  function  oi  both  I  and  hence,  bearing  accuracy  can 
vary  as  the  elevation  angle-o i-arr ival  changes  even  though  polarization 
tilt  is  constant..  The  time  variations  of  bearing  errors  caused  by 
polarization  rotation  are  a  function  of  elevation  anglc-of-arrival  and 
the  polarization  rotation  rate,  if  the  tilt  angle  r  varies  linearly  as 
a  function  of  time,  (2)  becomes 


0(C) 


tan  (sin  1= 


tan  Kt) 


(3) 


where  t  is  time  in  seconds  and  K  is  the  time-rate-of-change  of  i.  The 
time  derivative  of  (3)  is 


dl  (t) 
dt 


K 


sin  v 


2  2  2 

cos  Kt  +  sin  1  sin  Kt 


(4) 


This  function  is  plotted  in  Figure  2  for  selected  elevation  angles- 
of-arrival.  Figures  1  and  2  indicate  that  continuous,  uniform  polariza¬ 
tion  rotation  which  occurs  during  polarization  fading  causes  continuous 
changes  in  bearing  error  and  rate-of-change  and  that  these  variations  are 

determined  by  both  <t  and  T .  Each  elevation  angle  produces  different  0 
,  dO 

and  - —  curves. 
dt 

0  O 

The  time-rate-of-change  of  0  is  smallest  when  T  is  near  0  or  180  , 
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which  co rresponds  lu  minimum  bearing  ct  rur  polarization,  and  largest 
when  i  is  near  d ’U  or  90‘  which  coi  i  e s ponds  lu  maximum  bearing  error 
polar  izat  ion  For  low  elevation  angle  signals  (0  1  -  U'j  )  the 

variations  in  ~  are  very  pronoun. ed  as  ■  rotates  over  J60‘  ;  signals 
arriving  at  me  higher  elevations  angles  produce  less  change  in  j-  as 
i  rotates  over  JoU'  For  the  case  I  90',  tin-  bearing  relation  rare 
is  constant  implying  that  the  '  rue  bearing  ‘.annul  bt  defined  using  -j- 
as  a  measure 

Polarization  tilt  also  decreases  the  sharpness  or  distinctness  of 
tiie  bearing  null  by  decreasing  die  t  ime-rate-o  f -change  of  the  induced 
voltage  as  the  antenna  rotates  through  the  null  The  change  in  sharp¬ 
ness  is  much  greater  for  the  lower  elevation  angle  signals  Higher 
elevation  angle  signals  (60°  -  9U  )  have  less  than  a  15  per  cent  decrease  in 
the  time-rate-of-change  of  induced  voltage. 

Polarization  errors  are  particularly  prevalent  when  the  sky-wave 
signal  is  experiencing  Faracay  rotation  in  the  ionosphere  and  when 
multipath  effects  are  minimal  This  normally  occurs  near  the  MUF,  es¬ 
pecially  on  a  one-hop  path,  Faraday  rotation  leads  to  time-decorrelated 
fading  between  the  vertical  and  horizontal  components  of  the  incident 
signal  with  topical  fade  rates  in  the  order  of  several  per  minute. 

To  summarize,  bearing. polarization  rotation  errors  can  be  as 
large  as  +  90°  from  the  GCB.  Inaccurate  bearings  are  evidenced  by 
swinging  bearings  and  display  broadening  with  the  slowest  swing  and 
minimum  broadening  near  the  GC13,  However,  as  the  elevation  angle  in¬ 
creases,  the  swing  rate  and  bearing  display  characteristics  become  more 
uniform  over  a  complete  rotation  cycle 


1  od 

.1  vu  -  Ln  le  r  t  e  rence  (.Mul  Lip.ilti)  In  rs 

l’o  lar  ujL  iun  errors  are  an  inlierenl  inadequacy  ul  a  loop  antenna 
and  •..,in  be  reduced  by  tiding  an  entirely  diiierent  DF  antenna  system, 
sueii  as  spaced  loops  1 1 0  J  or  an  Adcock  [l’/J,  However,  wave  interference 
errors  have  their  sources  entirely  external  to  the  loop  and  are  a  source 
of  bearing  error  Lor  all  DF  systems.  The  Loop  is  more  susceptible  to 
wave  interference  (or  multipath)  errors  because  i :  is  a  small-aperture 
device. 

In  practice,  more  tlian  a  single  wave  front  at  tlie  same  frequency 
is  incident  on  a  DF  antenna.  In  fact,  most  of  t he  energy  is  normally 
contained  witliin  discrete  modes  tliat  arrive  from  many  liorizontal  and 
vertical  angles,  The  results  of  investigations  using  wide-aperture , 
special-purpose  antennas  (sucli  as  tlie  Wullenweber  array)  have  delineated 
the  complex,  multimode  structure  of  an  ionospherically-propagated  signal 
[3],  These  studies  have  indicated  that  the  discrete  modes  in  the  inci¬ 
dent  multimode  signal  are  not  only  often  of  comparable  amplitude,  but 
also  arrive  at  different  bearing  angles.  Also,  the  relative  phases  be¬ 
tween  the  modes  are  time-variables.  The  resultant  of  all  these  varia¬ 
tions  is  a  corrugated  phase  front  whose  pattern  is  changing  as  a  func¬ 
tion  of  time.  With  a  small-aperture  loop  antenna,  tlie  observed  bearings 
at  any  instant  will  be  normal  to  the  piiase  surface  at  that  point. 
Therefore,  under  multipath  effects,  the  bearings  are  constantly  sliifting 
and  swinging.  A  small-aperture  loop  follows  every  corrugation  and  shift 
of  the  equiphase  surface  since  the  aperture  size  is  too  small  for  sig¬ 
nificant  averaging  across  corrugations.  Also,  large  (+  90°)  bearing 
error  can  occur  due  to  multipath  propagation  even  thougli  no  horizontal 
i omponents  exist  in  the  Incident  field, 
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Multi  pa  tli  propagation  produces  ,»  variable  amplitude  (fading  I  signal , 

* 

and  there  exists  .1  relationship  between  the  amplitude  variation  and  bear¬ 
ing  swing  rate  and  etret  Maximum  hearing  error  and  swing  rate  or.ur  at 

> 

the  lade  null;  min ituum  bearing  error  and  swing  rate  occur  on  the  fade- 
crest 

in  general,  multipath  propagation  creates  time-correlated  fading 
between  the  h-tield  vertical  and  horizontal  components  whereas  polariza¬ 
tion  rotation  produces  t ime -decorrelated  lading  However  in  Loth  cases, 
minimum  bearing  error  and  swing  rate  occur  on  the  maximum  of  the  vertical 
E-field  component  of  the  fading  signal. 

Multipath  fading  is  the  more  common  type;  however,  the  effects  of 
polarization  fading  are  worse  for  a  loop  DE  system  since  a  rotating 
bearing  is  displayed  with  little  or  no  pattern  change 

III.  FUNDAMENTAL  APPROACH  TO  ERROR  REDUCTION  [18-21] 

The  fundamental  approach  to  error  reduction  was  based  on  the  deter¬ 
mination  of  the  specific  characteristics  of  the  incident  signal  that  pro¬ 
duce  minimum  bearing  errors  and  the  development  of  techniques  for  using 
these  characteristics  to  either  (1)  restrict  or  unblank  the  display  of 
bearing  information  to  only  these  conditions  that  produce  minimum  bear¬ 
ing  error,  or  (2)  create  distinct  and  interpretable  differences  between 
accurate  and  erroneous  bearing  displays 

I V .  SPECIFIC  APPROACHES 

For  both  multipath  and  polarization  errors,  minimum  bearing  errors 
occur  when  the  vertical  E-field  component  is  at  a  maximum,  The  fade 


crest  detection  method  oi  error  reduction  uses  a ‘ dual-channel  receiver 


>i. vent  i can  1  lo  jj)  u!  ;  it.*.-  oilier  uses  a  vertical 


One  iiannel  l  -  used  a-> 
probe  Lo  n  ili<  magnitude  .'I  >  Ik-  '/*.- i  i  i  *  t*  1  *  ou.p  neiii  and  supply  dis- 

plav  unb  l  JuK  l  Hi-,  pulse  •  at  at  near  Lite-  lade  •  tests 

l-'-i  b-'.li  unit  ip  on  ana  po  lot  i  it.it  ion  triors,  minimum  bearing  errors 


Oc'.L'l  wile 

u  »  'n« 

i>-  it  i  iijj 

,  :iW  i 

ng  i 

lie  i 

s  i  ta  i  u  i  uu  in  I  in-  swing  rat 

c  disci,  imi 

nut  ion  ie 

t  ii*'d 

t.1  i 

tia* 

be  a 

■  inti 

swing  rate  md  unblanks  the 

d  i  s  ]j  i  a  y 

as  a  l  .-tic 

l  U‘ii 

ol  t  Ik  i 

i  iic* 

talc 

-el  v- 

li.inge  -i  t lie  -tearing  swing 

I  h  i  s 

t  echni quo 

>.  in 

btj  list'd 

w  i  i  h 

a  *, 

lu^lc 

-receiving,  channel 

t'onveni  i  ’uu  1  loop  [)|-  systems  display  all  bearings  and  the  operator 
rejects  ot  a-oept-.  hearings  bused  *>n  subje-.iive  evaluations;  l  he  fade 
*  rest  Jetts  '  i  >n  md  swing  rate  discrimination  Ji-n  -ibed  above  remove  the 
more  ei  ron--  ."is  i'eatines  I  ro:'  *  he  operator's  vie”*'  and  display  only  the 
no  it1  i-  *  11  r  1 1  e  be  a  i  i  ngs 

1  he  tilted  loop  technique  lot  bearing  e- rot  -eduction  does  not  use 
unblinking  bto  *•  re.i  t  >.*«  distinct  and  i  u*  rpi  o  t  able  differences  between 
accurate  and  erroneous  beat  i tig  displays  Also,  the  technique  can  be 
used  to  obtain  a  measure  ot  elevation  angle  This  technique  can  also 
oe  used  with  a  single-re  eiving  channel 

V.  FUNCTIONAL  BLOCK  DIAGRAM 

Figure  3  is  a  functional  block  diagram  ot  the  final  version  of  the 
txpl  tv  system.  (Ihe  swing  rate  d  i  sc  t  imina  •  -r  i  not  shown,  having 

been  replaced  by  the  tilted  le-’p  <e  linique  which  per  lottos  essentially 
tiie  same  function,,'  Figure  '>  is  a  front  panel  view  of  the  post-IF  signa 
processor  which  includes  till  the  post-IF  circuitry  except  the  display 
and  sine/rosine  potentiometer 


Figure  3.  Simplified  Functional  Diagram  of  Post-IF  Circuitry 


two  antennas  are  used-  a  i  op  -'a  It 


1(j  j 

a i eat  fur  Channel  1  the 
KACAL  b  J  )  j  h ,  and  a  "eiiieal  ptobc  for  Channel  J  figure  a  is  a  lane1  lit  na  1 
diagram  ot  the  ant  emu  assembly;  iiguie  0  is  an  overall  view  Hie  Joup 
is  meehan i ea 1 1 y  lotJted,  scan  iaie  is  e anL innous 1 y  variable  C r cm  u  tc 
JjO  rpm  the  loop  plane  t_an  be  ecnl  inuonsl  y  tilled  from  0J  t.o  9  0U  ,  about 
a  horizontal  axis.  lilt  iate  Is  u  per  seeond  The  vertitcl  probe  is 
9  inches  long  and  is  mounted  directly  in  the  center  ol  the  loop  ihe 
effective  heights  ot  the  1  top  and  probe  were  experimentally  matched  at 
1U  MlTz  to  avoid  excessive  c  t  oss- coupl  i  ng  problems 

The  loop  IF  signal  contains  the  scan  modulation  which  can  be  either 
envelope  or  synchronously  detected. 

Convent ionel ,  s ing le -channe 1  loop  DF  systems  use  incoherent  deter  • 
cion  of  the  loop  scan  modulation  by  an  envelope  detector  Envelope  de¬ 
tection  has  several  disadvantages  for  DF  applications-  One,  the  output 
signal  represents  the  magnitude  of  the  scan  modulation  and  is,  hence, 
unipolar  with  a  high  harmonic  content.  The  low-pass  filtering  follow¬ 
ing  the  envelope  detector  must  have  a  frequency  response  much  greater 
than  the  basic  scan  fundamental  in  order  to  pass  h igher-order  harmonies 
and,  hence,  preserve  the  integrity  of  the  bearing  nulls-  This  wide 
bandwidth  adds  noise,.  Two,  the  envelope  detecto”  further  deteriorates 
low  SNR  signals,  e  g  ,  an  envelope  detector  input  SNR  of  -10  dB  yields 
an  output  SNR  of  about  - li  dii  Three,  the  envelope  detector  contains  a 
180c  ambiguity  in  the  bearing  information  in  that  a  bearing  null  at  0^ , 
for  example,  cannot  be  distinguished  from  a  true  null  at  6^+  r.  Four, 
the  envelope  detector  has  a  limited  dynamic  range  and  saturation  effects 
degrade  the  bearing  display. 


VERTICAL  PROBE  ANTENNA 


HORIZONTAL 
CENTER  LINE 


Figure  5.  Basic  Electrical  and  Mechanical  Configuration  of  Til  table  Loop 
Antenna  System. 


Synchronous  (coherent )  detection  reduces  the  disadvantages  oi 
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velope  t  Lurolie  renl )  detection  and  was  a  natuiai  addition  Lo  the  sysLem 
due  co  the  presence  ot  the  channel  2  ouLput  ,  This  outpuL  is  coherent 
(or  can  be  made  coherent)  wiLn  tlie  loop  channel  output  except  for  the 
desired  scan  moduiation  containing  die  bearing  information  The  chan¬ 
nel  1  output  is  used  as  a  coherent  demodulation  signal  for  the  synchro¬ 
nous  detector. 

The  synchronous  detection  output  is  a  sinusoid,  bipolar  signal  with 
the  zero-crossings  containing  the  bearing  information.  Resultant  ad¬ 
vantages  over  the  envelope  detection  are:  (1)  narrow-band  post-detection 
filtering  can  be  used  to  improve  SNR  since  the  synchronous  detector  out¬ 
put  is  a  simple  sinusoid  at  the  antenna  scan  fundamental;  (2)  low  SNR 
signals  are  not  degraded;  (3)  sense  information  is  inherent  in  the  bi¬ 
polar  signal;  and  (4)  dynamic  range  is  improved, 

In  the  exploratory  system,  synchronous  detection  improved  the  sen¬ 
sitivity  by  about  10  dB.  Figure  7  illustrates  tiie  improvement  in  sensi¬ 
tivity  and  display  characteristics  compared  to  envelope  detection. 

The  IF  output  of  channel  2  is  modulated  by  fading.  The  fading  en¬ 
velope  is  detected  by  a  conventional  envelope  detector  and  low-pass  filter 
circuitry  and  applied  to  the  fade  crest  detectors,  where  a  positive  gat¬ 
ing  pulse  is  generated  at  the  fade  crest.  For  optimum  detection,  two 
types  of  fade  crest  detectors  may  be  selected  depending  upon  the  type  of 
fading  (polarization  or  multipath)  encountered.  One  type,  the  averaging 
and  comparison  technique,  was  designed  primarily  for  polarization  fading 
conditions  where  a  good  signai-to-noise  ratio  and  reasonably  periodic 
fade  cycles  exist.  The  other,  the  peak  sample  and  hold  technique,  was 
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dcs  igned  ptimniiLy  lor  ii.ul  t  tpjt  h  l  ad  inn  'undiliuus  wilh  J  signal  lc-vi-l-, 
and  irregular,  er*  itu,  nanper i od ic  amp  1  i  L  ode  vat  ilium-; 

The  channel  J  do  l  t  ed  signal  is  also  supplied  l.o  clutter  re-due  l  ion 
l  i  r  c  u  i  L  r  y  that  1  e  moves  display  iIuUli-  duo  Lo  on-oil  transitions  ol  the 
o  a  r  r  i  e  r 

fhe  cireuitry  uses  i  he  output  ol  channel  2  as  a  gating  signal  Lo 
allow  only  i.he  true  bearing  nulls  to  be  displayed  and  gate-out  the  false 
bearing  nulls  due  to  carrier  interruptions.,  Channel  2  output  does  not 
contain  the  Loop  scan  modulation  but  does  contain  the  on-off  modulation; 
channel  1  output  contains  both  scan  and  modulation  nulls.  Clutter  re¬ 
duction  is  based  on  the  simple  concept  that  true  bearing  nulls  occur 
only  when  the  channel  2  output  is  high  and  the  channel  1  output  is  si¬ 
multaneously  low,  i„e„,  in  a  null.  The  clutter  reduction  circuitry  ap¬ 
plies  deflection  voltages  to  the  beating  display  only  when  the  ■"liannel  2 
output  is  above  an  established  threshold  and,  hence  carrier  on-off 
transitions  are  removed  from  the  display  leaving  only  the  true  bearing 
null  Figure  8  illustrates  the  effects  of  clutter  reduction. 

In  the  Loop  channel,  the  pos t -detection  circuitry  (display  shaping 
and  deflection  generator)  processes  the  bearing  information  for  display 
on  the  storage  scope  (HP141A  with  variable  persistence).  The  display 
drive  signals  are  applied  through  gate  circuits  prior  to  the  scope! 
Conventional  sensing  techniques  require  display  deflection  plate  switch¬ 
ing;  however  the  synchronous  detection  sensing  technique  obviates  the 
need  for  deflection  plate  switching  allowing  for  more  flexibility  in  the 


selection  of  displays. 
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VI.  l)i:i'All_Kl)  LH.St'K  1  I’  1  ION  01-  l  kKulx  KLDI  L  i  luM  1  bLUA  1‘,'LT.j 

tail.-  llLSt  He  III  l  ion 

file  J  Vi:  r  ag  L  ng  mb  i"!'|iot  iioii  mode  (figuic  ')  )  op-  rat  c*-  by  (I)  aver¬ 
aging  tlu*  lading  enve  L  .’{>*.•  ,  and  (J)  omparihg  tin.  average  value  with  u 
weighted  envelape  An  ouLput  •.»  l'us  only  i!  the  weighted  input  exceeds 
the  average  value  [he  weighting  lacier  i.alled  the  lade  targin,  deter¬ 
mines  how  tat  abo"e  t  In;  avetage  value  the  input  iu"st  r  i  s  *  before  an  out¬ 
put  jams.  For  example,  with  a  1  dll  lade  margin  an  output  occurs 
when  tlie  input  is  1  dB  oi  gt  eater  above  the  avetage. 

An  averaging  time  of  approximately  lb  seconds  was  chosen  based  on 
an  assumed  minimum  rate  at  four  fade  cycles  per  minute  The  averaging 
tine  muse  be  at  least  equ  ;1  to  rbo  fade  period  since  it  requires  at 
least  one  lade  cycle  to  indicate  short- term  behavior 

The  fade  margins  of  0,  1,  1.  b  and  10  dB  were  based  on  the  antici¬ 

pated  range  of  the  maxima  relative  to  the  overage  for  both  multipath  and 
polarization  fading.-  Multipath  may  be  character ized  by  a  Rician  distri¬ 
bution  [10]  which  allows  for  the  finite  probability  of  having  a  fade 
crest  of  infinite  amplitude  However,  10  dB  was  chosen  as  a  practical 
upper  limit  since  the  percentage  of  time  that  this  level  is  exceeded  for 
Rician  statistics  is  less  than  0^01%.  Based  on  theoretical  ('.onsideva- 
tions.  1  dB  was  considered  to  be  about  the  highest  fade  margin  that  could 
be  reliably  used  for  polarization  fading  (Operational  evaluation 
proved  this  to  be  true.)  Figures  10  and  11  illustrate  the  fade  detector 
input  (bottom)  and  output  (top)  waveforms  for  fade  margins  of  0  to  5  dB, 
respectively  The  time  base  was  2  sec/cm  (The  5  dB  was  later  changed 
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ratio  is  i  nc  r  eas  i  ng :  nub  1  atik  !  ng  dviog  reus  t  ng  sign.-)1  -  1  ■.'"noise  ratio 
is  rad'i  t'i  Tes  t  s  si:  -wed  thu'  I'jiU  n  t  displays  to  periods  of  in- 
reading  s i gnu  1  - 1 o - in  i se  ratios  reduied  tin-  number  of  spurious  bearings 
and  in'  reused  •••-<'  v  »  especially  when  !  he  maximum  si gnal ~ ' o-no ise 
rat  io  i-  1  .w  Opera'  ie”  '  •>  dep:<  1  ed  in  Figure  It  with  -a  t  ini"  base  of 
l  s*-«  /  m  It  may  be  no* ed  that  I  be  unblanking  time  is  proportional  to 
too  magnitude  <  f  the  fade  ires'  amplitude  Th<>  1  a1' go r  the  fade  crest 
level  'be  longer  unb I  unking  occurs  providing  automat  Jo  amplitude  weight 
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Figure  13.  Peak  Sample  and  Hold  Fade  Cr.  « 
0.3  v/ cm)  and  Output  (Top  True 
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Inc  ;>RU,  Kiguie  1  ,  uub1  nRy  'in.-  display  di  p«.-nd  i  ng  on  '.iic  tine 
t  j'.o-i.'l-ciuugi-  cl  Luc  bearing  swings  Inc  L»<  arings  ct  iMiivc-rtc.1  to 
digital  form  L  e  a  pulse  jc.iirij  at  ee  l.  bearing  null  As  '.lie  bear¬ 
ing  swings,  L  tie  Lime  interval  between  'be  b>-.n  iug  pulses  varies;  a  sLeady 
bedring  prodm.es  a  constant  time  interval  between  pulses  (about  100  msec 
tor  a  d  rps  loop  scan )  When  a  jiredet  eriuined  number  (eight  in  the  ex¬ 
ploratory  svstein)  ot  consecutive  bearing  azimuths  occurs  within  a  preset 
number  of  degrees  lro m  one  another,  tile  Ski)  unblanks  'he  display  'ltie 

preset  allowable  swing  region  can  be  varied  from  about  5"  to  1'0:  ’I he 

SRI)  converts  the  bearing  pulses,  wbii.h  are  in  the  form  ot  a  pulse-position 
modulated  signal,  to  a  DC  analog  form  with  the  DC  level  indicating  the 
bearing  in  degrees  relative  to  reference  azimuth  The  bearing  pulses 
sample  a  ramp  pulse  which  also  starts  a  timer  The  DC  level  cov -espond- 
iug  to  tne  first  bearing  pulse  is  stored  and  compared  with  the  DC  level 
of  seven  succeeding  bearing  pulses.  If  no  comparisons  are  out  cf  *-hr 
voltage  range  corresponding  to  the  allowable  bearing  swing  rate,  the 
timer  completes  its  2  6  second  cycle  and  tiie  next  bearing  display  is  un¬ 
blanked  If  excessive  swing  occurs,  the  timer  and  sampling  circuits  are 
reset  to  start  another  comp-'  ■  i son  sequence 

V  HI.  T 1 LT  ED  _  L  OOP  _  TDCiiN  IQUE 

A  major  problem  with  loop  J)F  systems  is  that  the  bearing  null  in¬ 
formation  and  consequently  the  display  for  both  accurate  and  inaccurate 
bearings  remain  essentially  unchanged  --  except  for  bearing  shift  --  as 
polarization  changes  lienee,  the  ability  to  separate  accurate  from 
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inaccurate  bearings  based  on  display  characteristics  docs  not  exist  in 
amvt'UL ional  systems  The  exploratory  system  developed  by  Georgia  lech 
uses  a  LiiLed  loop  technique  that  crcaLes  distinct  and  int erpre table  dif¬ 
ferences  between  accurate  and  erroneous  bearing  displays.  Also,  the 
technique  can  be  used  to  obLain  a  measure  of  the  elevation  angle  of  tne 
incident  sky-wave  signal  the  to! lowing  describes  how  Lliis  is  accomplished. 

The  expression  oi  the  tree-space  induced  voltage,  E^ ,  in  a  small 
loop  rotating  about  a  vertical  axis  is  given  by 


cos  T  sin  6  -  sin  1  sin  4>  cos 


where  K  =  a  constant  that  is  a  function  of  loop  area,  frequency,  and 
number  of  turns, 

’»  =  the  incident  signal  polarization  tilt  angle  relative  to  ver¬ 
tical, 

9  =  the  angle  between  the  direction  of  propagation  and  tne  normal 
to  the  plane  of  the  loop, 

and  f  =  the  elevation  angle  of  the  incident  signal  referenced  to  the 
horizontal . 

Equation  (5)  can  be  simplified  to 


E  =  KA  [cos  (6  -  6)]  (6) 

L  o 

where  A  = 
o  ■ 

A^  =  cos  Y 

=  -sin  T  sin  <t 

and  6  =  tan 

Equation  (5)  shows  that  the  response  is  a  function  of  both  signal 


depolnr  l /at  ion  ami  elevu . ion  angle,  and  (b)  Indicates  Llial  the  response 
as  a  function  ol  it  is  a  sinusoidal,  iigure-eighL  response  for  all  com¬ 
binations  ot  t  and  !  Al curate  bearing  indications  are  obtained  only 
whenever  K  =  U  and  d  =  U°  or  18U°  it  can  be  seen  from  (5)  that  this 
condition  occurs  only  when  t  -  L)°  or  18(J°  Then 


K  -  K  sin 


(/) 


For  nonzero  values  of  i  and  1,  (b)  indicates  two  inaccurate  but.  sym¬ 
metrical  nulls  180°  apart.  The  null  locations,  however,  do  not  occur 
at  the  correct  hearing  along  the  direction  of  signal  propagation.  Un¬ 
fortunately,  these  inaccurate  bearings  exhibit  the  same  shape  and  null 
symmetry  as  accurate  bearings,  Figure  15a  illustrates  the  type  of 
bearing  display  obtained  from  a  vertical  loop.  The  tips  of  the  lobes 
correspond  to  the  null  locations.  A  double-nulled,  symmetrical  response 
>ccurs  for  all  signal  conditions,  but  the  bearing  display  rotates  as  the 
polarization  rotates  "nd  will  indicate  an  accurate  bearing  only  when 
1  O'  for  sky-wave  signals 

Tilting  is  performed  by  positioning  the  plane  of  the  loop  at  a 
tilt  angle  a  relative  to  horizontal  while  maintaining  rotation  about 
the  vertical  axis.  Figure  16  depicts  the  relationship  between  a  and 
l  The  tilt  angle  u  is  variable  from  0C  t/o  90"1 ,  The  free  space  in¬ 
duced  voltage,  E  ,  for  the  tilted  loop  is 

L 


E  =  K  jcos  T  sin  6  cos  a  -  sin  V  sin  cos  0  cos  a 


+  sin  ?  cos  T  sin  o 


which  can  be  reduced  to 


Figure  15.  Bearing  Displays  Illustrating  (a)  Double-Nulled  Symmetrical, 

(b)  Double-Nulled  Asymmetrical,  and  (c)  Single-Nulled  Responses 
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Equation  (S)  shows  that  the  tilted  loop  response  is  also  a  function  of 
:  and  but  (9)  indicates  that  the  response  is  a  figure-eight  pattern 
providing  accurate  bearings  only  when  the  incident  sky-wave  signal  is 
vertically  polarized  ('r  0°).  If  i  /  0s,  the  response  is  not  a  figure^ 

eight  pattern  due  to  a  nonzero  value.  The  response  nulls  are  sepa¬ 
rated  by  less  than  180°,  and  the  bearing  display  exhibits  asymmetrical 
patterns  indicating  inaccurate  bearing  information.  Hence-,  accurate 
bearings  can  be  separated  from  erroneous  bearings  using  the  symmetry 
of  the  bearing  display. 

The  effect  of  tilting  on  display  characteristics  is  illustrated  by 
Figure  15.  If  a  is  much  less  than  4>,  the  display  is  as  shown  in  Figure 
15a  with  the  bearing  rotating  as  T  rotates.  In  (9),  is  much  less 
than  Bq  and  the  response  is,  for  all  practical  purposes,  a  symmetrical 
figure-eight  pattern  for  any  V  value.  As  a  is  increased  and  approaches 
1 ,  the  display  continues  to  rotate  with  T  rotation  but  develops  discernible 
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as  yuui.  el  i  y  lor  ii.uirero  values  ol  i  as  illustrated  by  Figure  1  i.b .  With 
i  t  lT,  B  in  ('■))  increases  as  ■  increases  and  begins  In  uffei  t  Liie  re¬ 
sponse  as  i  •  I.  Since  B^  is  a  1  uni.  Lion  oi  *,  the  degree  of  symmetry  is 
a  timet  ion  oi  t.  When  i  90“  ,  asyiiuuetry  is  largest.  When  :  0",  Lucre 

is  symmetry,  Lhus  indicating  an  accurate  bearing,  (The  above  is  true 
ior  any  tilt  angl<_.)  since  i  ,  B^  is  always  less  than  B^>  and  tile- 
response  is  double  nulled.  however,  with  .i  ;  note  from  (8)  that  one 
oi  the  response  nulls  remains  fixed  while  the  oLher  null  rotates  with 
and  the  two  nulls  merge  to  form  a  single  null  when  i  =■  90°  as  illus¬ 
trated  by  Figure  15c.  Here  =  B^  and  a  eardioid  response  results. 
Single-nulled  responses  occur  only  when  a  _  •!>.  It  lias  been  experi¬ 
mentally  determined  that  the  best  tilt  angle  tor  the  collection  of  bear¬ 
ing  information  is  with  «  4>, 

The  transition  from  a  double-nulled  to  a  single-nulled  response  as 
i  is  increased  also  allows  for  a  measure  of  elevation  angle-of-arrival 
For  t  i ,  a  double-nulled  response  occurs  for  any  f  value.  For  a 
a  single-nulled  response  occurs  as  T  rotates  through  90° .  Hence,  0  can 
be  estimated  by  increasing  >  until  a  single-nulled  response  is  obtained 
at  •  =  90"  .  In  practice,  most  ionospherically-propagated  signals  ex¬ 
perience  depolarization  as  large  as  90° ■  In  fact,  many  signals  undergo 
cyclic  polarization  rotation.  Therefore,  an  estimate  of  elevation  angle 
of  arrival  can  be  obtained  on  many  ionospherically-propagated  signals 

Another  method  of  analyzing  tne  tilted  limp  response  is  to  resolve 
the  tilted  loop  into  two  components  -  one  vertical  and  one  horizontal 
The  horizontal  component  produces  the  third  term  (B^)  in  (8) .  The 
first  and  second  terms  in  (8)  are  due  to  the  vertical  component.  Signal 
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pickup  in  Lhe  her  izonLul  I'ompoueiil,  whicu  occurs  only  when  1  Is  nonzero, 
combines  with  Lhe  vertical  oompone n L  pickup  to  produce  a  response  wiin 
asymmetrical  nulls.  When  i  0"  and  tiie  incident  signal  is  vertically 
polarized,  the  horizontal  component  Signac  pickup  is  zero,  and  symmetri¬ 
cal  nulls  occur  defining  an  accurate  bearing  situation.  The  horizontal 
and  vertical  components  result  tram  Lhe  .same  physical  loop  and,  lienee, 
have  the  same  constant  K  in  the  response  equations.  If  separate  loops 
were  used  tot  the  horizontal  and  vertical  components,  the  constants  in 
the  response  equations  may  not  Lrack  as  a  function  of  frequency, 

1 X .  OPERATIONAL  EVALUATION  -  TAPE  CREST  DETECTION 

Tes t  Methods 

Field  tests  were  conducted  over  a  five-month  period  during  the 
winter  of  1966-67  „  All  tests  were  performed  on  snort-wave  broadcast 
stations  of  known  beating  over  propagation  paths  ranging  from  600  km  to 
F600  km.  Essentially  all  the  tests  were  taken  on  transmissions  at  or 
below  the  calculated  F2MUF.  The  nominal  test  period  for  each  test  was 
five  minutes. 

Three-hundred  and  eighty- four  tests  were  performed  using  sixty- 
seven  different  frequencies  in  the  5-20  MHz  range  originating  from 
twenty-three  separate  locations.  Two-hundred  and  thirty-six  of  the 
tests  were  performed  during  the  1100-1500  GMT  (dawn)  and  2200-0200  GMT 
Husk)  periods  of  relatively  high  ionospheric  variation. 

The  test  phase  included  an  operational  comparison  of  the  explora¬ 
tory  systems  performance  with  that  of  the  AN/PRD-5  HF/DF  system 


DaLa  Reduction 


liu'  biMi  t  ug  .1  iia  I  t  urn  cadi  test  were  reduced  to  the  following  biiiiib- 
tical  parameters  : 

l  he  average  heating  deviation  f  tom  llie  great-circle-  bearing 
(GCJi > .  A  positive  value  indicates  a  clockwise  deviation;  a 
negative  value,  a  counterclockwise  deviation 

Si'  ■  i’ne  total  t/i:>iuLlial  spieail  oi  Lite  bearing  data. 

SD  ;  Standard  deviation 

For  malysis  purposes,  the  reduced  data  from  individual  tests  were 
grouped  and  turther  processed  to  obtain  the  ioilowing  parameters: 

A'  :  Tlie  average  ot  the  averaged  bearing  deviations, 

A"  •  i'he  average  ol  Llie  magnitudes  of  A",  ,  A1''  L 

SP  '  The  average  ol  (.be  spreads 

SD  :  The  average  ot  the  standard  deviations, 

RMS  Error;  Defined  as ~\J  (SD)^  + 

Comparison  oi  Fade  Crest  Detector  Modes 

The  initial  operational  tests  compared  Che  averaging  and  compari¬ 
son  technique  (Mode  1)  with  the  peak  sample  and  hold  technique  (Mode  2), 
Comparisons  were  made  by  using  both  modes  in  sequence  on  the  same  trans¬ 
mission.  Two  series  of  tests  were  performed  for  a  wide  variety  of  trans¬ 
mission  characteristics  with  the  test  results  presented  in  Table  I 
Mode  1  test  results  are  for  fade  margins  providing  smallest  A0,  Table  I 
shows  that:  (1)  the  average  bearing  deviation  and  standard  deviation 
are  about  10'  tor  both  modes;  (2)  the  total  spreads  are  restricted  to 
less  than  one-half  of  a  quadrant;  (3)  Llie  peak  sample  and  hold  technique 
performs  best  in  the  presence  of  a  low  SNR;  and  (A)  the  performance  of 
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both  modes  is  essentially  llie  same  I  or  both  .series  ol  LesLs  even  motion 
ditlerenl  transmissions  and  Limes  were  used.  Alsu,  the  results  show  Li; 
peak  sample  and  hold  technique  perlormance  iur  a  low  SNR  is  essentially 
identical  to  that  tor  a  good  SNR,implving  t ha L  this  mode  operates  inue- 
pendently  of  the  average  value  ot  the  signal. 


The  effect  of  fade  margin  increase  on  performance  was  also  evalua¬ 
ted.  Tests  were  performed  on  thirty-six  transmissions  having  good  SNR 
and  essentially  periodic  lading  using  fade  margins  of  0,  3,  6,  and  10 
dli  in  sequence  for  each  test.  The  reduced  data  are  shown  in  figure  17. 

TABLii  1 

Comparison  of  Fade  Detector  Modes 


Test  Series  Wo. 

.  Tests 

Mode 

Total  Wo. 
Bearings 

of 

!mI 

deg 

SP 

deg 

SD 

deg 

1 

15 

1 

456 

9 

33 

9 

15 

2 

635 

11 

44 

13 

2 

38* 

1 

1733 

7 

32 
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38* 

2 

1813 

12 

42 

10 
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Figure  17.  Test  Results  Depicting  Effects  of  Fade  Margin  on  Bearing  Error 
and  Dispersion. 


1 V) 


healings  wi'ii'  displayed  at  0  and  3  dli  Jade  max')', ins  on  ail  tests ; 
however,  Jive  tests  and  twenLy  lest;.  produced  no  hearings  at  iade  mar¬ 
gins  ot  a  and  Id  dli,  respectively.  lion  e,  ii  appears  tiiat  a  i  dJi  lade 
margin  is  about  the  lurgcsL  value  LiiaL  shouLd  be  used  in  any  pi  veil  opera¬ 
tional  situation,  especially  it  the  observation  time  is  limited.  Figure 
17  shows  monoLonic  decreases  in  bearing  tpread,  standard  devia'ion,  and 
RMS  error  with  increasing  lade  margin  and  supports  the  contention  lhaL 
the  more  accurate  bearings  occur  at  Lite  crests. 

Comparisons  between  the  3  and  b  dll  iade  margins  were  furLher  in¬ 
vestigated  by  an  additional  series  oi  twenty-two  tents  which  produced 
essentially  the  same  results  as  depicted  in  Figure  17  for  these  two 
fade  margins. 


Overall  Exploratory  Systems  Test  Results,  Using  Optimum 


Fade  Crest  Detector 


A  series  of  tests  was  performed  using  the  optimum  type  of  fade 
crest  detector  for  the  particular  fading  conditions.  Mode  1  was  se¬ 
lected  for  the  transmission  with  polarization  fading  and  a  good  SNK; 

Mode  2  was  selected  for  the  transmissions  with  multipath  fading  and/or 
a  poor  SNR.  The  results  are  depicted  in  Figure  18. 

The  Mode  I.  data  were  obtained  from  eight  separate  transmitter  lo¬ 
cations  in  the  North  American  continent;  however,  a  majority  of  the  tests 
were  on  VOA  (Greenville,  N.C.)  transmission  near  the  F2MUF.  Most  of 
the  Mode  2  data  were  obtained  using  transmissions  from  Europe,  Africa, 
and  South  America,  North-South  and  East-West  propagation  paths  were 


present  during  dawn,  dusk,  and  daytime  ionospheric  conditions.  Therefore, 
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Overall  Test  Results  Using  Optimum  Fade  Crest  Detector. 
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the  test  results  represent  a  wide  variety  oi  signal  and  propagation  lud- 
ui  t ions . 

i'lie  results  tor  boLii  modes  are  very  similar  and  show  average  bear¬ 
ing  errors  ol  L(JU  or  iess  with  a  majority  of  the  bearings  restricted  to 
within  +  lb'1  of  Lhe  CCii  by  the  unblanking. 

Data  from  seven  tests  wiLh  heavy  Kl-'l  and/or  very  low  SNR  were  in¬ 
cluded  in  the  overall  Mode  2  data.  'i'hese  data  were  also  reduced  sepa¬ 
rately  and  are  shown  under  the  Mode  2  data.  As  expected,  these  data 
exhibit  larger  dispersion  and  iess  accuracy  than  the  tests  with  beLter 
signal  conditions;  however,  these  data  indicate  that  it  is  possible  to 
obtain  fairly  accurate  bearing  information  using  Mode  2  even  for  very 
poor  signal  conditions. 

X.  COMPARISON  OF  EXPLORATORY  SYSTEM  WITH  THE  AN/PRD-5  HF/bP  SYSTEM 

Fade  Crest  Detection 

The  exploratory  system  (ES)  and  the  AN/PRD-5  system  were  compared 
by  simultaneous,  side-by-side  testing. 

The  AN/PRD-5  operates  from  0.5-20  MHz  using  a  two-foot  loop  rotat¬ 
ing  at  1  rps .  The  display  unit  is  a  3-inch  CRT  with  P7  phospher  that 
displays  a  completely  symmetrical  two-lobe  propeller  pattern  for  each 
detected  null  of  the  rotating  loop  (i.e,,  two  complete  displays  are 
generated  for  a  single  360°  scan).  The  exploratory  system  produces  only 
one-half  of  a  propeller  pattern  for  each  null,  i.e.,  only  one  display 
is  generated  for  a  single  360°,  This  type  of  display  more  clearly  in¬ 
dicates  the  presence  of  accurate  bearings  through  pattern  symmetry. 


1‘L‘ 

Figuie  1  *)  shows  the  Lost  results  loi  both  Modes  1  and  2.  bssciit  ia  J  ly 
all  ot  the  Mode  1  data  were  obtained  on  s  i  ngLu-hop  transmissions  at  or 
neat  the  FJMIF  with  polarisation  lading  douiintniL  A  i  dii  iade  margin  was 
used  on  a  large  majority  ol  the  Mode  1  Lests  Mode  '2  data  were  obtained 
on  multiple-hop  transmissions  wiLh  a  low  SDR;  in  itiosL  rases,  tlie  signal 
was  peaking-up  out  ot  the  iie>ise  eir  interlerence. 

Figure  id  illustrates  the  better  performance  of  the  exploratory 
system  in  that  the  bearing  spread,  RMS  error,  and  standard  deviation 
values  lor  tne  AN./PRD-5  are  2  to  J  times  greater  titan  those  for  the 
exploratory  system  For  both  systems,  the  bearing  errors  and  spreads 
are  larger  in  the  presence  of  polarization  fading  (Mode  1),  The  field 
tests  revealed  that  it  is  very  difficult  to  obtain  bearings  using  the 
PRD-5  in  the  presence  of  polarization  fading  With  constant  polariza¬ 
tion  rotation,  the  PRD-5  bearing  display  rotates  360c  contv m  o« ■  sly  with¬ 
out  any  noticeable  change  in  pattern  or  swing  rate.  In  order  to  obtain 
data,  a  highly  subjective  operating  technique  using  the  sense  antenna 
was  devised  for  the  PRD-5;  uowever,  this  technique  required  preliminary 
procedures  prior  to  each  test.  The  exploratory  system's  automatic  un¬ 
blanking  capability  allows  for  immediate  bearing  data  even  for  consis¬ 
tent  polarization  fading 

Uperac ional  Evaluation  -  Swing  Rate  Discriminator 

The  SRD  was  designed  primarily  as  an  error  reducing  technique  for 
time-varying  polarization  errors,  hence,  a  large  majority  of  the  field 
tests  were  performed  using  transmissions  with  consistent  polarization 


fading.  The  SRD  was  used  both  alone  and  with  Mode  1 
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i  mil'll-  11  presents,  the  results  lor  .sixty-live  separate  tests.  '1  lit 
AM/PKD-d  data  were  obtained  during  simultaneous,  side-by-side  testing. 
Table  11  shows  slightly  smaller  AO  and  |Ao|  values  for  the  exploratory 
system;  however,  the  must  significant.  i'mproveinent  is  in  the  standard 
deviation  and  RMS  error  values  when  the  SRD  is  used  with  Mode  1.  There 
is  approximately  ;i  J  In  1  difference  Be,.. 'cun  these  parameters  for  the 
exploratory  system,  Ski)  with  Mode  1,  and  the  J’RD-5. 

Bearing  Distribution  Using  SRD  Data 

The  SRD  provides  a  continuous  DC  analog  of  bearing  positions  ver¬ 
sus  time.  Polarization  fading  produces  a  monotonic  ch.  nge  in  the  ana¬ 
log  level  over  one  complete  polarization  relation.  The  DC  analog  data 
can  be  used  to  obtain  histograms  of  bearing  versus  observed  relative 
frequency.  The  number  of  bearings  that  lie  in  specified  azimuth  inter¬ 
vals  is  determined  over  the  entire  range  of  bearing  indications  for  a 
number  of  complete  polarization  rotation  periods.  Figure  20  shows  a 
histogram  obtained  from  the  analog  data  using  20°  azimuth  intervals. 

If  uniform  polarization  is  assumed  (T  =  Kt)  and  the  elevation  angle- 
of--arrival  known,  (3)  may  be  used  to  obtain  a  theoretical  histogram.  An 
example  is  shown  in  Figure  21  for  an  elevation  angle  of  57°  -  the  value 
calculated  for  a  set  of  DC  analog  data  [18] . 

The  calculated  and  observed  histograms  are  very  similar.  Both  show 
maxima  in  the  interval  containing  the  GCB  and  minima  in  the  interval 
which  differs  from  the  GCB  by  90°.  Figure  20  shows  that  the  angular 
rate-of-bearing-swing  is  lowest  near  the  GCB  since  the  largest  number 
of  bearings  occurred  in  the  interval  containing  the  GCB. 


TOTAL  NUMBER  OF  BEARINGS  741) 


DATA  OBTAINED  OVER  7  COMPLETE 
ROTATIONS 


rite  close  similarities  between  the  two  histograms  substantiate  the 


validity  ot  (3)  as  a  model  Lor  describing  bearing  errors  caused  by 
polarization  lading  Also,  the  measured  histogram  data  could  be  used 
t  '  determine  elevation  angle  since  each  elevation  angle  produces  a 
unique  histogram  For  example,  ii  the  elevation  angle  had  been  90', 
the  liistecram  would  have  been  Hat  across  the  .180''  range  whereas  a  0° 
elevation  angle  would  have  a  histogram  with  ail  hearings  in  the  120°  - 
140°  interval , 


Tab!  I!  I  I 


Comparison  of  exploratory  System 
Using  SKD  with  the  AN/PRD-5 


Exploratory  with 

Exploratory  with  SRI)  Mode 


PRl>  5 

SKD  Mode 

and  Mode 

cl’J 

-7’ 

-3° 

-  ] 

11“ 

9 

9 

S  D 

4  3° 

25  = 

13’ 

RMS  Error 

44° 

25  = 

133 

No  •  of 

Bearings 

83  7 

365 

653 

XI  •  OPERATIONAL  EVALUATION  -  TILTED  LOOP 

Field  tests  were  performed  over  a  four-month  period  with  emphasis 
on  sky-wave  signals  originating  from  distances  less  than  600  km  A 
majority  of  the  tests  were  on  one-hop,  high  elevation  angle  (•  30p) 
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signals  dominated  by  polarization  lading  due  to  Faraday  rotation.  The 
tests  eerit'ied  Lhat  the  Lilted  loop  response  is  essentially  as  predicted 
by  theory.  The  experimental  system  was  compared  with,  the  PKU-3  during 
side-by-side  testing.  Results  show  that  the  tilted  loop  technique  can 
provide  about  4  :  1  and  D  :  1  improvement  in  bearing  accuracy  and  dis¬ 
persion,  respectively,  on  transmissions  with  polarization  fading  in 
some  cases,  the  tilted  loop  system  provided  bearing  information  where 
none  could  be  obtained  using  the  conventional  system.  Significant  im¬ 
provements  were  also  observed  on  lower  elevation  angle  transmission 
dominated  by  multipath  effects. 

Figure  22  illustrates  test  results  for  the  exploratory  system  A 
majority  of  the  test  transmissions  were  ICW  signals  wi tli  polarization 
fading  —  a  worst-case  combination  for  loop  DF  systems.  The  instantan¬ 
eous  sensing  significantly  reduced  the  time  required  for  bearing  acqui¬ 
sition  relative  to  the  PRD-5.  The  average  reduction  was  about  2:1; 
however,  in  some  cases,  no  sense  could  be  obtained  with  the  PllD-5  even 
after  tens  of  minutes.  Typical  test  time  was  two  minutes. 

Measures  of  elevation  angle  obtained  from  the  field  tests  showed 
good  correlation  with  the  expected  angles,  which  were  calculated  assum¬ 
ing  E  -  and  F  -  layer  propagation  and  predicted  average  layer  heights. 
The  measured  angles  occurred  within  the  expected  region  on  53  of  60 
tests.  No  ionogram  data  were  available  and  the  expected  angles  could 
not  be  accurately  defined.  However,  the  data  imply  a  +  10°  accuracy, 
Under  conditions  of  exceptionally  good  polarization  fading,  the  ele¬ 
vation  angle  could  be  measured  co  within  several  degrees. 


—  2  +2°  +4  +6  4-8°  +10 

BEARING  ERROR 

Test  Results-Tilted  Loop  and  ICW  Signals. 
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XII.  CON  C I  .JJlCLON  S 

1'echn iques  lo  improve  simple  iuop  11F/DF  system  performance  on  iuno- 
spherical  ly-propugatod  signals  have  been  developed,  implemented  in  an 
exploratory  ilF/DF  system,  and  proven  effective  by  operational  evalua¬ 
tion.  Significant  improvements  in  performance  relative  to  a  contempo¬ 
rary  system  have  been  obLained. 

Bearing  error  reduction  techniques  include  fade  crest  detection, 
swing  rate  discrimination,  and  the  tilted  loop  technique.  The  latter 
two  techniques  provide  for  a  measure  of  elevation  angle-of-arrival . 
System  improvements  include  clutter  reduction  and  synchronous  detection, 
which  increases  sensitivity  and  allows  for  instantaneous  sensing.  Bear¬ 
ing  acquisition  time  has  also  been  increased,  and  operation  cn 
continuous  and  interrupted  carrier  transmissions  is  possible. 

Georgia  Tech  submits  that  many  of  these  techniques  can  be  applied 
to  larger-aperture  DF  systems  to  both  reduce  bearing  errors  and  improve 
over-all  performance. 
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DIGITAL  l' EC  I  IN  IQULS  FOR  RADIO  DIRECTION  FINDING 

1 .  I  NT  KUDUCTI  ON 

The  digital  computer  has  had  and  is  continuing  to  have  a  marked 
impact  on  the  field  of  Radio  Direction  Finding.  Althougn  RDF  was  a 
relatively  mature  Field  oi  study  and  application  by  the  time  the  digital 
computer  made  iLs  appearance,  the  impact  oi  the  digital  computer  on  this 
field  is  a  considerable  one.  The  influence  of  digital  methods  has  been 
most  pronounced  in  three  facets  of  Radio  Direction  Finding: 

1)  Studies  of  RDF  systems  by  means  of  digital  simulation, 

2)  The  acquisition  and  analysis  of  experimental  data  by  digital 
methods , 

3)  Performing  some  of  tire  functions  of  an  RDF  system  by  means  of 
digital  hardware. 

Although  these  three  facets  have  developed  more  or  less  together,  they 
will  be  presented  and  discussed  somewhat  separately. 

II.  SIMULATION  STUDIES 

An  essential  requirement  for  a  radio  direction  finding  system  to 
be  studied  by  means  of  digital  simulation  is  that  it  must  be  possible  to 
describe  the  system  under  study,  or  the  essential  characteristics  of  the 
system  under  study,  by  means  of  a  mathematical  model .  For  purposes  of 
developing  the  needed  mai hematical  model,  it  is  useful  to  divide  the  RDF 


system  into  six  parts: 


JUT 

I)  1'  lit?  signal  which  impinges  upon  the  KDF  sysLcm, 

J)  J  lie  antenna  eienienLs  ot  Llie  allay  (including  any  estimates  that 
can  be  made  of  tile  el  lecl  el  local  site  conditions), 

3)  the  means  by  which  the  signals  1  rom  Llie  antenna  elements  are 
combined  to  give  directional  characteristics, 

■* )  l'he  receiver  portion  ol  the  RDi'  system, 

b)  The  detectors  which  extracL  information  from  Llie  signal  pre¬ 
sented  aL  the  ouLput  of  Lhe  receivers, 

b)  The  display  and  data  presentation  portion  of  Lhe  RDF  system. 

The  considerable  amount  of  work  LhaL  has  been  done  in  RDF  makes  it 
possible  to  write  a  mathematical  description  for  most  parLs  of  the  RDF 
system.  In  addition,  the  linear  nature  of  Lhe  pertinent  characteristics 
of  the  system  (with  the  exception  of  the  detector)  simplifies  consider¬ 
ably  the  model  which  must  be  constructed  and  the  nature  of  the  computa¬ 
tion,  In  the  simulation  process,  the  values  of  the  signal  (or  other 
variables)  are  determined  for  successive  intervals  of  time.  This  allows 
one  to  determine  any  particular  value  at  the  output  or  other  point  in  the 
RDF  system  as  a  function  of  time  under  conditions  of  fixed  system  para¬ 
meters.  It  is  then  possible  to  allow  one  or  more  of  the  parameters  in 
the  system  to  assume  new  values  and  compute  a  new  set  of  values  for  the 
variable  under  study.  Once  a  digital  simulation  model  of  a  system  has 
been  prepared,  it  is  relatively  easy  to  determine  the  behavior  to  be  ex¬ 
pected  from  the  system  under  any  specific  combination  of  values  of  the 
various  portions  of  the  system.  In  one  sense,  the  simulation  model  may 
be  considered  a  means  of  answering  several  different  types  of  questions 
concerning  the  RDF  system.  One  type  of  question  is,  "What  will  happen 
with  this  particular  RDF  system  if  I  have  a  particular  (assumed)  type 
of  input  signal  condition?" 


A  second  type  of  question  is,  "What 
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wiLL  happen  ii  i  ii  ■  RDF  system  itscit  is  made  diliorunL  than  it  is  at 
present  in  some  way?''  So  LhaL  some  idea  may  be  obtained  of  the  useful¬ 
ness  oi  this  means  ol  studying  RDF  systems,  JeL  us  consider  several  typ¬ 
ical  examples  of  simulation  studies. 

Several  studies  have  been  conducted  to  determine  various  character¬ 
istics  ol  t.he  circularly  disposed  antenna  array  (Wullenwcber)  which  is 
used  lor  experimental  work  at  the  University  oi  Illinois..  As  part  of 

(y ) 

the  studies  made  by  Broeker,  it  was  helpful  to  determine  the  response 
oi  tiie  sum  output  oi  the  scanner  to  a  signal  which  consisted  of  two  in¬ 
terfering  wa^ea.  These  studies  showed  'hat  for  interfering  waves  which 
arrived  at  azimuthal  angles  very  close  together,  the  result  was  a  dis¬ 
torted  and  broadened  azimuthal  scan  pattern;  whereas  for  waves  which  ar¬ 
rived  at  angles  separated  by  several  degrees,  the  two  signals  were  ob¬ 
served  as  separate  on  the  azimuthal  scan  pattern.  Typical  output  results 
determined  by  digital  simulation  are  shown  in  Figures  1  and  2.  Another 
means  for  determining  direction  of  arrival  with  the  Wullenweber  system 
compares  the  output  from  the  sum  pattern  and  the  output  from  the  dif¬ 
ference  pattern.  In  the  presence  of  a  single  signal,  the  phase  differ- 

/ 

ence  between  the  right-bank  signal  and  the  left-bank  signal,  determined 
by  using  the  sum  and  the  difference,  can  be  interpreted  in  terms  cf  the 
deviation  of  the  direction  of  arrival  from  the  beresite  of  the  array. 

As  shown  in  Equation (1),  the  direction  of  arrival  and  the  phase,  5,  between 
the  right  bank  and  the  left  bank  is  related  by  a  constant,  k. 


(i) 


(4) 

Figures  3  and  4  are  from  the  studies  made  by  Rosenbaum  and  show 


Figure  1.  SUM  Magnitude  for  Two  Close  Signals. 
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certain  oharacloi  isi  ics  ol  Lhe  sum  and  diflurencu  pattern  as  dcUmuiiiLd 
by  simulation  studies,  wliile  Figure  b  gives  the  value  of  the  constant  k 
in  liquation  (1)  as  a  l  unction  ot  1  requeney  lor  a  fixed  value  of  elevation 
angle  of  arrival.  An  added  question  concerning  this  particular  Lype  of 
system  cent iguiv.t ion  is,  "What  does  this  indicate  concerning  the  direc¬ 
tion  ot  arrival  ot  Lhe  signal  in  the  presence  of  wave  interference?" 
Rosenbaum  used  digital  simulation  to  study  the  question  and  some  of  his 
resulLs  are  shown  in  Figures  b  and  7,  which  give  Lhe  orientation  of  the  el¬ 
lipse  axis  as  viewed  on  the  cathode  ray  oscilloscope  and  the  ellipticity 
of  the  pattern  as  a  function  of  the  azimuth  and  the  relative  magnitude 
of  the  interfering  signal.  It  m  of  interest  to  note  that  strong  inter¬ 
fering  signals  close  (in  azimuth)  to  the  dominant  or  primary  signal  can 
cause  large  azimuthal  errors.  For  example,  a  perturbing  wave,  separated 
by  two  degrees  from  the  primary  wave,  could  cause  an  apparent  shift  in 
the  bearing  by  as  much  as  six  degrees,  i.e.,  to  a  bearing  completely  away 
from  either  of  the  two  signals.  Figures  8  and  9  are  from  the  studies  of 
Medill^^  and  give  a  comparison  of  the  magnitude  of  the  sum  patterns  and 
the  differential  phase  as  obtained  from  the  right  bank  and  left  bank  of 
the  Wullenweber  array.  Additional  studies  of  the  differential  phase,  6  , 

as  a  function  of  frequency,  produced  the  results  shown  in  Figure  10.  In 

(2) 

an  earlier  work,  Crush  showed  that  the  apparent . direction  of  arrival 
as  determined  by  a  Watson-Watt  RDF  system  in  the  presence  of  wave  inter¬ 
ference  could  exhibit  large  variations  in  value,  and  that  these  depended 
upon  the  relative  phase  between  the  interfering  signals.  This  is  shown 


in  Figure  11. 
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Lenweber  array  bad  a  180°  phase  shill  eiieiie  tisly  introduce  I  Figure 
which  shows  part  >»l  Lhe  results  i  ram  taal  sLudy.  indicates  LnaL  at  a 
1  lequei’.cy  ot  9  .  >  Mil/,  the  errors  Lo  be  exported  because  el  Lhis  defect 
are  quite  sma L l  at  the  cophasal  angle  lor  which  the  array  was  designed, 
but  that  they  may  become  as  large  as  2  degrees  lor  tiier  el  vution  angles 
ot  arrival.  Similarly,  one  may  postulate  UiaL  the  cables  connecting  the 
antennas  to  the  equipment  building  are  not  all  the  same  electrical  length, 
but  vary  in  some  pattern  suen  as  a  sinusoidal  variation  over  a  portion 
of  the  cables  with  a  random  variation  superimposed  upon  the  systematic 
variation.  it  would  be  expected  that  this  would  be  the  cause  of  bearing 
error,  at  least  for  some  directions  of  arrival  More  recently,  studies 

have  been  made  of  the  error  in  determining  the  direction  of  arrival  when 
a  large  number  of  adjacent  elements  are  assumed  to  be  missing,^7 ^  The 

plots  of  Figure  13  show  the  error  in  determining  direction  of  arrival  as 
a  function  of  the  difference  in  azimuth  between  the  center  of  the  defec¬ 
tive  Dortion  ef  the  array  and  the  true  direction  of  arrival.  Of  par¬ 
ticular  interest  is  the  fact  that  for  this  large  defect  in  the  array,  the 
error  expected  to  occur  does  not  appear  to  be  large.  Figure  14  shows 
the  sum  patterns  expected  under  these  conditions.  Three  sum  patterns  are 
shown:  a)  the  sura  pattern  far  away  from  the  defect  in  the  array;  b)  the 

sum  pattern  for  the  signal  arriving  from  a  direction  corresponding  to  tiie 
center  of  the  delect;  c)  the  sura  pattern  corresponding  to  the  signal  ar¬ 
riving  from  the  direction  which  yields  the  maximum  error 
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Figure  12 .  Effects  of  180°  Phase  Shift  in  One  Antenna 
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Figure  13.  Error  Pattern  for  the  39  Antenna  Defect  at  4  MHz 


Ill  any  RD  F  sysLern  in  which  di  i  lurui.  L ial  phase  is  tu  be  measured, 
i  L  is  usually  assumed  iliaL  the  receivers  lor  the  iwo  channels  are  matched, 
L . e* ,  ,  tiie  gain  and  phase  shilL  through  the  two  receiver  enannels  are  iden¬ 
tical.  As  identical  receivers  are  difficult  L o  design  and  maintain,  it 
is  important  to  know  the  error  that  will  be  conLribuLed  if  the  two  chan¬ 
nels  are  not  identical.  Figures  15  and  lb  from  Llie  work  of  Henderson  ^  ^  ^ 
indicate  the  error  that  may  be  expected  as  a  consequence  of  phase  and 
gain  mismatch  of  tile  two  receiver  channels  of  a  Watson-Watt  type  RDF 
sys  tem . 

Another  modified  system  of  interest  is  a  system  designed  to  yield 
improved  performance One  of  the  methods  which  has  been  postulated 
for  reducing  the  effect  of  mismatch  error,  as  discussed  in  the  preceding 
paragraph,  is  to  perform  a  channel,  interchange,  i.e.,  exchange  receiver 
channel  A  for  receiver  channel  B  and  take  as  the  correct  value  for  the 
indicated  bearing  the  average  between  the  two  values  determined  under 
the  interchange  and  noninterchange  conditions.  Figure  17  shows  the  re¬ 
sults  of  a  study  which  inquired  what  the  effect  of  such  a  process  would 
be  on  errors  in  determining  the  direction  of  arrival  with  a  Watson-Watt 
system.  It  can  be  seen  that  the  error  will  be  reduced  considerably  but 
will  not  be  removed.  It  can  also  be  seen  that  if  one  knew  the  differ¬ 
ence  in  phase  shift  between  the  two  receiver  channels  and  the  difference 
in  gain,  it  would  be  possible  to  use  a  digital  computer  at  the  output  to 
cori'prt  for  this.  One  way  by  which  this  may  be  done  is  to  connect  a  com¬ 
mon  signal  to  che  two  channel  inputs  at  regular  intervals.  The  relative 
magnitude  of  the  output  and  the  differential  phase  of  the  two  outputs  are, 
then,  a  direct  measure  of  the  imbalance  of  the  two  amplifiers.  This 


Phase,  degrees 


Relative  Phase,  degrees 
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into  rraul  toil  can  Lhcn  be  used  tu  correct  Lhe  values  ol  Llie  direction  u  I 
arrival  which  ate*  determined  wiLh  lhe  nonideal  unbalanced  syslem.  Fig- 
ui'e  IS  shows  lhe  error  that  would  be  expecLed  under  Lhese  conditions  il 
lhe  relative  phase  and  relalive  gain  between  Lhe  Lwo  channels  were  allowed 
lo  drill  as  a  function  ot  Lime,  and  corrections  were  made  aL  intervals 
ot  twenty  successive  readings.  The  drilL  is  posLulaLed  as  being  more 
rapid  at  first  than  later.  IL  is  seen  that  even  for  large  values  oi 
phase  shift  and  gain  imbalance,  Lhe  magnitude  of  the  error  in  determining 
the  direction  of  arrival  is  relatively  small. 

One  of  the  troublesome  aspects  of  the  cathode  ray  tube  display  in 
the  Wullenweber  system  is  that  the  pattern  observed  is  considerably  af¬ 
fected  by  die  modulation  and  fading  characteristics  of  the  signal.  One 
method  postulated  for  reducing  this  difficulty  was  to  utilize  for  the 

/  Q  \ 

displayed  quantity  the  ratio  of  two  signals.  If  the  signals  were 

modulated  by  the  same  value,  then  the  ratio  should  be  constant.  Figures 
19  and  20  are  plots  of  the  sum  and  difference  patterns  which  have  been 
determined  by  means  of  digital  simulation  studies  for  unmodulated  sig¬ 
nals  and  modulated  signals,  respectively.  Figure  21  shows  the  results 
of  dividing  the  sum  pattern  by  the  difference  pattern  and  the  difference 
pattern  by  the  sum  pattern  ,  whereas  Figure  22  indicates  the  sum  and  dif¬ 
ference  patterns  which  would  be  obtained  if  both  were  divided  by  the  volt¬ 
age  on  a  fixed  antenna  which  received  the  signal  but  did  not  contain  any 
directional  information  which  varied  as  a  function  of  the  scan.  A  sys¬ 
tem  which  utilized  analog  signal  processing  to  carry  out  these  functions 
was  constructed  at  a  later  date  and  found  to  perform  essentially  as  pre¬ 


dicted  . 


l'lie  University  ol  Illinois  Wul  Icuwebet  .winy,  as  originally  cun- 


Ligured,  did  not  yield  iuiorm.iL  ion  concerning  elevation  angle  of  arriv¬ 
al.  One  of  the  suggestions  made  lor  adding  this  capability  to  the  sys¬ 
tem  was  to  combine  certain  ol  Llie  antenna  elements  so  as  Lo  form  phantom 
antennas  aligned  in  Lhe  azimuthal  direction  toward  the  signal  source^"*'* 
as  shown  in  Figure  2J  One  of  Lhe  difficulties  that  is  expected  with 
such  a  system  is  that  errors  may  lie  produced  as  a  consequence  of  Lhe 
inability  to  precisely  measure  the  relative  differential  phase  of  the 
phantom  interferometer  pair  One  means  of  simulating  such  an  error  is 
to  simpLy  truncate  the  precision  witli  which  the  differential  piiase  is 
known.  Figure  2a  shows  the  error  to  be  expected  for  various  elevation 
angles  of  arrival  for  two  frequencies  (4  MHz  and  6  Mliz)  if  the  differ¬ 
ential  phase  is  truncated  to  a  precision  of  .1  degree-  It  can  be  seen 
that  for  low  elevation  angles  (5  degrees  or  less)  noticeable  error  may 
be  expected,  even  for  a  relatively  high  precision  differential  piiase 
measurement  such  as  has  been  postulated. 

One  of  the  principal  areas  in  which  digital  methods  have  had  a  pro¬ 
nounced  impact  on  Radio  Direction  Finding  is  in  the  use  of  digital  hard¬ 
ware  to  perform  functions  similar  to  those  performed  by  analog  methods 
either  alone  or  with  human  operators.  The  most  widely  studied  of  these 
is  the  use  of  digital  computation  methods  which  can  be  translated  into 
digital  hardware  for  determining  the  direction  of  arrival.  In  this  type 
of  system  modification,  it  is  postulated  that  the  voltages  available  at 
certain  parts  of  the  system  are  processed  in  some  way  such  as  determining 
the  magnitude  or  the  relative  phase.  These  quantities  are  then  used  as 
input  for  a  digital  computation  of  the  direction  of  arrival  One  of 
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the  principal  advantages  o) lured  by  this  method  oi  operation  is  Li." 
removal  of  Lhe  operator  irom  Llie  direcL  line  e> i  observation  That  is, 
it  is  possible  to  obtain  estimates  of  tile  direction  of  arrival  without 
requiring  the  judgment  and  interpretation  ol  a  human  operator.  The 
method  which  has  been  used  tor  automatic  bearing  computation  at  the 
University  of  llfinois  Wullenweber  system  is  based  on  the  work  of  Smith^^ 
and  is  known  as  the  "Center  of  Gravity"  method.  Thu  centroid  ("Center 
of  Gravity")  of  the  major  lobe  of  Lhe  array  response  as  a  function  of 
azimuthal  position  of  the  scanner  is  determined.  The  displacement  of  the 
centroid,  C,  from  an  arbitrary  reference  is  given  by  liquation  (3)  where 
x  is  the  azimuthal  displacement  from  the  given  reference  where  the 

n 

Z  xi>F  (xi) 

c  =  A:i _  (3) 

n 

I  F  <*.) 
i  =  1 

values  of  the  sum  pattern,  F(x^),  are  known.  The  principal  problem  in 
implementing  this  bearing  computation  procedure  is  to  devise  a  means  by 
which  those  points  which  are  to  be  considered  as  the  major  lobe  of  the 
response  can  be  determined.  In  the  method  currently  used,  the  maximum 
signal  level  is  first  determined  and  the  first  point  to  the  right  and  to 
the  left  of  the  maximum  which  falls  below  10  per  cent  of  the  maximum  is 
noted.  (See  Figure  25.)  Those  points  between  the  10  pet  cent  level  de¬ 
fine  the  major  lobe.  If  no  point  falls  below  the  10  per  cent  point  before  a 
minimum  is  reached,  then  the  first  minimum  on  either  side  of  the  maximum 
is  used  to  define  the  limit  for  the  major  lobe.  One  of  the  errors 


associated  with  this  me L hud  of  determining  Ur  direction  of  arrival  is 
that  the  sum  pattern  is  sampled  at  regular  intervals.  Thus  an  error  may 
be  found  which  depends  upon  Lhe  relative  positions  of  the  direction  ol 
arrival,  and  the  sampling  azimuLhs.  Figure  2b  is  a  plot  of  Lhe  errors 
determined  through  digital  simulation  sLudies  of  Lhe  center  of  gravity 
method  ot  bearing  computation.^  Not  shown  here  is  Lhe  considerable 
error  which  comes  about  when  Lhe  number  of  samples  per  revolution  is 
small;  or,  what  amounts  to  the  same  tiling,  Lhe  number  of  samples  that  is 
included  within  the  major  lobe  is  small,  digital  simulation  studies 
have  shown  that  when  the  number  of  samples  included  in  the  major  lobe 
becomes  less  than  twelve,  the  error  increases  noticeably. 

As  was  noted  ’  earlier ,  the  difference  between  the  azimuthal  direction 
of  signal  arrival,  a,  and  the  boresite  position  of  the  scanner  is  re¬ 
lated  to  the  differential  phase  between  the  right- bank  voltage  and  the 
left-bank  voltage  by  Equation  (4  ).  This  is  shown  in  Figure  27, 


^YR  yl) 


For  a  fixed  azimuthal  direction  of  signal  arrival,  Equation  (5)  -.hows  the 
relation  between  the  differential  phase  and  the  scanner  boresite  azi- 


(yR  -  Yl)  =  k  x  +b  (5) 

muthal  position,  x.  As  b  is  the  differential  phase  for  a  zero  value  of 

azimuth,  the  scanner  boresite  position  v;hen  the  differential  phase, 

( f  -  y  )  ,  is  equal  to  0  is  the  direction  of  arrival  of  the  signal  as 
R  L 


given  by  Equation  (6) ■ 
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Figure  26.  Errors  in  Automatic  Hearing  Computer 
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I'lie  differential  phase  is  measured  ai  regular  intervals  of  scanner  bore- 
site  position  and  a  least  squares  fit  Lo  a  sLraighL  line  made  tc  Lhese 
data.  i'lie  slope  and  intercept  lor  the  straight  iine  will  correspond  to 
the  value  ot  k  and  tile  value  of  b,  respectively.  From  this,  the  direc¬ 
tion  ot  arrival  of  a  signal  may  be  readily  determined. ^ 

in  an  earlier  work,  Grusl/“)  allowed  that  a  digital  computer  could 
be  used  to  determine  the  direction  of  arrival  in  a  Watson-Watt  KDF  system 
by  measuring  the  amplitude  and  relative  phase  of  the  two  signals  which 
are  used  to  form  the  display.  The  algorithm  for  this  computation  is  that 
of  equation  (7)  in  which  e^  and  are  the  magnitudes  of  the  two  voltages 
and  u  is  the  differential  phase  between  them. 


B  =  j  Tan  1 


ie.  e„  cos  0 


el  '  e2 


An  alternate  form  is  shown  in  liquation  (8)  in  which  e^  is  the  magnitude  of 
the  phasor  sum  of  e^  and  e Thus,  if  measurements  of  the  magnitude  of 
e^ ,  and  e^  are  available  and  used  as  input  for  the  computational  pro¬ 
cess,  the  direction  of  arrival  may  be  readily  determined. 


„  1^-1 
B  =  —  Tan 


2,2  2 
es  -  (al  +  e2  ) 
2  2 
el  "  e2 


It  should  be  clear  that  if  one  assumes  the  availability  of  a  digit¬ 
al  computer  and  the  ability  to  determine  by  direct  measurement  the  values 
of  appropriate  quantities  in  the  RDF  system,  and  uses  these  values  as  in¬ 
put  for  the  digital  computer,  the  digital  computer  can  be  made  to  func¬ 
tion  as  a  digital  bearing  computer  and  yield  estimates  of  the  direction 
of  arrival  on  an  automatic  basis. 
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In  many  arcus  ol  research,  iL  lias  been  lound  Lo  be  desirable  to 
pre'eess  the  experimental  daLa  wiLh  a  digiLai  computer  One  of  the  prin¬ 
cipal  reasons  Lor  doing  this  is  Lin*  speed  and  accuracy  with  which  a  large 
volume  of  experimental  data  can  be  processed  .  Studies  in  RDF  and  in  other 
areas  in  which  direction  finding  p.ays  an  important  part  are  no  exception 
The  large  volume  of  daLa  makes  it  extremely  impractical  to  consider  .ac¬ 
quiring  data  by  any  other  than  automatic  methods.  Further,  if  Lne  data  are 
to  be  introduced  to  a  digital  computer  for  computation,  it  is  quite 
desirable  tnat  the  data  be  coded  into  a  digital  format  suitable  for  direct 
input  to  a  digital  computer  at  the  earliest  possible  stage. 

Although  a  wide  variety  of  data  can  be  recorded  for  input  to  a 
computer,  the  experimental  studies  conducted  by  the  Radiolocation  Re¬ 
search  Laboratory  have  tended  to  record  those  quantities  which  can  be 
directly  measured  with  somewhat  conventional  instrumentation.  Thus, 
rather  than  perform  a  computation  to  determine  the  direction  of  arrival 
on  an  on-line  basis  and  record  the  direction  of  arrival  thus  obtained, 
the  voltages  and  differential  phases  which  could  be  used  to  determine 
direction  of  arrival  have  been  recorded  and  used  as  input  to  the  computer 
from  which  the  direction  of  arrival  computation  has  been  subsequently 
made.  The  subject  of  data  acquisition  for  digital  computer  processing 
is  a  topic  of  its  own.  It  is  important  to  recognize  that  when  using  an 
automatic  data  acquisition  system,  one  is  bound  by  previous  decisions 
that  have  been  made  concerning  the  data  that  are  to  be  acquired.  For 
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Lhis  reason,  cons  idcL able  thought  should  be  given  to  tlie  design  ol  the 
data  acquisition  system  and  Lo  Liie  manner  in  which  iL  is  to  be  used  be¬ 
fore  (.lata  taking  stalls  Some  ol'  Lhe  items  Lo  be  considered  include: 

a)  What  data  should  be  recorded, 

b)  What  is  lhe  precision  needed  for  the  various  data, 

c)  What  is  the  format  to  be  used  for  the  recorded  data, 

d)  Ik’w  rapidly  should  ttie  data  be  recorded. 

It  the  two  data  acquisition  systems  used  by  the  Radiolocation  Research 
Laboratory  are  considered,  the  answers,  for  these  two  cases,  may  be  seen. 
The  data  acquisition  system  (DAS)  used  to  acquire  data  from  the 

/  I  U\ 

interferometer  RDF  systeni  is  shown  in  the  block  diagram  of  Figure  28. 
The  data  and  the  precision  needed  which  were  determined  to  be  desirable 
for  this  system  are  as  follows: 

Three  values  of  differential  phase,  12  bits  (0.1°)  for  each  value; 

Three  values  of  signal  amplitude,  12  bits  for  each  value; 

Che  value  of  calibration  phase  and  one  value  of  a  calibration 
amplitude,  12  bits  for  each  value; 

Time  (BCD)  hours,  minutes,  seconds; 

Frequency  (BCD)  Megahertz  (to  39),  Kilohertz  * 

The  format  for  recording  these  data  utilizes  8  channel  punched  paper 
tape.  Six  channels  are  used  for  data  bits,  one  channel  for  a  parity  bit, 
one  channel  for  control  and  possible  expansion.  Thus,  each  phase  or  am¬ 
plitude  measurement  requires  two  characters  and  by'  squeezing  the  time  and 
frequency  data  together,  these  data  may  be  contained  in  22  bits  each, 
thus  utilizing  all  but  2  bits  of  4  characters.  By  proper  encoding,  these 
2  bits  can  indicate  whether  the  4  character  group  is  time  datum,  frequency 
data,  or  some  other  kind  of  data  as  indicated  in  Table  1.  Thus, 


Figure  28.  System  block  diagram. 
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Table  1 

Encoding  Pattern  lor  Time  and  Frequency  Data 
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Llif  lutdi  record  oi  data  loiis  Lsl  ing  ul  Lime  or  1  tc-queiicy ,  plus  Liu-  ik-li-s- 
sary  phase  and  ampliLude  measurements,  consists  ui  30  data  characters 
plus  an  end  at  record  ciiaracler  —  a  LoLal  ol  21  ctiaracters  as  shown  in 
Table  2.  L'he  speed  wiLli  whicn  Lire  data  should  be  taken  Is  determined  by 
the  rate  of  variation  of  Lire  daLa  sampled  However,  the  maximum  speed  or 
tiie  maximum  rate  is  determined  by  hardware  limitations  For  tills  system, 
the  limitation  is  the  speed  oi  the  paper  Lape  punch  whicn  lias  a  maximum 
of  2aU  punched  characters  per  second  This  implies  a  minimum  Lime  be¬ 
tween  data  records  of  a  little  bit  less  than  90  milliseconds 

The  other  system  tor  acquiring  data  for  direct  input  to  the  computer 

is  the  high  speed  data  system  (USDS),  shown  in  Figure  '29,  located  at  the  Bond- 

(19) 

viile  Road  Field  Station.  As  this  system  is  designed  to  acquire  data 

from  the  Wullenweber  RDF  system,  it  is  designed  to  take  samples  of  the  sum 
and  difference  amplitudes  as  well  as  the  differential  phase  between  the 
rignt  and  left  bank  voltages  at  the  output  of  the  scanner  for  several  po¬ 
sitions  during  rotation  of  the  scanner.  The  number  of  positions  at  which 
these  three  quantities  may  be  sampled  per  scan  varies  and  may  be  as  large 
as  90,  although  30  points  per  scan  (one  sample  per  degree)  have  been  used 
for  most  of  the  experimental  work.  In  addition,  values  of  eleven  other 
quantities  are  sampled  once  per  scan,  making  a  total  of  101  data  samples 
for  eacn  rotation  of  the  scanner,  Eacn  data  sample  has  two  parts-  one 
part  is  a  code  to  identify  which  of  the  pieces  of  data  the  second  part 
(whicn  is  the  value  of  the  data)  represents.  As  there  are  16  possible 
channels  of  input  data,  the  4  bit  code  is  used  for  encoding  the  channel 
number.  A  14  bit  binary  code  is  used  for  encoding  the  value  of  the  quan¬ 
tity  recorded.  Although  the  information  is  recorded  on  a  nine  track 
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Format  ul  Data  Record 


21-Character  Block 
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No. 


31-CharacLer  Block 


Frequency  or 
1  Line  Data 


Phase  (2,1) 


Amplitude  (2) 


Phase  (1,3) 


Amplitude  (1) 


Phase  (3,2) 


Amplitude  (3) 


Phase  (common) 


Amplitude  (common) 


Longitudinal  Parity  and 
End  of  Block  Marker 


Characters  1  to  20 
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the  21  Character  Block 


Phase  (2,4) 


Amplitude  (4) 


Phase  (4,5) 


Amplitude  (5) 


Phase  (4,2) 


Longitudinal  Parity  and 
End  of  Block  Marker 


Figure  29.  Simplified  Block  Diagram  for 


Thus  , 
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digital  magnetic  tape,  only  six  tracks  arc*  used  lor  recording  daLa. 
three  diameters  are  needed  Lo  encode  a  single  sample  of  data.  This  system 
is  capable  ol  recording  at  a  much  more  rapid  raLe.  A  complete  data  sample 
may  be  recorded  in  between  100  and  20U  microseconds.  This  allows  Lhe  101 
data  points  per  scan  to  be  recorded  in  abouL  10  Lo  20  milliseconds. 

Although  the  recorded  data  are  in  a  format  suitable  for  direct  input 
Lo  the  computer,  they  are  not  used  for  .irect  computation.  Rather,  the 
data  are  edited,  errors  are  eliminated,  and  the  error-free  data  are  stored 
in  a  revised  format  in  a  data  file  which  may  be  more  readily  accessed  for 
subsequent  calculations  than  can  the  original  data.  All  the  data  in  a  com¬ 
plete  record  (that  needed  for  one  calculation  of  direction  of  arrival  of 
the  interferometer  or  that  for  one  scan  with  the  Wullenweber  system)  are 
maintained  in  the  file  as  a  single  data  record.  The  data  record  has  as 
identifying  keys  the  date  and  time  at  which  a  record  was  made.  If  more 
than  one  data  record  has  the  same  date  and  time,  a  sequence  number  is  given 
to  indicate  the  different  times  to  be  assigned  to  each  of  the  records. 
Figure  30  shows  three  plots  produced  by  a  digital  computer  from  the  data 
recorded  at  the  BRFS.^^ 
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DIGITAL  DF  APPLICATIONS 

A  growing  segment  of  the  DF  community  has  begun  to  apply  digital 
techniques  to  radiolocation  problems.  Computer  processing  and  digital 
logic  are  applicable  to  receiver  tuning,  bearing  calculations,  remote 
direction  finder  operation,  and  other  tasks  requiring  sophisticated 
data  processing.  In  this  report,  we  will  describe  digital  DF  exper¬ 
iments  and  hardware  completed  in  our  laboratory  for  special  radio¬ 
location  problems. 

Multichannel  receivers  are  necessary  for  certain  radiolocation 
applications.  For  fast  response  time  or  careful  amplitude  and  phase 
measurements,  at  least  two  phase  and  gain  matched  receiver  channels 
are  used.  Ultimately,  as  the  demands  for  short  response  time  and 
sophisticated  processing  increase,  receivers  with  three  channels  or 
more  will  be  used.  The  history  of  the  development  of  twin  channel 
receivers  shows  significant  technical  difficulties  in  achieving  gain 
and  phase  match  of  the  order  of  1  dB  and  10°  over  the  HF  band. 

The  conventional  amplitude  comparison  twin  channel  direction 
finder  or  Watson-Watt  system  uses  a  CRT  with  the  bearing  given  by 
the  display  ellipse  major  axis.  Figure  1  shows  the  ellipse  orientation 
equation.  When  the  amplitude  and  phase  difference  between  the  two 
channels  are  measured,  the  major  axis  is  readily  computed,  for  example, 
using  a  digital  computer  without  the  need  for  CRT  display.  More 
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important  tluu  Lius,  howevot,  i  a  tlic  capability  lot  i  ompinsa  L  i  up  1  ,r 
phase  aiul  gain  mismatch  bet  ween  Lite  channels  using  simple  receiver 
calibration  data.  Figure  2  shows  a  simple  lot  muiat  ion  wtiere  Uie 
amplitude  and  pliase  mismaLch  is  determined  by  the  injection  oi  equal 
in-phase  signals.  The  mismatch  raLios  (K)  are  applied  as  cumpeusul i ng 
iactors  in  Lite  bearing  equation.  This  simple  approach  oilers  Ltie 
capability  for  tolerating  receiver  mismatch  conditions  when  Liiey  arc- 
known  .  Figure  3  shows  a  block  diagram  of  a  simple  verification  of  Lhis 
compensation  technique  using  a  laboratory  DVM  and  pliase  meter. 

Figure  A  shows  a  comparison  of  calibration  curves  obtained  using 
crossed  spaced  loop  antennas  for  vertical  polarization.  The  bearing 
is  obtained  by  amplitude  comparison  of  the  sine-cosine  related  antenna 
response.  The  calibration  curve  at  the  bottom  is  the  conventional 
analog  calibration  curve  obtained  reading  the  CRT  ellipse.  The  calibra¬ 
tion  curve  at  the  top  was  obtained  using  the  major  axis  computation 
measured  from  the  digitized  amplitudes  and  phase  difference.  The 
difference  between  these  curves  is  essentially  determined  by  the 
residual  precision  of  the  two  measurements,  and  amounts  fco  approximately 
+1  1/2°,  Figure  5  shows  an  analog  calibration  curve  at  the  bottom 
taken  with  3  dB  receiver  gain  mismatch  and  30°  phase  mismatch.  These 
mismatch  values  were  deliberately  set  into  the  receiver  and  the 
resulting  calibration  curve  shows  i8°  octantal  error  as  expected. 

The  top  calibration  curve,  however,  is  obtained  after  numerical 
compensation  of  the  receiver  channels.  The  bearing  performance  is 
essentially  equivalent  to  that  obtained  previously  with  matched  channels. 
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1'lic  must  important  implication  ol  this  technique  is  the  clear 
possibility  LhaL  multiple  channel  receivers  wiLh  nominal  phase  and 
gain  match  can  bo  compensated  by  digital,  calibration  data  updated 
as  required  in  the  1)1'  compuLaLion.  This  digiLai  compeusac  ion  technique 
should  be  as  tolerant  as  b-lU  dll  mismatch,  and  may  s  ign i  1  icant ly 
increase  the  feasibility  oi  multichannel  receivers  A  multiplicity 
oi  applications  would  immediately  arise  when  they  become  available. 

Figure  b  shows  the  phase  reiaLion  between  vertically  polarized 
sine-cosine  antennas  for  sense  finding.  For  a  system  such  as  the 
spaced  loop,  the  phase  reverses  in  each  octant  correspond  to  the  20 
amplitude  response.  At  the  bottom  of  Figure  6,  the  quadrant  and  phase 
data  are  summarized  into  a  truth  table  where  the  (0)  implies  an  in- 
phase  condition  and  the  (1)  implies  180°  phase.  The  sense  function 
around  the  circle  thus  is  given  by  binary  words,  each  word  corresponding 
to  a  different  octant  for  the  spaced  loop  or  quadrant  for  the  simple 
loop.  Figure  7  shows  a  block  diagram  of  a  four-channel  receiver  system 
for  determining  octant  sense  by  instantaneous  phase  comparison  of  the 
spaced  loop  and  simple  loop,  DF  and  sense  antennas.  It  should  be 
emphasized  that  in  sense  finding,  it  is  not  necessary  to  make  a  precise 
p.tase  measurement  but  rather  to  determine  whether  the  DF  and  sense  signal 
are  more  nearly  in-phase  than  out-of-phase.  The  nominal  tolerance  is 
t60°.  Figure  8  shows  a  conceptual  implementation  of  this  technique 
where  the  sense  logic  lights  the  correct  octant  of  the  alidade, 

Figure  9  shows  a  digital  sense  block  diagram  of  a  circuit  built  and 
tested  for  a  three-channel  receiver  (since  the  four  channels  were  not 
available)  where  fast  acting  logic  gates  sense  positive  going  or 
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UL'gJLiVL'  going  slopes  ol  the  IF  wavuionn  Lo  deLerm  Lne  relative  in-phase 
or  out-ol-phase  conditions  across  L lie.  array.  A  pair  ol  DF  and  sense 
antennas  are  sampled  simultaneously  Lhus  requiring  Lwu  steps  to  acquire 
a  complete  set  of  pliase  data.  The  sampling  sequence  is  controlled  by 
the  monopole  reference.  One  complete  set  ol'  phase  data  can  be  obtained 
in  less  than  10U  microseconds.  This  digital  approach  to  sense  finding 
eliminates  the  CRT  sense  display  and  permits  sense  finding  simultaneous 
with  the  DF  calculation. 

Another  application  of  digital  control  circuits  is  programmed 
azimuth  scan  of  a  DF  array.  Figure  10  shows  an  AN/F11D-10  System, 
bach  of  three  racks  in  the  center  are  180  Stage  Commutators  which 
perform  DF  azimuth  scan.  Figure  11  shows  typical  DF  displays  obtained 
with  the  AN/FRD-10  DF  system  connected  to  a  circular  array  of  Beverage 
antennas.  Sum  and  difference  polar  patterns  are  shown  as  we.  1. 
the  inverted  difference.  The  thing  to  emphasize  is  that  the  azimuth 
scan  is  performed  by  a  180-stagn  solid  state  ring  commutator  which 
permits  simultaneous  digital  selection  of  an  arbitrary  number  of 
antennas  at  arbitrary  locations  in  the  array.  The  selected  array  is 
then  sequenced  in  azimuth  around  the  azimuth  circle.  Cogging  and 
other  pattern  deteriorating  characteristics  are  absent  in  this  design 
The  scan  rate  is  10  cycles  per  second  chosen  equal  to  the  conventional 
Wullenweber  scan  rate  but  is  limited  only  by  the  receiver  bandwidth. 
Figure  12  shows  a  photograph  of  the  digital  commutator  and  the  azimuth 
readout  strobe.  The  circular  sweep  is  generated  electronically  with¬ 
out  a  mechanically  rotating  clock.  Either  the  strobe  or  the  DF  pattern 
can  be  rotated  to  read  out  the  bearing.  Because  the  Beverage  antenna 
elements  are  directive,  simple  summing  was  provided  in  the  RF  switch. 
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switch.  Each  ol  tlie  I8U  antennas  is  1  ed  in  a  balanced  transistor 
gate.  Each  gate  provides  8  dli  ampi  il  teat  ion  matched  to  il  d  li  witli 
phase  m.iLcl.  better  than  'b3  The  gale  is  turned  on  and  oil  by 
application  ot  a  square  wave  Lu  Lhe  switeliing  signal  input*  The 
balanced  design  eliminates  DC  oi'iset  switching  transients  at  the  gate 
ouLpuL.  A  10  dll  (S+N)/d  ratio  is  obtained  tor  an  input  of  less  than 
1^  volt  RMS. 

The  use  of  digital  switching  permits  very  generalized  computer 
control  of  the  azimuth  scan  process  for  circular  or  other  arrays. 

These  techniques  can  be  applied  to  the  spatial  filter  type  (Wullenweber ) 
array  or  interference  type  direction  finding  arrays  such  as  the 
interferometer.,  Figure  14  shows  the  block  diagram  of  the  Army  inter¬ 
ferometer  system  in  operation  in  our  laboratory  emphasizing  the  use 
of  digital  circuits  and  an  on-line,  real-time  digital  computer  for 
bearing  calculation.  The  multiple  phase  cycle  ambiguities  in  the 
long  baselines  are  resolved  by  an  appropriate  algorithm  using  phase 
measurements  from  the  short  baselines.  The  8K  DDF  516  computer 
provides  sufficient  core  storage  for  basic  on-line  direction  finding 
and  single  site  location  programs  with  operator  interactive  features, 
fn  addition,  the  98K  disc  file  to  be  installed  shortly  with  the 
computer  will  remit  rapid  access  input  and  output  of  other  computer 
programs  and  storage  of  up  to  75K  words  of  data  for  subsequent  off¬ 
line  processing. 

Data  editing  is  performed  by  digital  logic  in  two  ways.  The 
threshold  circuit  inhibits  the  A/D  converter  except  when  a  set 
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threshold  is  I’xo’cilud. 


1  n  Lliis  way,  interrupted  carrier,  Sbii,  and 
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other  such  modulations  i  an  be  liandled  to  eliminate  A/ii  conversion 
when  no  signal  is  present.  The  plane  wave  test  circuit  provides  an 
additional  threshold  criterion  by  requiring  equal  amplitudes  within 
a  set  tolerance  before  the  A/D  converter  samples  the  phase  meter 
At  present  the  tolerance  is  no  more  restrictive  than  ldii. 

Figure  15  shows  the  comparison  of  snap  bearings  obtained  with 
and  without  data  editing.  Where  signal  fading  and  interrupted  carrier 
occur,  this  data  editing  feature  is  indispensable.  As  currently 
designed,  the  threshold  circuit  permits  the  A/D  converter  to  take  phase 
samples  (at  the  rate  of  1  per  120  microseconds)  for  a  32  ms  interval  on 
each  sequencer  state.  These  phase  readings  are  averaged  for  the  phase 
angle.  If  the  threshold  continues  to  be  satisfied  after  32  ms,  the 
sequencer  switches  to  the  next  state  and  so  on  such  that  eight  frames 
of  DF  data  can  be  obtained  per  second.  A  specified  number  of  azimuth 
and  elevation  frames  rs  then  averaged  for  a  single  answer  typed  out  to 
the  operator. 

A  similar  five-element  interferometer  is  being  assembled  for  the 
Navy  at  the  simulated  ship  site  at  our  laboratory-  The  simulated  ship 
at  the  center  of  the  site  is  a  full-scale  electrical  mock-up  of  a  typi¬ 
cal  destroyer  superstructure.  Figure  16  shows  a  plan  view  of  the 
simulated  ship  site  and  the  interferometer  arrangement  for  the  first 
siting  configuration. 
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READY: 

READY  :SY 
SEQUENCE  B: 

FREQUENCY  (MHZ)  8.  500:6 
RUN  NO.  70  22MY69: 


With  Plane  Wave  Editing 
Elevation  Spread  -  1.0 

Azimuth  Spread  -  9.  3 


Without  Plane  Wave  Editing 
Elevation  Spread  -  1.  5° 

Azimuth  Spread  -  18.2 


With  Plane  Wave  Editing 
Elevation  Spread  -  0.7° 
Azimuth  Spread  -  10.  5° 


FIGURE  15 

BANDERA  TARGET,  FREQUENCY  8.  5  MHz,  RANGE  55  km 
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The  long  baseline  spaciugs  are  1(J‘j  leeL  ami  2/0  leeL.,  respec  Lively 
The  receiver  and  phase  meLer  are  located  in  l lie  Laboratory  area  at 
the  site.  The  interferometer  computer  is  located  approximately  3/4  mi  le 
northwest  ol’  the  ship  site..  The  remote  DDP  bib  computer  Is  used  to 
perform  on-line  remote  processing  of  the  ship  site  Interferometer  data. 
The  standard  telopnone  dial  system  was  chosen  as  the  most  convenient 
medium  Lor  transmitting  information  Lo  and  from  the  computer  and  remote 
interferometer.  in  addition,  these  techniques  provide  a  versatile 
approach  for  remote  computer  control  and  data  processing  at  other 
sites . 

Figure  17  shows  a  block  diagram  of  the  remote  interferometer 
system.  Four  phase  angles  are  measured  using  the  antenna  located  at 
tne  apex  of  the  L  for  reference.  The  equipment  used  here  employs 
standard  laboratory  instruments  to  minimize  set-up  time  for  the  system. 
The  IF  outputs  of  the  twin  receiver  are  fed  to  an  analog  phase  meter. 

A  commercial  digital  voltmeter  performs  the  necessary  A/D  phase  angle 
conversion  with  BCD  output  for  three  decimal  digits.  The  sequencer 
assembly  provides  the  sample  and  hold  instructions  to  the  DVM.  After 
encoding  the  phase  angle  and  state  number,  a  standard  shift-register 
converts  the  parallel  input  to  serial  output  and  feeds  the  acoustic 
telephone  coupler.  The  Data  Set  at  the  remote  computer  converts  the 
received  tones  to  serial  digital  data  for  input  to  the  516  computer. 

Interferometer  control  is  provided  from  a  teletype  unit  at  the 
ship  site.  The  TTY  terminal  operates  at  10  characters  per  second; 
therefore,  it  requires  approximately  1.7  seconds  to  complete  a  four- 
phase  angle  frame  measurement,.  The  time  averaged  azimuth  and  elevation 
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answers  are  ted  from  the  computer  to  the  control  TTY. 
system  is  ties  Lgned  ior  rudimentary  skvwave  perlormunce  in  tlie  sliip 
environment.  Tile  sophistication  available  in  the  Army  system  can 
be  incorporated  as  desired;  however,  ior  the  high  speed  data  transfers 
required,  a  broadband  d4(JU  bil/see  phone  circuit  is  needed. 

We  have  summarized  a  number  of  applications  to  radiolocation 
systems  for  digital  computation  and  logic.  Ultimately,  time  sharing 
software  available  witli  modern  system  computers  makes  possible  tne 
simultaneous  operation  of  a  number  of  remote  DF  systems  and  experiments. 
L)F  systems  can  be  remoted  not  only  across  the  campus  but  to  any 
compatible  computer  using  radio  digital  communication  links  or  even 
the  direct  dial  long  distance  telephone.  Applications  requiring  this 
technology  are  sure  to  increase  in  the  coming  years. 
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L.  J.  MILLER 

J’Lpt wtint'.nl  of  Eloo',  vioai  Eiiijinoorinj 
L'n:  Ly  oj  lll:noit}y  Uvbuna ,  Illinois 

ELEVATION  AN OLE  MEASUREMENTS  UN  THE  UNIVERSITY  UE  ILLINOIS 
WULLENWEBER  ANTENNA  ARRAY 

I .  INTRODUCTION 

la  an  attempt  Lo  measure  elevation  angles  using  the  Wullenweber 

Antenna  Array  at  the  University  of  Illinois,  two  methods  have  been  used. 

The  first  method  was  used  by  Westlund  to  obtain  elevation  angle  data 

(2) 

while  the  second  method  is  based  upon  a  system  proposed  by  Bailey. 

II.  WESTLUND  METHOD 


The  Westlund  method  of  determining  the  elevation  angle,  as  applied  to  the 

Wullenweber  Antenna  Array,  is  shown  in  Figure  1.  This  method  is  an  appli- 

(3) 

cation  to  the  Wullenweber  system  of  the  Wilkins-Minnis  method  of  ele¬ 
vation  angle  determination.  It  consists  of  two  vertically  stacked,  ver¬ 
tically  polarized  antennas  positioned  at  heights  z^  and  above  the 
earth.  The  lower  antenna,  at  height  z^,  is  one  °f  the  folded  monopole 
antennas  of  the  Wullenweber  Antenna  Array.  The  second  antenna,  at  height 
z^,  is  a  dipole  antenna  which  is  supported  off  one  of  the  poles  used  to 
support  the  reflecting  screen.  This  dipole  antenna  is  mounted  directly 
over  the  folded  monopole  antenna. 

For  this  configuration,  it  can  be  shown  that  for  the  ideal  case  the 
corresponding  antenna  terminal  voltages  are: 
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incident  magneLie  field  intensity 
intrinsic  impedance  ul  free  space 
elevation  angle  measured  from  the  ground  plane 
wavelength  aL  Lite  operating  frequency 

vertical  distance  above  the  ground  plane  of  Lite  dipole  and 
monopole,  respectively 

horizontal  distance  from  the  vertical-wire  reflecting  screen. 


From  these  equations,  it  is  seen  that  the  antenna  terminal  voltages 
have  the  same  time  phase.  It  is  further  seen  that  the  amplitude  of  the 
voltage  at  each  of  these  antenna  terminals  is  a  function  of  the  elevation 
angle  of  arrival,  the  frequency  of  operation,  and  the  respective  height 
of  these  antennas  above  ground.  Since  for  any  given  elevation  angle 
measurement,  the  frequency  of  operation  and  the  height  of  each  of  these 
antennas  above  ground  are  known,  then  the  amplitude  of  each  of  these  an¬ 
tenna  terminal  voltages  is  dependent  upon  the  elevation  angle  of  arrival 
6.  Therefore,  to  determine  the  elevation  angle  of  arrival  by  this  method, 

i 

some  means  of  extracting  this  information  from  the  two  signal  amplitudes 
is  required. 

One  method  of  extracting  this  information  is  to  apply  the  output  of 
each  antenna  terminal  to  one  channel  of  a  matched  twin-c.hannel  receiver. 
At  the  IF  output  of  the  receiver,  prior  to  detection,  these  signals  are 
then  applied  to  a  conventional  X-Y  display  oscilloscope.  This  will  give 
a  line  or  ellipse  on  the  display  at  an  angle  y  such  that: 
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by  measuring  the  angle  y  from  the  display  scope,  llio  elevation  angle 
0  may  be  determined.  This  is  the  method  used  by  WesLiund  to  obtain 


elevation  angle  data. 


Displaying  and  determining  elevation  angle  information  in  Lliis 
manner  is  acceptable  if  rather  infrequent  or  manual- Lype  measurements 
are  to  be  employed  to  record  this  information.  However,  if  an  auto¬ 
mated  data  recording  system  is  to  be  employed,  sucli  as  the  High-Speed 
(4) 

Data  System  (USDS)  which  is  normally  used  to  record  data  obL^-ned  by 
use  of  the  University  of  Illinois  Wullenweber  Antenna  Array,  then  some 
other  method  of  extracting  this  information  from  the  signal  amplitudes 
is  required.  To  record  these  data  using  the  HSDS,  one  approach  is  to 
detect  the  antenna  voltages  E  ^  and  E  at  the  output  of  the  twin-channel 
receiver.  These  voltage  amplitudes  are  then  recorded  by  the  HSDS  and 
the  elevation  angle  3  is  determined  by  programming  the  G-20  computer  to 
solve  the  equation; 


2ttz^ 

cos ( — : -  sin  3) 
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2itz2 

cos  ( — - —  sin  0) 


(4) 


At  the  present  time,  v  crk  is  in  progress  to  reinstrument  this 
method  of  elevation  angle  measurement.  The  data  obtained  will  be 
recorded  by  the  HSDS  and  processed  by  the  G-20  to  determine  the 
elevation  angle  6, 

This  method  of  determining  elevation  angle  when  used  in  conjunction 
with  the  U  of  I  Wullenweber  Array,  or  any  other  Wullenweber  type 
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ot  antenna  array,  has  two  disadvantages.  The  llrst  disadvantage  is 
that  in  using  this  method  iL  is  noL  possible  to  obtain  Lhe  required 
intoriuation  by  use  oi  just  those  antennas  which  already  make  up  the 
antenna  array.  Additional  dipoles  or  other  types  of  antenna  elements 
are  required.  The  second  disadvantage  is  Lhat  for  measuring  elevation 
angles  at  different  azimuths,  these  additional  antenna  elements  must 
be  located  on  the  array  at  those  azimuths  of  interest. 

ill.  BAILEY  MbTiiOL) 

The  Bailey  method  of  determining  elevation  angle  is  directly 
applicable  to  the  Wullenweber  Antenna  Array  and  is  based  upon  the 
technique  of  measuring  the  differential  phase  between  two  horizontally 
displaced  antenna  elements.  Theoretical  studies  of  this  method  by 
Jones,  Schlicht,  and  Ernst^  introduced  the  concept  of  "Phantom  Antennas" 
and  this  system  is  shown  in  Figure  2 „  This  is  the  system  which  is 
presently  being  used  to  measure  elevation  angles  on  the  University  of 
Illinois  Wullenweber  Antenna  Array, 

Measurement  of  elevation  angle  by  the  Phantom  Antenna  concept 
is  based  upon  the  measurement  of  the  differential  phase  between  a  Front 
Phantom  and  one  or  more  Rear  Phantoms.  The  Front  Phantom  is  formed 
by  selecting  two  adjacent  antennas  on  the  array  such  that  the  direction 
of  arrival  of  the  desired  incoming  signal  passes  between  these  two 
antennas.  The  output  from  each  of  these  antennas  is  combined  by  an 
RF  combiner  to  create  the  Front  Phantom.  By  the  phasor  addition  of 
voltages  from  other  symmetrically  disposed  pairs  of  antennas,  the 
Rear  Phantoms  are  formed  and  they  are  located  directly  behind  the  Front 
Phantom.  For  elevation  angle  measurements  so  far  conducted  on  the 
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University  of  Illinois  Wullenweber,  at  feast  two  Rear  Phantoms  have  been 
in  use  and  sometimes  three.  To  avoid  ambiguity  problems,  the  horizontal 
spacing  between  the  Front  and  Rear  Phantoms,  d,  is  restricted  to  a  dis¬ 
tance  of  one-half  wavelength  or  less  for  any  given  operating  frequency. 

For  this  configuration,  the  elevation  angle,  0,  can  be  determined 
by  solving  the  equation: 

2nd 

,j,  cos  o  cos  a 

A  (5) 

where  d  =  horizontal  spacing  between  respective  Front  and  Rear  Phantoms 
\  =  wavelength  at  the  operating  frequency 

a  =  difference  between  the  angle  of  arrival  of  signal  and  the 
azimuthal  alignment  of  Phantom  Antennas 

U  =  differential  phase  angle  measured  between  Phantom  Antennas 

c  =  elevation  angle  measured  from  the  ground  plane 

To  determine  the  elevation  angle  0,  the  differential  phase  between 

the  Front  and  Rear  Phantoms,  1,  is  measured  by  a  phasemeter.  The  output 

from  the  phasemeter,  an  analog  voltage  which  is  proportional  to  this 

phase  difference,  is  recorded  by  the  HSDS,  By  means  of  a  G-20  computer 

program,  the  value  of  for  use  in  equation  (5)  is  calculated,  d  and  \ 

are  already  known  for  any  given  operating  frequency,  and  the  elevation 

angle  9  can  be  determined.  For  the  studies  so  far  conducted^,  the 

value  of  a  in  the  above  equation  has  been  zero  or  nearly  zero  degrees; 

therefore,  this  quantity  has  not  been  used  in  any  of  the  elevation 

angle  calculations  which  so  far  have  been  made. 

IV.  CONCLUSIONS 

By  means  of  these  two  methods  which  involve  measurement  techniques 
which  are  completely  independent  of  each  other,  it  is  hoped  that  some 
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measure  of  elevation  angle  may  be  obtained.  it  is  further  hoped  LhaL 
these  two  methods  of  measurement  will  yield  elevation  angle  values 
which  are  somewhat  correlated.  This  has  not  necessarily  been  Lrue 
£or  elevation  angles  obtained  by  use  of  different  Phantom  Antenna 
pairs  and  some  adjustment  of  this  method  appears  necessary. 

Based  upon  the  daLa  so  far  obtained,  it  would  appear  LhaL  for 
elevation  angles  above  15  degrees,  either  or  both  of  these  methods 
wilL  give  some  indication  of  elevation  angle.  However,  for  elevation 
r.ngies  from  15  degrees  down  to  zero  degrees,  both  methods  of  measure¬ 
ment  become  increasingly  insensitive  and  the  results  of  measurements 

in  this  range  are  suspect. 
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OBLIQUE  INCIDENCE  LONOSONDE  AMl’LITUDE  DATA 

Engineering  changes  to  the  University  of  Illinois  ionosonde  equip¬ 
ment  include  the  addition  of  the  "A"  scan  data  and  automatic  gain  con¬ 
trol  (AGC) . 

The  "A"  scan  data  are  recorded  directly  on  the  ionogram  at  the  fre¬ 
quency  of  interest;  the  "A"  scan  is  deflected  at  a  45°  angle.  The  single 
line  scope  trace  is  then  shifted  to  form  a  "gap";  the  "gap"  does  not 
represent  missing  frequencies.  The  length  of  the  line  of  the  "A"  scan 
trace  is  an  indication  of  the  relative  signal  strength  Erom  each  mode 
of  transmission.  The  45’  deflection  (as  shown  in  Figures  1  and  2) 
increases  the  resolution  of  the  "A"  scan  trace  length. 

The  ionosonde  AGC  adjusts  for  the  difference  in  day-night  signal- 
to-noise  ratio.  The  AGC  samples  the  background  noise  level  and  adjusts 
the  gain  to  a  preset  optimum  noise  level.  A  complete  report  on  the  AGC 
system  will  be  available  in  the  near  future.^"  Figures  3  and  4  are  exam¬ 
ples  of  ionograms  obtained  without  (Fig.  3)  and  with  (Fig.  4)  AGC- 


"An  AGC  System  for  a  Step  Frequency  Ionosonde  Receiver,"  R.  Baumgartner, 
E.  W.  Ernst,  F.  0.  Fahlsing,  W.  W.  Wood,  RRL  Publication  No.  356,  Tech¬ 
nical  Note  No.  10,  Radiolocation  Research  Laboratory,  Dept,  of  Electrical 
Engineering,  University  of  Illinois,  Urbana,  Illinois. 


Figure  2. 


Figure  1. 


,  . .  .  , 

•  t  *  ^  . L/ 

«  mi'll. i!ii|«lilllltll||l‘MMRii|,,,,,,>l|l<>ll>  ]  J 

■  -  i  .  i  •  iJ-  '  » 

i«imuuiiMU»uu*iiiu^ni-'llpummu4<un  j 


liMUf lUUJUJUtJI  I  •*** 
. . 


,'VV 


PRECEDING  PAGE  BLANK  -  NOT  FILMED  ^ 


2'Jl 

GliKAl.D  A.  SV i  i'll 

.v,:i.i2i£Mi,n:  w>;’  A'.hv:  tinj 

L'livcrsi  :y  J'f 

W A VK FRONT  TESTING  AND  ITS  APPLICATIONS 

the  triple  interferometer  lias  been  used  tor  several  years  in 
radioiucat ion  work.  its  ability  to  resolve  angle  of  arrival  in 
elevation'  and  azimuth  is  a  great  asset.  However,  due  lo  the  naLure 
of  the  interferometer,  it  looks  in  all  directions.  This  "open  window" 
allows  signals  from  all  directions.  One  problem  this  brings  is  the  inabilit 
to  shut  out  unwanted  signals.  These  unwanted  signals  can  be  placed 
in  two  categories.  The  first  occurs  when  two  transmitters  are  at  tiie 
same  frequency  and  at  different  locations.  This  does  not  pose  too 
much  of  a  problem  because  it  is  not  likely  the  transmitters  will  be 
close  enough  in  frequency  (i.e.  within  one-half  cycle  per  second)  to 
prevent  accurate  determination  of  direction  of  arrival, 

The  unwanted  signal  that  poses  the  greatest  problem  is  the 
multipath  propagated  signal.  This  signal,  transmitted  from  one  source, 
travels  via  two  paths  to  the  receiver.  This  then  appears  as  two 
signals,  at  slightly  different  frequencies  at  the  receiver.  The 
difference  may.be  the  order  of  one  cycle  per  minute  to  one  cycle  per 
second  The  closer  the  two  frequencies,  the  more  difficult  it  is  to 
resolve  the  incidence  and  azimuth  of  either. 

The  multipath  propagation,  referred  to  as  modes,  modifies  the 
interferometer  output  in  such  a  way  as  to  render  the  angle  of  arrival 
measurements  useless  unless  means  for  resolving  the  modes  are  used. 


Assuming  two  modes  are  present,  it  has  been  proposed  that  one 
might  find  the  correct  angles  of  arrival  if  there  were  a  way  of  knowing 


JL  what  instant  in  Lime  llie  second  path  (or  Llie  iirsL)  were  nul  present 
it  is  known  i_hat  Lor  a  single  mode  wave,  Llie  niagniLude  ul  Llie  signal 
presenL  at  each  antenna  is  equal.  This  suggesis  LliaL  each  anLenna 
should  be  checked  and  ii  Llie  magnitude  ol  Llie  signals  is  indeed  equal, 
(implying  a  single  mode  signal)  then  a  measurement  should  be  made. 

une  method  of  accompl isliing  this  is  to  determine  the  magnitude 
of  the  signal  at  each  anLenna,  compare  these  signals  (three  in  the 
case,  of  the  triple  interferometer)  and  if  they  are  equal,  display  the 
outputs  so  that  the  answer  might  be  computed. 

A  preliminary  investigation  at  the  University  of  Illinois  has 
shown  that  logical  circuits  will  perform  the  necessary . functions  required. 

Some  problems  have  been  encountered;  however,  these  are  unique 
for  our  particular  installation.  First,  we  use  a  time -shared  ,twin- 
chanael  receiver  to  sample  the  phase  difference  and  amplitude  of  three 
pairs  of  antennas.  This  would  require  a  memory  device  to  store  the 
magnitude  information  and,after  the  complete  cycle,  make  the  magnitude 
comparison.  This  results  in  information  from  the  previous  cycle 
controlling  the  display  for  the  present  cycle.  With  a  triple -channel 
receiver,  the  memory  would  be  unnecessary. 

Second,  the  storage  oscilloscope  that  we  use  for  the  display 
requires  AC  unblanking;  therefore,  an  additional  circuit  must  be 
added  to  furnish  the  unblanking  signal  . 

Wavefront  testing,  in  our  case,  has  been  limited  to  testing  the 
antenna  pair  being  sampled  and, if  the  magnitudes  are  equal,  displaying 
the  data.  This  method  displays  the  results  of  only  one  antenna  pair 
or  two  pairs,  or  all  three.  Should  either  of  the  first  two  occurrences 
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i'Ci'uI  quickly,  one  would  bo  let!  to  believe  that  a  valid  solution 
was  present  when  it  was  not.  iL  has  also  been  noted,  in  our  ease, 
that  when  Lite  two  modes  are  in- phase  or  1 80  degrees  out -ol -phase .  the 
magnitudes  ot  antenna  pairs  can  be  equal,  which  leads  again  to  a 
false  belief  that  a  valid  solution  is  present.  As  a  result,  lour 
solutions  for  the  Lwo  mode  eases  are  possible. 

At  present  we  have  been  able  to  record  equal  magnitude  signals 
only  when  the  two  signals  are  in-phase  or  ISO  degrees  out-of-phase 
This  is  not  to  say  that  the  waves  are  never  "single  mode"  during 
interference;  it  has  been  noted  to  happen  but  not  recorded.  it  means 
that  at  the  times  pictures  were  taken,  the  double  mode  propagation 
did  not  disappear. 

Only  a  limited  amount  of  wavefront  testing  has  been  performed, 
and  this  on  a  highly  restricted  basis  which  prohibits  any  firm  conclu¬ 
sions  from  being  reached,  Wavefront  testing  shows  some  promise  and 
also  some  pitfalls.  It  may  be  necessary  to  wait  2  minutes  for  5  seconds 
of  display,  for  there  is  no  assurance  that  one  of  the  two  modes  will 
disappear  completely  or  rapidly.  Also,  three  or  more  modes  may  be 
present.  One  must  then  wait  for  all  but  one  to  disappear, 

There  is  some  thought  that  the  "false  solutions,"  when  interfer¬ 
ence  is  in-phase  or  180  degrees  out-of-phase,  will  be  removed  when 
the  decision  to  display  data  is  based  on  all  three  antenna  magnitudes 
rather  than  two.  We  do  know  that  errors  caused  by  displaying  solutions 
when  one  or  two  pairs  are  displayed  will  be  eliminated. 

There  may  be  value  in  the  fact  that  the  amplitudes  are  equal 
when  we  have  interference  at  the  in-phase  and  out-of-phase  conditions 
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because  this  gives  die  major  diagonals  oi  Lhe  polygon  described  by 
the  two  modes.  This  information  eouid  cause  Lhe  data  acquisition 
system  to  record  these  diagonals  and  then  Lhe  computer  could  recon¬ 
struct  the  polygon  to  extract  useful  information. 
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